Neste trabalho são apresentadas soluções numéricas para o escoamento laminar e turbulento em um canal aletado com material poroso, para um escoamento de fluido monofásico em regime permanente, no qual a região porosa é modelada como sendo homogênea e isotrópica. A condição de periodicidade espacial ao longo do canal é empregada. As equações do movimento e continuidade de massa são integradas em um volume elementar representativo acarretando em um único conjunto de equações governantes válido para todo o domínio computacional. Estas equações são discretizadas pelo método de volume de controle e o sistema de equações resultante é resolvido pelo método SIP [Strongly Implicity Procedur

In this study the influence of physical properties on heat transfer between the fluid and solid phases for laminar and turbulent flow in a square cavity filled with porous material, heated by one side and cooled by the opposite one is investigated. In order to simulate flow and heat transfer between the phases models, one and two-energy-equation models are used. The transport equations are discretized using the control volume method and the system of algebraic equations is relaxed via the SIMPLE algorithm. Validations were perfomed first for the case of clean cavity, using Laminar and High Reynolds Turbulent model, and the results are in agreement with the existing literature. In the case of porous cavities, simulation runs were performed using the 1EEM and 2EEM. Overall, this study showed that as porosity increases, the value of the average Nusselt at hot wall, for the same Ram, decreases. It is also shown that the average Nusselt decreases, from a critical Rayleigh , as the ratio between the thermal conductivities of the phases increases. Comparisons were made between both energy models and 2EEM shown smaller results for average Nusselt as ks/kf increases. Also found that the critical Rayleigh depends directly on the ratio ks/kf, being larger as it also grows. Still object of this work, it has been discovered that as the porosity decreases and the ratio between thermal conductivities increases, the turbulence within cavity increases as well. In general, the results and analysis obtained by this work can be used in real engineering situations where the porous cavities may be used.

As velocidades de chamas laminares pré misturadas de etanol e ar para diferentes razões de equivalência e concentrações de água no combustível foram medidas na pressão atmosférica e temperatura de.423K. A concentração de água em volume no combustível variou entre 0,2% (etanol anidro) e 30% em volume, valor máximo no qual uma chama laminar pré-misturada se manteve estável, e a razão de equivalência variou em 0,1 no intervalo entre 0,9 e 1,3. Um queimador tipo Bunsen e a técnica óptica da quimiluminescência do radical OH foram utilizados para a obtenção do perfil da chama. A velocidade de chama laminar, SL, foi calculada utilizando o método da área da chama cônica, e os resultados experimentais foram comparados com os existentes na literatura e com os teóricos obtidos utilizando o modelo de cinética química para oxidação do etanol em altas temperaturas desenvolvido por Marinov utilizando software CHEMKIN-PRO. As simulações numéricas tiveram como dados de entrada as propriedades físicas do combustível no estado de vapor e do ar atmosférico, ambos a 423K e pressão atmosférica, e velocidade teórica da mistura ar e combustível no duto do queimador. Os resultados obtidos mostram se satisfatórios para a faixa de concentração de água de até 20% no combustível. A adição de água ao combustível produziu uma redução da velocidade propagação da chama em cerca de; 5,8% para a concentração de 5% de água no combustível, 9,0 para 10% de água no combustível, 15% para 15 de água no combustível e 21,7% para 20% de água no combustível.

Em linhas gerais, esta dissertação de mestrado envolve o estudo de padrões de escoaomento durante a ebulição convectiva em microcanais. Resultados experimentais foram levantados para um tubo com diâmetro de 2,32 mm durante a evaporação convectiva dos refrigerantes R134a e R245fa. Para a investigação, técnicas experimentais e de análise foram desenvolvidas. A caracterização dos padrões de escoamento envolveu o tratamento simultâneo de sinais provenientes dos seguintes dispositivos: um par diodo/sensor-laser tendo um tubo transparente entre eles no interior do qual ocorre o escoamento bifásico; um transdutor de pressão piezo-elétrico de tamanho reduzido com o objetivo de determinar a variação local da pressão do escoamento; e de um micro-termopar em contato com o fluido refrigerante. A técnica de tratamento de dados utilizada envolve a aglomeração progressiva de dados que apresentem características médias similares através do algoritmo k-means. Os sinais de pressão, intensidade de radiação e temperatura foram adquiridos simultaneamente a uma freqüência de 25 kHz. Imagens simultâneas do escoamento bifásico a uma velocidade de captura em torno de 10.000 imagens/s foram levantadas através de uma câmera de filmagem rápida (até 100.000 imagens/s), e os padrões de escoamentos observados contrastados aos resultados fornecidos pelo método proposto. Baseado nesta análise, mapas de escoamento foram propostos, os quais incorporaram não apenas critérios subjetivos como a visualização, mas também objetivos como as variações transientes da pressão local do escoamento e da morfologia do escoamento através do seu efeito na dispersão da radiação emitida pelo foto-diodo. Os resultados previstos pelo método objetivo apresentam concordância razoável com os dados caracterizados com base em visualizações. Adicionalmente, características de bolhas alongadas foram determinadas. / The present research has been focused on the study of flow patterns inside channels of diameter less than 3 mm during the convective evaporation of refrigerants such as R134a and R245fa. For the investigation of such topics, experimental techniques and methods of analysis of results were developed. A broad database was gathered in an experimental test facility. The characterization of flow patterns involved the simultaneous processing of signals from the following devices: a pair diode / laser-sensor having a transparent tube between them, within which occurs the two-phase flow; a micro piezoelectric pressure transducer to determine the local variation of pressure of the flow and a micro-thermocouple fixed within the fluid. The technique used in data processing involves the gradual agglomeration of data having similar average characteristics; this method was developed based on the k-means clustering algorithm. The signals from the transducers were acquired simultaneously at a frequency of 25 kHz. The program for the acquisition and for processing of the signals was developed using LabView. Simultaneous images of two-phase flow at a speed of capture around 10,000 images / s were obtained through a high speed camera and the observed flow patterns were contrasted to the results provided by the objective method. Based on this analysis, flow maps were proposed, which incorporate not only subjective criteria such as visualization, but also objective criteria like the transient variations of local pressure of the flow, temperature of the fluid and the effect of the flow morphology based on the dispersion of light which effect was captured by the photo-diode. The maps obtained by the objective method were compared against flow pattern segregated based on visualization and a reasonable agreement was obtained. Besides the elongated bubble characteristics were determined.

Orientador: Luiz Roberto Carrocci / Banca: Ivonete Ávila / Banca: Maurício Zanardi Araújo / Banca: Marco Antonio Rosa do Nascimento / Banca: Valesca Alves Correa / Resumo: O leito fluidizado é uma importante operação unitária de grande interesse na indústria química, geração de energia, extração e refino do petróleo, entre outros. As muitas aplicações do leito fluidizado na indústria e a demanda constante de melhorias na eficiência do processo de fluidização, tem produzido um aumento no número de trabalhos envolvendo simulações numéricas do escoamento gás-sólido em leito fluidizado através de CFD. Na simulação numérica, ainda não existe consenso na seleção de um modelo de arrasto e na definição de um valor adequado para o coeficiente de restituição. Nesta tese investiga-se a influência de diferentes modelos de arrasto e do coeficiente de restituição; sobre a fluidodinâmica do escoamento gás-sólido em leito fluidizado, através de simulações numéricas. Para isto é usado um modelo baseado na abordagem Euleriana-Euleriana considerando um leito fluidizado bidimensional, o qual será resolvido usando o codigo livre MFIX. Para avaliar os modelos de arrasto, os resultados numéricos são validados por meio da determinação dos dados experimentais de Taghipour et al. (2005), o qual apresenta parâmetros de projeto de leito fluidizado, tais como a queda de pressão, expansão do leito, velocidade mínima de fluidização e distribuição das fases. Estudou-se também o efeito do coeficiente de restituição sobre a fluidodinâmica do escoamento, através do desenvolvimento de uma correlação em função da fração de vazio. Os resultados demonstram que o modelo de arrasto têm uma influência muito grande na determinação das variáveis de projeto de leito fluidizado. No estudo, um grupo de modelos representam adequadamente a fluidodinâmica, sendo que o modelo de arrasto de Gidaspow (1994) foi o que mais se aproximou aos resultados experimentais. Foram desenvolvidas simulações ... ( Resumo completo, clicar acesso eletrônico abaixo) / Abstract: Fluidized bed is an important unit operation of great interest in petroleum, chemical, metallurgical and energy industries. The large number of applications of fluidized bed in the industry and constant demand for improvements in the efficiency of the fluidization process, has produced an increase in the number of studies involving numerical simulations of gas-solid flow in fluidized bed through CFD. In numerical simulation, there is not yet an agreement on the selection of a drag model and in defining a suitable value for the coefficient of restitution. The purpose of this work was to study the influence of different drag models and restitution coefficient on the hydrodynamics of gas-solid flow in fluidized-bed, through numerical simulations. It is used a model based on the Eulerian-Eulerian approach considering a two-dimensional fluidized bed, which will be resolved using the open source MFIX code. To evaluate drag models, numerical results will be validated through the determination of the experimental data of Taghipour et al. (2005), which presents fluidized bed design parameters, such as fluidization curve, bed expansion, minimum fluidization velocity and void fraction profiles. It was also studied the effect of restitution coefficient on the hydrodynamics, through the development of a correlation in function of void fraction. The results demonstrated that drag model have an great influence in fluidized bed design parameters determination. In this study, a group of models were able to represent adequately the hydrodynamics, however Gidaspow (1994) drag model was the most accurate to represent the experimental data. Simulations were carried out considering restitution coefficient as a function of void fraction and using Gidaspow (1994) drag model. The simulation results obtained was able to improve the prediction of void fraction profile / Doutor

O presente trabalho tem como objetivo fundamental a construção de um sistema de identificação capaz de diagnosticar em tempo real as diferentes configurações de escoamentos bifásicos horizontais. É importante ressaltar que o desenvolvimento deste know-how é capital para a operação eficaz de instalações de manipulação e ou transporte de fluidos multifásicos, e representa, hoje, um dos grandes desafios nas indústrias do petróleo e termonuclear. O princípio de funcionamento do sistema proposto baseia-se nos sinais captados por um sensor de pressão flutuante de resposta rápida, e no seu pósprocessamento com auxílio da transformada de Gabor e de uma rede neural convenientemente treinada. A implementação é tal que a operação de diagnóstico pode ser feita on-line, desde a aquisição dos sinais até o pósprocessamento. Resultados experimentais foram obtidos no circuito experimental do NETeF - Núcleo de Engenharia Térmica e Fluidos da USP - Universidade de São Paulo, para uma secção de testes horizo ntal com 12 m de comprimento e diâmetro interno de 30 mm. Em específico foram ensaiados os seguintes padrões de escoamento ar-água: estratificado liso, ondulado, intermitente, anular e a bolhas. Os resultados mostram que, dependendo dos limites de detecção pré-estabelecidos, todos o principais padrões de escoamento bifásico horizontal são identificados corretamente. / The fundamental objective of this work is the construction of an identification system capable of diagnosing in real time different configurations of horizontal two-phase flow patterns. It is important to emphasize that the development of this know-how is capital to the efficient operation of facilities for manipulation and transportation of multiphase fluids, and represents, today, one of the most important challenges in the oil and thermonuclear industries. The working principle of the proposed system is based on the signals acquired by a rapid response fluctuating pressure sensor, and on its post processing through Gabor Transform and on a conveniently trained artificial neural network. The implementation is accomplished in way that the diagnosis operation is performed on-line, from signal acquisition to post-processing. Experimental results were obtained on the experimental circuit at NETeF - Núcleo de Engenharia Térmica e Fluidos of USP - Universidade de São Paulo, at São Carlos, using a horizontal test section, with 12 m length and 30 mm internal diameter. Experiments were done with the following air-water flow patterns: stratified smooth, wavy, intermittent, annular, and bubbly. Results show that, depending on the preset detection limits, all the main horizontal two phase flow patterns were correctly identified.

Made available in DSpace on 2014-08-13T14:50:37Z (GMT). No. of bitstreams: 0 Previous issue date: 2013-11-01Bitstream added on 2014-08-13T18:01:37Z : No. of bitstreams: 1 000741253_20150101.pdf: 141727 bytes, checksum: bf2e02cb75454b0d8c9536e890144bde (MD5) Bitstreams deleted on 2015-01-05T11:00:41Z: 000741253_20150101.pdf,Bitstream added on 2015-01-05T11:01:44Z : No. of bitstreams: 1 000741253.pdf: 2376805 bytes, checksum: d8aaa4bb45713ee7c22ca08294c74e07 (MD5) / O leito fluidizado é uma importante operação unitária de grande interesse na indústria química, geração de energia, extração e refino do petróleo, entre outros. As muitas aplicações do leito fluidizado na indústria e a demanda constante de melhorias na eficiência do processo de fluidização, tem produzido um aumento no número de trabalhos envolvendo simulações numéricas do escoamento gás-sólido em leito fluidizado através de CFD. Na simulação numérica, ainda não existe consenso na seleção de um modelo de arrasto e na definição de um valor adequado para o coeficiente de restituição. Nesta tese investiga-se a influência de diferentes modelos de arrasto e do coeficiente de restituição; sobre a fluidodinâmica do escoamento gás-sólido em leito fluidizado, através de simulações numéricas. Para isto é usado um modelo baseado na abordagem Euleriana-Euleriana considerando um leito fluidizado bidimensional, o qual será resolvido usando o codigo livre MFIX. Para avaliar os modelos de arrasto, os resultados numéricos são validados por meio da determinação dos dados experimentais de Taghipour et al. (2005), o qual apresenta parâmetros de projeto de leito fluidizado, tais como a queda de pressão, expansão do leito, velocidade mínima de fluidização e distribuição das fases. Estudou-se também o efeito do coeficiente de restituição sobre a fluidodinâmica do escoamento, através do desenvolvimento de uma correlação em função da fração de vazio. Os resultados demonstram que o modelo de arrasto têm uma influência muito grande na determinação das variáveis de projeto de leito fluidizado. No estudo, um grupo de modelos representam adequadamente a fluidodinâmica, sendo que o modelo de arrasto de Gidaspow (1994) foi o que mais se aproximou aos resultados experimentais. Foram desenvolvidas simulações ... ( Resumo completo, clicar acesso eletrônico abaixo) / Fluidized bed is an important unit operation of great interest in petroleum, chemical, metallurgical and energy industries. The large number of applications of fluidized bed in the industry and constant demand for improvements in the efficiency of the fluidization process, has produced an increase in the number of studies involving numerical simulations of gas-solid flow in fluidized bed through CFD. In numerical simulation, there is not yet an agreement on the selection of a drag model and in defining a suitable value for the coefficient of restitution. The purpose of this work was to study the influence of different drag models and restitution coefficient on the hydrodynamics of gas-solid flow in fluidized-bed, through numerical simulations. It is used a model based on the Eulerian-Eulerian approach considering a two-dimensional fluidized bed, which will be resolved using the open source MFIX code. To evaluate drag models, numerical results will be validated through the determination of the experimental data of Taghipour et al. (2005), which presents fluidized bed design parameters, such as fluidization curve, bed expansion, minimum fluidization velocity and void fraction profiles. It was also studied the effect of restitution coefficient on the hydrodynamics, through the development of a correlation in function of void fraction. The results demonstrated that drag model have an great influence in fluidized bed design parameters determination. In this study, a group of models were able to represent adequately the hydrodynamics, however Gidaspow (1994) drag model was the most accurate to represent the experimental data. Simulations were carried out considering restitution coefficient as a function of void fraction and using Gidaspow (1994) drag model. The simulation results obtained was able to improve the prediction of void fraction profile

[pt] A produção de óleo diminui e a produção de água aumenta com o passar do tempo na vida de um reservatório. A mistura de óleo e água é geralmente produzida na forma de uma emulsão. A formação da emulsão começa no escoamento bifásico no interior do reservatório e sua estrutura muda na medida em que os líquidos escoam através de tubulações, bombas e válvulas até as instalações de superfície. Durante a produção, as gotas maiores da fase dispersa quebram-se mudando a distribuição do tamanho das gotas. É importante conhecer a distribuição de tamanho de gota da fase dispersa da emulsão a fim de projetar as unidades de separação e prever as quedas de pressão ao longo do escoamento. O objetivo do presente trabalho é o estudo do processo de quebra de gotas em capilares retos e válvulas agulha a fim de prever o tamanho das gotas resultantes em função das condições de escoamento. O principal desafio é entender como os diferentes parâmetros operacionais de escoamento afetam o processo de quebra. Duas bancadas experimentais foram utilizadas para o estudo do processo. Na primeira, foi realizada uma análise paramétrica de formação de emulsões em um escoamento laminar através de capilares retos. Os experimentos foram realizados utilizando duas seringas conectadas através de um capilar. A emulsão óleo-emágua foi forçada a uma passagem de ida e vinda através do capilar. O diâmetro médio de gota e a superfície específica da fase dispersa foram obtidos em função da vazão de injeção, taxa de cisalhamento, tempo de residência e trabalho dissipado na parede do capilar. Como esperado, o diâmetro médio da fase dispersa diminui com o aumento do tempo de cisalhamento atingindo um valor assintótico e possui uma grande dependência com a taxa de cisalhamento na parede do capilar. Na segunda bancada, foi realizada uma análise paramétrica de um escoamento turbulento de emulsão óleo-em-água através de uma válvula agulha. Os experimentos foram realizados utilizando uma bomba helicoidal para controlar a vazão através da válvula. O diâmetro médio de gota e a superfície específica da fase dispersa a montante e a jusante da válvula foram obtidos em função da queda de pressão, vazão e taxa de dissipação de energia na válvula. O diâmetro médio da fase dispersa diminui e a superfície específica aumenta com o aumento da queda de pressão na válvula agulha até atingir um valor assintótico. / [en] Oil production decreases and water production increases as time goes by in the life of a hydrocarbon reservoir. The mixture of oil and water is usually produced as an emulsion. Emulsion formation starts in the two-phase flow inside the reservoir. The emulsion structure changes as it flows through pipes, pumps and valves up to the surface facilities. During all stages, large drops of the dispersed phase break up leading to smaller drops. It is important to know the droplet size distribution of the dispersed phase in order to design separation units and predict the pressure drop along the flow. The aim of the this work is to study the droplets break-up process that takes place in capillaries and in a needle valve in order to make predictions of the size of the resulting droplets that emerge from this process. The main challenge is to understand how the different operating flow parameters affect the break up process. In order to achieve this goal, two laboratory scale experimental set-ups have been used. In the first experiment, we conducted a parametric analysis of oil-water emulsion formation in laminar flow through straight capillaries. The experiments were carried out using two syringe pumps connected by a double-hubbed capillary pipe. The oil-water emulsion is forced back and forth through the pipe. The mean diameter and the specific surface area of the dispersed phase were obtained as a function of flow rate, shear rate, residence time and rate of energy dissipation at the capillary wall. As expected, keeping all other variables fixed, the dispersed phase mean diameter decreases with the shearing time, reaching an asymptotic value, which was a strong function of the shear rate at the capillary wall. Secondly, we conducted a parametric analysis of turbulent oil-in-water emulsion flow through a needle valve. The experiments were carried out using a helicoidal pump to control the flow rate through the needle valve. The mean diameter and the specific surface area of the dispersed phase upstream and downstream of the valve were obtained as a function of the pressure drop in the valve, flow rate, and rate energy dissipation of the flow. The dispersed phase mean diameter falls and the specific surface area rises with the pressure drop in the valve until reaching an asymptotic value.

Submitted by Marco Antônio de Ramos Chagas (mchagas@ufv.br) on 2017-02-08T16:17:33Z No. of bitstreams: 1 texto completo.pdf: 6182702 bytes, checksum: b0debc9cf8e3c2e2e9811d5905d8e473 (MD5) / Made available in DSpace on 2017-02-08T16:17:33Z (GMT). No. of bitstreams: 1 texto completo.pdf: 6182702 bytes, checksum: b0debc9cf8e3c2e2e9811d5905d8e473 (MD5) Previous issue date: 2002-12-20 / Conselho Nacional de Desenvolvimento Científico e Tecnológico / Desenvolveu-se uma metodologia que permite realizar o balanço hídrico de bacias hidrográficas e obter, a partir do modelo digital de elevação do terreno, o hidrograma de escoamento superficial, a vazão máxima e seu tempo de ocorrência, e o volume escoado superficialmente, considerando escoamento em encostas ou em canais. Para a utilização do modelo, o usuário deve fornecer imagens relativas ao modelo digital de elevação do terreno, à rede de drenagem e às direções de escoamento na bacia hidrográfica e de informações relativas à precipitação, à interceptação pela cobertura vegetal, ao armazenamento superficial, à rugosidade do terreno e do solo. A taxa de produção de vazão de cada célula da bacia é obtida pela diferença entre a precipitação e a infiltração, obtida por meio da equação de Green-Ampt modificada por Mein-Larson. A modelagem do escoamento superficial é realizada utilizando o modelo de ondas cinemáticas. Uma vez informado o local para o qual se deseja obter o hidrograma, o modelo identifica se o escoamento é em canais ou na encosta. O hidrograma resultante de cada célula é obtido pela soma dos hidrogramas advindos da área de contribuição com o hidrograma gerado na mesma. Para a avaliação do modelo, instrumentou-se uma bacia hidrográfica na localidade de Palmital, Viçosa-MG, monitorando-se a precipitação e o escoamento superficial. Visando facilitar a utilização da metodologia elaborou-se um software no ambiente de programação Delphi. A fim de se realizar a análise de sensibilidade do modelo, procedeu-se a simulações considerando variações nos valores de interceptação pela cobertura vegetal, armazenamento superficial, rugosidade do terreno, condutividade hidráulica do solo saturado, umidade de saturação e potencial matricial na frente de umedecimento. Os valores obtidos nas simulações realizadas com o modelo foram comparados com os valores obtidos na bacia, com o Método do Número da Curva, com o Método Racional e com o WEPP. Os resultados das simulações apresentaram diferenças em relação aos valores obtidos na bacia. Os resultados obtidos mostraram a sensibilidade do modelo em relação aos fatores citados. / The developed methodology allows the obtaining of the water balance in watersheds, and, based on a digital model of ground elevation, of the runoff hydrogram, the maximum outflow, its duration, and of the surface runoff volume, considering runoff down slopes or concentrated in channels. To use the model, images relating to the digital model of ground elevation, the drainage system, and flow direction within the watershed are required as well as information on rainfall, interception of the canopy, surface storage, and the soil and ground roughness. The overland flow production rate of each watershed cell is obtained by the difference between rainfall and infiltration, calculated by the Green-Ampt equation, modified by Mein-Larson. The runoff is modeled with the kinematic wave model. Once information on the location of the required hydrogram is provided, the model identifies whether the runoff occurs in channels or down the slopes. The hydrogram of each cell is obtained by the sum of hydrograms derived from the catchment area by the hydrogram created by the same. For an evaluation run of the model, rainfall and runoff were measured in a watershed in Palmital, Viçosa, State of Minas Gerais. Software of the Delphi programming line was worked out to facilitate the use of the methodology. For the sensitivity analysis of the model, simulations were set up that considered variations in the values of interception by the canopy, surface storage, ground roughness, saturated hydraulic conductivity, saturation moisture and matric potential in the wetting front. The simulation values obtained by the model were compared to the values observed in the watershed, by the Curve Number Method, the Rational Method and the WEPP. Simulation results presented differences in relation to the values obtained in the watershed. Results confirmed the sensitivity of the model in relation to the mentioned factors. / Tese importada do Alexandria

Dissertação (mestrado) - Universidade Federal de Santa Catarina, Centro Tecnológico, Programa de Pós-Graduação em Engenharia Química, Florianópolis, 2013. / Made available in DSpace on 2014-08-06T17:09:28Z (GMT). No. of bitstreams: 1 323813.pdf: 923726 bytes, checksum: a8015d8b04c2edbedf60057ec1e4ebba (MD5) Previous issue date: 2013 / Esta dissertação de mestrado originou-se da necessidade de avaliação de estratégias numéricas que permitam uma adequada predição in silico dos fenômenos de transferência de calor e massa em sistemas multicomponentes e multifásicos. Em particular, o foco do trabalho consistiu na implementação de modelos robustos em ferramenta de CFD. Para tanto, as equações de Maxwell-Stefan em conjunto com a lei de Fick generalizada foram codificadas na linguagem de programação C e o código gerado foi dinamicamente acoplado ao código comercial ANSYS® CFD (FLUENT®), versão 14.0. Determinou-se, então, a composição e a temperatura de equilíbrio em um sistema vapor-líquido, sendo cada fase composta por uma mistura de quatro hidrocarbonetos puros (metano, n-pentano, n-hexano e n-octano). Levou-se em consideração todas as características inerentes a misturas multicomponentes (a exemplo da correção para altas taxas e da contribuição de todos os gradientes de composição para a taxa de transferência da espécie i), o que introduz uma complexidade considerável ao sistema quando comparada à solução de escoamentos envolvendo misturas binárias. Os resultados obtidos nos estudos com a ferramenta de CFD foram verificados com predições obtidas do simulador de processos comercial em estado estacionário denominado PRO/II®, versão 8.2. Um estudo preliminar com a mistura binária água/ar foi conduzido com o objetivo de validar os resultados obtidos com o código CFD por meio de comparação direta com dados fornecidos pela carta psicrométrica padrão. <br> / Abstract : This dissertation was originated from the need of the evaluation of numerical strategies that may allow an adequate prediction in silico of heat and mass transfer phenomena in multicomponent and multiphase systems. In particular, the aim of the work has been to implement robust models in a CFD tool. Therefore, the Maxwell-Stefan's equations in conjunction with the generalized Fick's law have been codified in the C programming language and the code generated has been dynamically coupled to the commercial code ANSYS® CFD (FLUENT®), version 14.0. It has been determined, then, the equilibrium composition and temperature in a vapor-liquid system, in which each phase has been composed by a mixture of four pure hydrocarbons (methane, n-pentane, n-hexane and octane). All the inherent characteristics of multicomponent mixtures (e.g., the high flux correction and the contribution of all composition gradients to the transport rate of species i) have been taken into account, which introduces considerable complexity to the system when compared to the solution of flow involving binary mixtures. The results obtained in the studies conducted with the CFD code have been verified with predictions obtained from the steady state process simulator named PRO/II®, version 8.2. A preliminary study with the binary mixture water/air has been carried out with the aim to validate the results obtained with the CFD code by means of direct comparison with data provided by standard psychrometric chart.