Made available in DSpace on 2014-06-11T19:23:38Z (GMT). No. of bitstreams: 0 Previous issue date: 2008-09-16Bitstream added on 2014-06-13T19:50:37Z : No. of bitstreams: 1 bornschlegell_as_me_ilha.pdf: 1868861 bytes, checksum: b2f9910d0ae7f6f842769165c54b95de (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Durante as últimas décadas, grandes esforços têm sido dirigidos ao desenvolvimento de técnicas numéricas capazes de tratar escoamentos envolvendo geometrias complexas e/ou móveis. Neste contexto, o chamado método de fronteira imersa, que permite a resolução desta classe de problema usando malha cartesiana fixa independentemente do tipo de geometria tratada, tem recebido grande atenção por parte da comunidade internacional de dinâmica dos fluidos computacional. No presente trabalho, o método de fronteira imersa, associado ao modelo físico virtual, foi implementado em um código computacional capaz de simular escoamentos isotérmicos, incompressíveis e bidimensionais, utilizando o método dos volumes finitos. Escrito em linguagem Fortran 90, este programa emprega o método SIMPLEC para o acoplamento pressão-velocidade, oferecendo ao usuário diferentes opções para o tratamento dos termos advectivos – Upwind, Diferenças Centrais, Power-Law e QUICK. Embora a linguagem Fortran 90 não possibilite a programação orientada a objetos, procurou-se, na execução deste trabalho, aplicar alguns dos conceitos envolvidos neste tipo de programação, tais como polimorfismo, herança e encapsulamento. Testes exaustivos foram executados para avaliar as funcionalidades (funções e sub-rotinas) introduzidas no código computacional e a performance global do programa. Neste contexto, dois problemas clássicos da literatura foram utilizados como benchmark na validação do programa: o escoamento em torno de um cilindro estacionário e o escoamento ao redor de um cilindro em rotação. Em ambos os casos, os resultados obtidos se mostraram fisicamente consistentes, apresentando uma concordância satisfatória com dados da literatura / In the last decades, considerable efforts have been directed to the development of numerical techniques for solving fluid flow past complex and/or moving geometries. In this context, the so-called immersed boundary method, which allows treating this class of problems by using Cartesian grids independently of the considered geometry, has received large attention from international community of computational fluid dynamics. In this work, the immersed boundary method associated to the physical virtual model has been computationally implemented into a two-dimensional finite volume for simulating incompressible isothermal flows. The SIMPLEC algorithm has been employed to the velocity-pressure coupling, while different schemes – second order upwind, central difference, power-law, and QUICK – were available to discretize the advective terms. The code has been written in FORTRAN 90. Although such a compiler does not provide support to object oriented programming, the concepts of polymorphism, inheritance and encapsulation have been applied during implementation of the virtual physical model. Exhaustive tests have been performed in order to evaluate not only the new functions and subroutines added to the main program, but also the whole program itself. For this, the flows past stationary and rotating circular cylinder have been considered as benchmarks. In both cases, the obtained results are physically consistent, presenting a satisfactory agreement with the available data in the literature

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Systemic sclerosis (SS) is a chronic and multisystemic disease characterized by fibrosis and vascular abnormalities in the skin and in visceral organs. Pulmonary complications are the most common cause of death and morbidity associated to the disease. The pulmonary function tests are abnormal in 67-95% of patients with SS. Some studies have demonstrated the potential of the forced oscillation technique (FOT) in the detection of respiratory mechanics modifications in several respiratory diseases, but there are no studies in the literature using this technique in patients with systemic sclerosis. Thus, the goals of the present study were: (1) to analyze changes in resistive and reactive properties of the respiratory system in patients with SS and (2) to investigate the diagnostic potential of this technique. This research consists of a controlled observational study. In the first part of the study, considering spirometry, 82 individuals were analyzed, 30 control subjects and 52 patients with diagnostic of SS, 22 with normal spirometry and 30 with restrictive alteration to spirometry. In a second study, based on evaluation of lung volumes, 50 individuals were evaluated, 21 control subjects and 29 with systemic sclerosis, 7 with normal result on volumes exam and 22 with restrictive ventilatory disturbance in this exam. All subjects performed FOT exam for the analysis of resistive and reactive properties of the respiratory system. In subsequent examinations, the individuals were submitted to spirometry and/or volumes exam. Considering the groups divided from the spirometry, the resistive and reactive parameters and the impedance module of the respiratory system at 4Hz (|Z4Hz|) showed significantly alterations (KW- ANOVA/ANOVA, p<0.004). In the comparison between groups, the restrictive group showed the highest changes in all study variables. Considering the division by lung volumes exam, the elastic parameters, |Z4Hz| and the slope of the resistance (S) showed significant changes (KW-ANOVA, p <0.006), while the others resistive variables (total and mean resistance, R0 and Rm, respectively) showed a tendency to increase. In this analysis, the restrictive group also showed greater impairment in oscillometric and volumetric parameters. The SS resulted in changes in resistive and reactive properties of the respiratory system. Analyzing the diagnosis power of FOT, with spirometry as reference, when comparing the normal exam group with restrictive group, the mean reactance (Xm) showed useful clinical value with an area under the curve (AUC) of 0.81. In the comparison between the control and restrictive groups, R0, dynamic compliance (Cdyn,rs) and |Z4Hz| showed the following AUC values, 0.81, 0.95 and 0.91, respectively, representing high accuracy for the last two variables. Using volumes exam as reference, comparing control and restrictive groups, the Xm and resonance frequency (fr) were noteworthy with values of 0.82 and 0.81 of AUC, respectively. Those results are consistent with physiopathological changes related to SS, which may include diffuse fibrosis that affects the lung, the bronchial wall and the alveolar septum with changes in parenchimal architecture, confirming the potential of FOT in the assessment of modifications in respiratory mechanics related to SS. / A Esclerose sistêmica (ES) é uma doença crônica, multissistêmica caracterizada por fibrose e anormalidades vasculares na pele e em órgãos viscerais. As complicações pulmonares são as principais causas de morte e morbidades associadas à doença. Os testes de função pulmonar apresentam-se anormais em 67-95% dos pacientes com ES. Alguns estudos têm demonstrado o elevado potencial da Técnica de Oscilações Forçadas (FOT) na detecção das alterações da mecânica respiratória de diversas doenças respiratórias, no entanto não existem trabalhos na literatura que utilizem esta técnica em pacientes com esclerose sistêmica. Desta forma, os objetivos desta pesquisa foram (1) analisar as alterações das propriedades resistivas e reativas do sistema respiratório em portadores de ES e (2) investigar o potencial de uso diagnóstico desta técnica. Trata-se de um estudo observacional controlado. Na primeira parte do estudo, realizada com base no exame espirométrico, foram analisados 82 indivíduos, 30 indivíduos controles e 52 pacientes com diagnóstico de ES, sendo 22 normais ao exame espirométrico e 30 com distúrbio ventilatório restritivo à espirometria. Em um segundo estudo, baseado na avaliação de volumes pulmonares, 50 indivíduos foram avaliados, 21 indivíduos controles e 29 com ES, sendo 7 normais ao exame de volumes e 22 com distúrbio ventilatório restritivo detectado por esse exame. Todos os indivíduos realizaram exames da FOT para análise das propriedades resistivas e reativas do sistema respiratório. Posteriormente aos exames da FOT, os indivíduos foram submetidos à espirometria e/ou exame de volumes. Considerando os grupos divididos a partir da espirometria, os parâmetros resistivos, reativos e o módulo da impedância do sistema respiratório em 4Hz (|Z4Hz|) modificaram-se significativamente (KW- ANOVA/ANOVA, p<0,004). Na comparação entre grupos, o grupo restritivo apresentou as maiores alterações em todas as variáveis estudadas. Considerando o exame de volumes pulmonares, os parâmetros elásticos, |Z4Hz| e a inclinação da curva de resistência (S) apresentaram alterações significativas (KW-ANOVA, p<0,006), já as demais variáveis resistivas (resistências total e média, R0 e Rm, respectivamente) demonstraram uma tendência à aumento. Nessa análise, o grupo restritivo também apresentou maior prejuízo nos parâmetros oscilométricos e volumétricos. A ES resultou em alterações nas propriedades resistivas e reativas do sistema respiratório. Na análise do poder diagnóstico da FOT, sendo a espirometria referência, na comparação entre o grupo normal ao exame e o restritivo, a reatância média (Xm) apresentou valor útil para uso clínico com valor de área sob a curva (AUC) de 0,81. Já na comparação entre os grupos controle e restritivo, R0, complacência dinâmica (Cdin,sr) e |Z4Hz| apresentaram os seguintes valores de AUC, 0,81, 0,95 e 0,91, respectivamente, representando alta acurácia para as últimas duas variáveis. Com o exame de volumes como referência, comparando os grupos controle e restritivo, a Xm e a freqüência de ressonância (fr) destacaram-se com valores de 0,82 e 0,81 de AUC, respectivamente. Esses resultados são coerentes com as alterações fisiopatológicas relacionadas à esclerose sistêmica, que podem incluir fibrose difusa envolvendo o interstício pulmonar, a parede brônquica e o septo alveolar com alteração da sua arquitetura, confirmando o potencial da FOT na avaliação das modificações na mecânica respiratória relacionadas à ES.

No presente trabalho avalia-se o comportamento da difusão e da convecção natural em sistemas térmicos que possuem armazenadores de calor, como os que envolvem coletores solares ligados a reservatórios de água quente. Através da solução das equações diferenciais que regem o problema, obtém-se os campos de velocidade e temperatura no interior destes armazenadores. Este último, especialmente, é importante para a determinação da performance térmica de todo o sistema. A intensidade da estratificação térmica depende do fluxo de entrada e saída de água no reservatório, da sua razão de aspecto, da posição da resistência elétrica de apoio, se ela existir, da espessura e difusividade térmica da parede e da espessura e material do isolamento térmico. A perda de calor para o meio ambiente causa um fluxo de líquido descendente, junto às paredes do reservatório, incorrendo em camadas limites hidrodinânlica e térmica. Além do papel destas camadas limites, também investiga-se a importância relativa da difusão de calor através da parede metálica do armazenador, na sua direção axial (altura). Não são investigados neste trabalho os efeitos de entrada e saída do líquido do armazenador, e também as correntes convectivas originadas do suprimento de calor por resistência elétrica. As equações da energia e do movimento são resolvidas de forma acoplada, em coordenadas cilíndricas axissimétricas, em regime transiente, empregando o Método dos Volumes Finitos. / In the present work the diffusion behavior as well as the natural convection in thermal systems that possess heat reservoirs, like the ones that involve solar collectors connected to hot water reservoirs, have been evaluated. Through the solution of the differential equations that describe the problem, the velocity fields and the internal temperature of these tanks have been obtained. The last one is specially important for the determination of the thermal performance of the whole system. Intensity of the thermal stratification depends on the entrance and exit water flow in the reservoir, on the tanks aspect ratio, on the position support of the electrical resistance, if it exists, on the thickness and thermal diffusion of the wall and on the material of thermal insulation. The loss of heat for the environment brings about a descending liquid flow, next to the walls of the reservoir, ending up by creating hidrodynamic and thermal boundary layers. Besides the role of these boundary layers, the relative importance o f the heat diffusion through the metallic wall of the reservoir in its axial direction (height), also has been investigated. The effects of the entrance and exit of liquid in the reservoir, as well as the convective currents originated from the heat supplied by electric resistance have not been investigated. The energy and moment equations are solved in a coupled way, in axissimetrical cilyndrical coordinates, in transient regime, employing a Finite Volume Method.

Orientador: Eugenio Spano Rosa / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-28T23:56:11Z (GMT). No. of bitstreams: 1 MelgarejoMorales_RigobertoEleazar_D.pdf: 8900310 bytes, checksum: 9f57b0d7b80e8430b0d0c8941dda51e1 (MD5) Previous issue date: 2000 / Resumo: Neste trabalho é realizado uma simulação numérica do escoamento com superficie livre completamente desenvolvido em um canal helicoidal de seção retangular, em regime laminar e turbulento. As equações de conservação foram escritas para um sistema de coordenadas ortogonais locais. O sistema de equações foi discretizado e resolvido pelo método de volumes fmitos com o esquema ln'brido de interpolação, implementados no PHOENICS-CFD. A interface livre é calculada utilizando o algoritmo HOL (Height-Of-Liquid). As simulações numéricas foram realizadas para fluidos com viscosidade de 1 a 100cP e Reynolds de 352 a 45500. A turbulência é modelada com o modelo de duas equações 1(-E para baixo números de Reynolds com uma correção devido ao efeito da curvatura. As simulações mostram que o método utilizado, quando comparado com dados experimentais, calcula satisfatoriamente a localização da interface. Apresenta-se os campos da velocidade média, do escoamento secundário, da energia cinética turbulenta e da dissipação e produção da energia cinética turbulenta. A partir dos resultados numéricos calcula-se o fator de atrito local e médio / Abstract: A numerical simulation of the fully developed flow in an open helical channel with rectangular cross section is accomplished, for the laminar and turbulent regime. The conservation equations are written in a local orthogonal coordinate system. The conservation equations are solved by the finite volume method with the hybrid scheme of interpolation, implemented in PHOENICS-CFD. The ftee surface is calculated using the HOL algorithm (Height-Of-Liquid). The numerical simulations were accomplished for fluids with viscosity varying ftom 1 to 100cP and Reynolds ftom 352 to 45500. The turbulence is modeled by two equations k-E model for low Reynolds numbers with a correction due to the effect of the curvature. The numerical results for the tree surface position agrees within :t4 mm against the experimental data. The averaged axial velocity and the secondary flow fields are presented, as well as the turbulent kinetic energyand the dissipation of the turbulent kinetic energy. Finally the local and the averaged ftiction factors are presented. / Doutorado / Doutor em Engenharia Mecânica

Orientador: Jose Ricardo Figueiredo / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-03T18:04:25Z (GMT). No. of bitstreams: 1 Vilela_CarlosAlbertodeAlmeida_D.pdf: 6593008 bytes, checksum: 706d78c356008d4b1f84aa2b735170e6 (MD5) Previous issue date: 2001 / Resumo: O esquema de discretização UNIFAES, é aqui avaliado e comparado com alguns dos mais tradicionais e utilizados esquemas de discretização para o método dos volumes finitos da atualidade. Diversos casos testes lineares do problema de transporte convectivo-difusivo em meio fluido e das equações de Navier-Stokes foram utilizados para comparação. Os esquemas utilizados para comparação foram o Central, Exponencial, LOADS e QUICK em duas versões diferentes de formulação. Os casos testes realizados para avaliação do comportamento numérico realizados em coordenadas generalizadas com malhas regulares e irregulares foram: Transporte convectivo-difusivo em, a) escoamento paralelo b) escoamento rotacional de fluido como sólido, c) escoamento em canto, d) escoamento em canto com perfil de entrada parabólico, escoamento em cavidade fechada inclinada e escoamento em canal divergente gradual. Os resultados obtidos nos testes demonstraram que o esquema UNIFAES, apresentou comportamento bastante promissor comparado aos comportamentos dos demais esquemas. Na grande maioria dos testes os resultados obtidos pelo esquema UNIFAES foram os melhores em questão de acuidade, estando próximo aos obtidos pelo esquema QUICK, que foi o segundo esquema que melhores resultados apresentou / Abstract: The UNIFAES discretization scheme is evaluate and compare with some classical and most used schemes for finite volume method in present time. A lot of linear convection-diffusion problems in fluid flow and the solution of Navier-Stokes equations were used for comparison. The compared used schemes were Central, Exponential, LOADS, and QUICK with two different formulations. The chosen test cases for comparisons, using regular and irregular grid arrangements in generalized coordinates formulations were: a) Convective-diffusive transport in parallel fluid flow, in rotational fluid flow, in fluid flow near a wall, in fluid flow near a wall with parabolic inlet profile of the transported property , b) Inclined side walls lid driven cavity, and c) Gradually expanding channel. The obtained results shown that UNIFAES scheme has a very promise behavior compared to the other schemes. In most studied cases, UNIFAES scheme demonstrated better accurate results, elose to QUICK, the second better one / Doutorado / Doutor em Engenharia Mecânica

The resource modulation hypothesis suggests that the influence of genes on cognitive functioning increases with age. The KIBRA single nucleotide polymorphism rs17070145, associated with episodic memory and working memory, has been suggested to follow such a pattern, but few studies have tested this assertion directly. The present study investigated the relationship between KIBRA alleles (T carriers vs. CC homozygotes), cognitive performance, and brain volumes in three groups of cognitively healthy adults-middle aged (ages 52-64, n = 38), young old (ages 65-72, n = 45), and older old (ages 73-92, n = 62)-who were carefully matched on potentially confounding variables including apolipoprotein epsilon 4 status and hypertension. Consistent with our prediction, T carriers maintained verbal memory performance with increasing age while CC homozygotes declined. Voxel-based morphometric analysis of magnetic resonance images showed an advantage for T carriers in frontal white matter volume that increased with age. Focusing on the older old group, this advantage for T carriers was also evident in left lingual gyrus gray matter and several additional frontal white matter regions. Contrary to expectations, neither KIBRA nor the interaction between KIBRA and age predicted hippocampal volumes. None of the brain regions investigated showed a CC homozygote advantage. Taken together, these data suggest that KIBRA results in decreased verbal memory performance and lower brain volumes in CC homozygotes compared to T carriers, particularly among the oldest old, consistent with the resource modulation hypothesis.

Un des ingrédients principaux pour la simulation numérique des écoulements des fluides (hydrocarbures, gaz naturel) en milieux poreux est la discrétisation des termes elliptiques anisotropes et hétérogènes. Dans l’industrie pétrolière, la nécessité d’améliorer la précision des simulations dans les régions proches des puits a suscité l’utilisation de maillages non structurés généraux et des tenseurs de perméabilité pleins. Notre défi a été de trouver des discrétisations consistantes et robustes des termes elliptiques anisotropes, hétérogènes sur maillages généraux. Notre recherche s’est focalisée sur des méthodes volumes finis qui soient consistantes, stables sur maillages polyédriques généraux, robustes par rapport à l’anisotropie et l’hétérogénéité du tenseur de perméabilité, qui mènent à des systèmes linéaires bien conditionnés pour lesquels des stratégies de préconditionnement optimales peuvent être conçues et qui donnent un stencil compact pour réduire les communications dans les implémentations parallèles. Pour répondre à cette recherche, nous avons proposé plusieurs schémas tels que MPFA O généralisé, G scheme, CG method, VFSYM, DIOPTRE. Nous avons prouvé également que toutes ces méthodes convergent sous des hypothèses adéquates à la fois sur le tenseur de perméabilité et le maillage / One of the key ingredients for the numerical simulation of Darcy flow in heterogeneous porous media is the discretization of anisotropic heterogeneous elliptic terms. In the oil industry, the need to improve accuracy in near wellbore regions has prompted the use of general unstructured meshes and full permeability tensors. Our effort has therefore been devoted to find consistent and robust finite volume discretizations of anisotropic, heterogeneous elliptic terms on general meshes. Our research was focused on finite volume methods which are consistent and coercive on general polyhedral meshes as well as robust with respect to the anisotropy and heterogeneity of the permeability tensor ; yield well-conditioned linear systems for which optimal preconditioning strategies can be devised ; have a narrow stencil to reduce the communications in parallel implementations. To answer to this search, we have proposed several scemes such that generalized MPFA O, G scheme, CG method, VFSYM, DIOPTRE. We proved also the convergence of all these methods under suitable assumptions on both the permeability tensor and the mesh

Orientador: Rubens Maciel Filho / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-03T15:03:24Z (GMT). No. of bitstreams: 1 Barreira_MardonnyNazareno_M.pdf: 5649911 bytes, checksum: 74e112ba543ced12ed31508868adfdde (MD5) Previous issue date: 2003 / Mestrado

Orientador: Eugênio Spanó Rosa, Marcelo Moreira Ganzarolli / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-03T15:19:17Z (GMT). No. of bitstreams: 1 Saito_MarceloBatista_M.pdf: 6699986 bytes, checksum: 845433a3bcb3aa43a33251e4350ebf40 (MD5) Previous issue date: 2002 / Resumo: O presente trabalho tem por objetivo determinar numericamente a queda de pressão e o fluxo de calor de um trocador de calor ar-ar utilizado para resfriamento de um ambiente selado que não pode ser contaminado pelo ar externo. Esta análise tem por finalidade otimizar a geometria que maximiza o calor transferido para uma determinada queda de pressão. O trocador é formado por dutos paralelos de seções retangulares com escoamento em contra-corrente no regime turbulento. As duas correntes são separadas por paredes sólidas e lisas. A análise paramétrica é efetuada resolvendo-se equações de conservação da massa, quantidade de movimento e energia, uti1izando o modelo de turbulência K-E padrão e considerando propriedades variáveis. A condução de calor nas paredes sólidas também foi incluída nas simulações. A solução numérica tridimensional é obtida através da discretização das equações de conservação em volumes finitos, sendo empregado um pacote numérico comercial (pHOENICS 3.3). As simulações são realizadas para número de Reynolds entre 8000 e 160000. Os resultados obtidos permitem avaliar o desempenho do trocador de calor para valores especificados da queda de pressão e determinar a geometria que maximiza o calor transferido / Abstract: The pressure drop and heat transfer flux on an air-air heat exchanger employed as a cooling device for a closed cavity with heat dissipation are determined numerically. The analysis is further extended to optimize the cross section dimensions to render the maximum heat transfer rate at a given pressure drop. The heat exchanger is of rectangular cross section ducts placed side by side in counter flow turbulent regime. The duct walls are of smooth and flat aluminum. The conservation equations of mass, momentum and energy are solved with the K-B turbulence model using the finite volume code embodied on the PHOENICS 3.3. The heat conduction along the solid walls is considered employing conjugate heat transfer treatment. The simulations are performed for Reynolds number between 8000 and 160000. The pressure drop and heat transfer rate are determined and expressed in terms of dimensionless parameters. The analysis is further advanced determining cross section aspect ratio which gives the maximum heat transfer rate for a given pressure dro / Mestrado / Termica e Fluidos / Mestre em Engenharia Mecânica

Orientadores: Maria Regina Wolf Maciel, Henry França Meier, Antonio Andre Chivanga Barros / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-03T15:31:22Z (GMT). No. of bitstreams: 1 Noriler_Dirceu_M.pdf: 4837242 bytes, checksum: ed71f25bb9219c76441c2042ff4f6a07 (MD5) Previous issue date: 2003 / Resumo: Os modelos para representação de colunas de destilação atuais baseiam-se no conceito de estágios de equilíbrio e de não equilíbrio, ambos considerando o equilíbrio hidrodinâmico em uma abordagem macroscópica. No entanto, em meados da década de 90, iniciou-se a aplicação de modelos microscópicos para subsidiar e incrementar a análise e otimização dos processos de refino, dando-se início à era dos códigos comerciais de Fluidodinâmica Computacional. O principal objetivo deste trabalho é desenvolver um modelo microscópico para descrever a fluidodinâmica do escoamento líquido-vapor em um prato perfurado de uma coluna de destilação, com base na conservação de quantidade de movimento sob condições de turbulência e aplicar e implementar o modelo matemático no código comercial de fluido dinâmica computacional CFX, da ABA Technology, para a realização de um estudo de experimentação numérica de modelos e métodos. Foram propostos três modelos distintos para representação do escoamento em um prato perfurado de destilação: um modelo bidimensional e monofásico, 2-DSM (Liu et ai. (1999)) , um modelo tridimensional bifásico e homogêneo, 3-DHomM, e por fim um modelo tridimensional bifásico e heterogêneo, 3-DHetM. Desta forma, os três modelos foram implementados no código comercial CFX 4.4 da ABA technology, o qual utiliza o método dos volumes finitos com variáveis co-localizadas resolvendo o sistema em coordenadas generalizadas. As malhas numéricas foram construídas por intermédio da técnica dos multiblocos, utilizando o pré-processador "BUILD" presente no código. Os principais resultados mostram que o modelo 3-DTHomM representa, qualitativamente, de forma satisfatória a hidrodinâmica de um prato perfurado de destilação. A metodologia proposta neste trabalho se mostrou adequada e apta a reproduzir tal comportamento fluidodinâmico / Abstract: The current models for predicting distillation column behavior are based on equilibrium and non-equilibrium stage concepts, but both use the hydro dynamical equilibrium in a macroscopic point of view. The objectives of this work are: developing the microscopic model to describe the gas-liquid flow on a distillation sieve tray, based on momentum equations under turbulence conditions, inter-phase momentum transfer and other modeling features; and applying the mathematical modeling in a commercial CFD code, for numerical experiments in order to verify the software and to validate the models. Three models were proposed: the first model named Two Dimension Single Model (2-DSM model), with two dimension and single phase assumptions; the second one, named Three Dimension Homogeneous Model (3-DHomM model), with two phases and three dimension with homogeneous flow, i. e., the liquid and the gas velocities are the same; and the last model is a model with two phases and three dimension assumptions, named Three Dimension Heterogeneous Model (3-DHetM model). These three models were applied in CFX 4.4 CFD commercial code from ABA Technology, by using the finite-volume methods with co-located variable in a generalized coordinate system, to obtain numerical solutions. The pre-processor BUILD were used to build up the numerical grid with the multi-block technique. The main results show that the 3-DHomM model, despite of its simplicity when compared with 3-DhetM, describes the fluid dynamics on a sieve distillation tray satisfactorily. The methodology proposed in this work also showed to be adequate, and with capability to describe the gas-liquid flow on a sieve distillation tray / Mestrado / Engenharia de Processos / Mestre em Engenharia Química