Fenômenos aeroelásticos podem levar à drástica redução na vida útil de uma aeronave ou ainda resultam em danos severos à estrutura. Para manter as respostas dinâmicas em níveis aceitáveis técnicas como as estruturas adaptativas têm sido aplicadas. Este conceito explora a integração entre os elementos ativos (atuadores e sensores) e o controlador à estrutura. Dentre os materiais próprios para uso em estruturas adaptativas estão os fluidos eletro-reológicos e magneto-reológicos que tem se mostrado como um dos mais promissores materiais ativos. Estes materiais apresentam rápidas mudanças nas suas propriedades reológicas devido à ação de um campo elétrico ou magnético. Para sua incorporação em uma estrutura é utilizada uma viga sanduíche que tem seu comportamento dinâmico modelado através do método GHM para incorporar a dependência da freqüência dos fluidos ER/MR em um modelo estrutural no domínio do tempo. Através do acoplamento deste modelo com o método da malha de vórtices, é possível estudar a resposta aeroelástica temporal. Também é analisada a eficiência dos fluidos ER/MR no atraso da ocorrência de flutter. Isto é feito com o auxílio do método PK que determina a velocidade crítica de flutter. / Aeroelastic phenomena can lead to a drastic reduction in the fatigue life of aircraft or result in severe structural damage. To keep the dynamical responses at acceptable levels techniques such as the so-called adaptive structures have been adopted. This approach integrates active elements and controllers (actuators and sensors) to the structure. Among the materials suitable for adaptive structures are the electro-rheological (ER) and magneto-rheological fluids which are some of the most promising active materials. This kind of materials presents change in their rheological properties due to action of an external field, such as electrical or magnetic. In order to integrate these kind of fluids in the structure a sandwich beam with ER/MR fluids core is studied. The dynamical behaviour is modelled through a GHM method to incorporate the frequency dependence of the ER/MR fluids in a structural time domain model. By coupling this model to a vortex lattice model, it is possible to study the aeroelastic response in time domain. The ER/MR fluids efficiency to delay the flutter occurrence is also studied by using a PK-method that determines a critical velocity of flutter.

As complexidades das configurações, da geometria e da topologia dos escoamentos encontrados em sistemas tais como câmaras de combustão de turbo-máquinas requerem métodos numéricos eficientes para modelar os escoamentos turbulentos reativos existentes. Este trabalho de modelagem numérica está relacionado com o desenvolvimento de ferramentas numéricas eficientes necessárias ao estudo dos escoamentos turbulentos reativos nessas aplicações. Na primeira parte desta tese, após relembrar os elementos de base que permitem a descrição de um escoamento turbulento reativo prée-misturado, são apresentadas as diversas etapas que presidem a elaboração do programa computacional. Inicialmente, introduz-se a técnica de compressibilidade artificial, que consiste em modificar as equações governantes acrescentando termos não-estacionários no pseudo-tempo, destinados a introduzir uma velocidade do som finita no domínio de cálculo. Em seguida é exposta a discretização por volumes finitos das equações governantes em malhas não-estruturadas bidimensionais híbridas, que contêm elementos triangulares e/ou quadriláteros. A adição de termos de dissipação artificial bem com o uso de um método do tipo Runge-Kutta para integrar os termos temporais no pseudo-tempo das equações governantes, que completam o método, são comentados em seguir. O programa de cálculo, desenvolvido na linguagem C++, foi validado considerando-se o problema de uma chama plana turbulenta, propagando-se livremente em uma turbulência "congelada". Este fenômeno foi modelado com um termo de produção química médio do tipo flamelet, no qual um corte representando extinção foi introduzido com a finalidade de selecionar uma só velocidade possível de propagação. Os resultados obtidos para uma chama lenta (0,5 m/s) ou rápida (10 m/s) foram comparados com aqueles obtidos na literatura. A precisão dos resultados é bastante satisfatória e os gradientes elevados de pressão estática, característicos do termo fonte empregado, são capturados corretamente, tanto no que se refere a sua intensidade quanto à sua localização, dentro da frente de chama média. Uma análise de sensibilidade e precisão da malha foi igualmente conduzida, assim como um estudo das características de convergência em função do valor do coeficiente de compressibilidade artificial. Tendo sido validado, o programa de cálculo foi, em seguida, utilizado com a finalidade de estudar esteiras inertes ou reativas encontradas à jusante de um obstáculo de seção transversal triangular, colocado no meio de um escoamento turbulento em um canal. Um modelo de turbulência do tipo k-e, acoplado a dois tipos de leis de parede foi utilizado. O desprendimento turbilhonar, típico dos escoamentos de esteiras inertes considerados, foi previsto corretamente com, em particular, uma freqüência de Strouhal de desprendimento com boa concordância com os resultados experimentais disponíveis. As esteiras reativas foram calculadas com um modelo de flamelet, onde os transportes turbulentos da variável de avanço da combustão foram calculados por aproximação do tipo gradiente. Três formas do termo de produção química médio, existentes na literatura, foram utilizados. Aparentemente, a estrutura da esteira reativa prevista é fortemente dependente da formulação adotada para o termo fonte médio. Embora os resultados obtidos sejam comparados àqueles apresentados por outro autores, em uma configuração de escoamento similar, a comparação com resultados experimentais mostra que é necessário uma melhoria na modelagem física de forma a prever este tipo de escoamento reativo de maneira satisfatória.

O presente trabalho consiste no estudo, implementação e análise de métodos não-oscilatórios de alta ordem de precisão em malhas não-estruturadas através da técnica de volumes finitos. A motivação deste trabalho surge da necessidade de se ter simulações acuradas para os escoamentos com altos números de Mach e, portanto, com fortes descontinuidades, que se pretende analisar. Nos últimos anos, o grupo de dinâmica dos fluidos computacional do Instituto de Aeronáutica e Espaço vem desenvolvendo ferramentas computacionais capazes de capturar com precisão as descontinuidades, ondas de choque e superfícies de contato, presentes nos escoamentos aerodinâmicos de interesse. O objetivo principal do presente trabalho de pesquisa consiste em estudar os esquemas essencialmente não-oscilatórios, do inglês essentially non-oscillatory (ENO), e os esquemas essencialmente não-oscilatórios ponderados, do inglês weighted essentially non-oscillatory (WENO). Os esquemas ENO e WENO foram desenvolvidos com o propósito de capturar com acurácia descontinuidades presentes em problemas governados por equações diferenciais parciais hiperbólicas. No desenvolvimento deste trabalho, vai-se utilizar uma ferramenta numérica de simulação de escoamentos compressíveis em duas dimensões já desenvolvida pelo atual grupo de pesquisa. Esta ferramenta resolve as equações governantes da mecânica dos fluidos através da técnica de volumes finitos com uma formulação que atribui os valores das propriedades conservadas aos centróides dos volumes de controle. O código numérico trata de malhas com qualquer geometria poligonal, como triângulos e quadriláteros, além de malhas híbridas compostas por ambos os elementos. As equações de Euler em duas dimensões representam os escoamentos de interesse no presente trabalho. O processo de reconstrução dos esquemas ENO e WENO é descrito em detalhes para polinômios lineares, quadráticos e cúbicos e, portanto, para esquemas numéricos de segunda, terceira e quarta ordem de precisão. Estes polinômios interpolam variáveis primitivas em pontos de quadraturas de Gauss utilizando moléculas de cálculo determinadas por vizinhanças compostas por volumes de controle adjacentes a um volume principal. Métodos upwind são utilizados para se aproximar os valores dos fluxos nas faces dos volumes de controle. As discretizações temporais são efetuadas por esquemas do tipo Runge-Kutta de alta ordem. Uma ferramenta multigrid é utilizada para acelerar a convergência dos esquemas para soluções de estado estacionário e uma técnica de refinamento adaptativo de malha é utilizada para melhorar a qualidade da malha em regiões com altos gradientes das propriedades do escoamento. Aplicações em escoamentos aerodinâmicos de interesse são realizadas e os resultados são comparados com a literatura e com os métodos previamente implementados pelo grupo de pesquisa.

The objective of this project is to study three dimensional aerodynamic flows over realistic aerospace configurations, such as the first Brazilian Satellite Launcher (VLS) in its first stage flight form. In order to achieve this objective, computational fluid dynamics and three dimensional adaptive mesh refinement techniques are used. The project focus in finite element and finite volume techniques to simulate the flows of interest. Both the Taylor-Galerkin finite element method and the Jameson finite volume method are used in the simulations. A new interpretation of MacCormack's method in an unstructured grid context is presented. The analysis, implementation and validation of the adaptive mesh refinement tools in the present context is part of the work developed. The use of meshes with hanging-nodes in a finite volume code is also described. The Spalart and Allmaras turbulence model is implemented in the finite volume code in order to account for turbulent effects. In the first part of this project, inviscid flows over the VLS ranging from subsonic to supersonic regime are studied using the finite element and finite volume techniques and, in the following part, turbulent viscous flow cases are analyzed using the finite volume code.

Neste trabalho foi realizado um estudo utilizando os modelos de turbulência Simulação das Grandes Escalas, Large Eddy Simulation (LES), e Simulação dos Vórtices Desprendidos, Detached Eddy Simulation (DES), para simular o escoamento compressível interno em um duto contendo válvulas controladoras na saída dos gases de combustão de um reator de craqueamento catalítico fluido, com o objetivo de comparar o desempenho numérico e computacional de ambas as técnicas. Para isso foi utilizado um programa comercial de dinâmica dos fluidos computacional, Computational Fluid Dynamics (CFD), que possui em seu código os dois modelos de turbulência. / In the present work a study was conducted using Large Eddy Simulation (LES) and Detached Eddy Simulation (DES) turbulence models in order to simulate the internal compressible flow in a duct containing the flue gas discharge control valves of a fluid catalytic cracking reactor so as to compare the numerical and computational behavior of both techniques. A commercial Computational Fluid Dynamics (CFD) software, which includes these turbulence models in its code, was used.

O concreto é um dos materiais de construção civil mais versáteis, podendo se adaptar a formas diversas quando em seu estado fresco e resistindo a grandes cargas de compressão em seu estado rígido. No entanto, em estruturas mais densamente armadas ou com geometrias mais complexas, pode-se apresentar dificuldades para a moldagem, causando falhas no preenchimento da forma, o que reduz a capacidade resistente da peça e sua vida útil. Neste trabalho foi utilizado o método de partículas lagrangeanas Moving Particle Semi-Implicit (MPS) como base para um simulador para o estudo do comportamento do escoamento de pastas e argamassas cujo comportamento pode se aproximado por modelos reológicos como Bingham e Herschel-Bulkley. Foram propostos módulos para a simulação da viscosidade não-newtoniana, variação térmica no processo de cura e modelagem de turbulência. Para modelar a variação de viscosidade de um fluido não-newtoniano foi utilizado o modelo de Herschel-Bulkley. O modelo de Herschel-Bulkley possui uma singularidade para taxas de deformação muito pequenas, que resulta em valores de viscosidade infinitas, dificuldade contornada pela solução proposta por Papanastasiou (1987). Na modelagem térmica foram analisados dois modelos de dissipação, sendo um original do método e outro calculado por meio do divergente do gradiente utilizando os modelos de partículas e que teria melhores resultados para o cálculo da dissipação térmica. Também foi modelada a convecção térmica, utilizando para isso a hipótese de Boussinesq. A reação de hidratação do concreto foi modelada utilizando uma equação do tipo Hill para representar a elevação de temperatura obtida por meio um ensaio adiabático. Para complementar as simulações, foi utilizado o modelo de turbulência Detached Eddy Simulation (DES), baseado no método Large Eddy Simulation (LES), que utiliza um modelo de parede para simular a interação do fluido. Para a implementação deste modelo de turbulência foi proposto um algoritmo para o cálculo da distância da partícula de fluido à parede. Este algoritmo utiliza estruturas de dados já existentes no método de partículas de modo que sua execução requer muito menos recursos que a busca binária. Apesar do trabalho ter se guiado pela simulação do concreto fresco, que é um material particularmente complexo, outros escoamentos encontrados na engenharia civil também podem ser beneficiados pelo método, como os estudos do escoamento em sistemas prediais de água e esgoto, do escoamento e prevenção de erosão em rios e córregos, do escorregamento de encostas, dos reatores para depuração de esgotos, entre outros. / The concrete is one of the more versatile civil construction materials, it can adapt to various forms when in its fresh state while resisting to high compression loads in its rigid state. However, in some cases like densely reinforced concrete structures and complex geometry concrete structures can present issues to the casting, and failure in proper form filling can occur. These failures can reduce the resistance and the lifetime of the structure. This work used a simulator based on the lagrangean particle method called Moving Particle Semi-Implicit (MPS) to study the concrete behavior in its distinctive characteristics. Also, this work proposed modules to simulate the non-Newtonian viscosity, the thermal process of concrete cure and to model the turbulent flow. To model the variation of viscosity of a non-Newtonian fluid, the Herschel-Bulkley model, which relates the shear stress with the strain rate, was applied. The Herschel-Bulkley model has a singularity at low strain rates, which results in infinite viscosities. To overcome this issue, Papanastasiou (1987) proposed a modification in the model in order to eliminate the singularity. For the thermal modeling, two models for thermal dissipation were compared, the original method and other calculated from the divergence of gradient using the differential operators for the particle model and that could present improved results for the thermal dissipation calculation. Also, the thermal convection was modeled using the Boussinesq hypothesis. The hydration reaction of the concrete was modeled using a Hill type equation to reproduce the temperature elevation. In addition, a Detached Eddy Simulation (DES) based turbulence model with a simple wall model in the interaction of wall and fluid was applied in the simulations. To improve the turbulence model, an algorithm to calculate the distance between fluid and the nearest wall particle was proposed. The algorithm uses already calculated information from particles so that the execution requires much less resources than a binary search. Although the work has been focused on to the modeling of fresh concrete simulation, which is a particularly complex material, other flows found in civil engineering can also be benefited by the method, such as studies of drainage in water and sewage systems, drainage and prevention of erosion into rivers and streams, prevention and damage mitigation of landslides, reactors for sewage treatment among many others.

O concreto é um dos materiais de construção civil mais versáteis, podendo se adaptar a formas diversas quando em seu estado fresco e resistindo a grandes cargas de compressão em seu estado rígido. No entanto, em estruturas mais densamente armadas ou com geometrias mais complexas, pode-se apresentar dificuldades para a moldagem, causando falhas no preenchimento da forma, o que reduz a capacidade resistente da peça e sua vida útil. Neste trabalho foi utilizado o método de partículas lagrangeanas Moving Particle Semi-Implicit (MPS) como base para um simulador para o estudo do comportamento do escoamento de pastas e argamassas cujo comportamento pode se aproximado por modelos reológicos como Bingham e Herschel-Bulkley. Foram propostos módulos para a simulação da viscosidade não-newtoniana, variação térmica no processo de cura e modelagem de turbulência. Para modelar a variação de viscosidade de um fluido não-newtoniano foi utilizado o modelo de Herschel-Bulkley. O modelo de Herschel-Bulkley possui uma singularidade para taxas de deformação muito pequenas, que resulta em valores de viscosidade infinitas, dificuldade contornada pela solução proposta por Papanastasiou (1987). Na modelagem térmica foram analisados dois modelos de dissipação, sendo um original do método e outro calculado por meio do divergente do gradiente utilizando os modelos de partículas e que teria melhores resultados para o cálculo da dissipação térmica. Também foi modelada a convecção térmica, utilizando para isso a hipótese de Boussinesq. A reação de hidratação do concreto foi modelada utilizando uma equação do tipo Hill para representar a elevação de temperatura obtida por meio um ensaio adiabático. Para complementar as simulações, foi utilizado o modelo de turbulência Detached Eddy Simulation (DES), baseado no método Large Eddy Simulation (LES), que utiliza um modelo de parede para simular a interação do fluido. Para a implementação deste modelo de turbulência foi proposto um algoritmo para o cálculo da distância da partícula de fluido à parede. Este algoritmo utiliza estruturas de dados já existentes no método de partículas de modo que sua execução requer muito menos recursos que a busca binária. Apesar do trabalho ter se guiado pela simulação do concreto fresco, que é um material particularmente complexo, outros escoamentos encontrados na engenharia civil também podem ser beneficiados pelo método, como os estudos do escoamento em sistemas prediais de água e esgoto, do escoamento e prevenção de erosão em rios e córregos, do escorregamento de encostas, dos reatores para depuração de esgotos, entre outros. / The concrete is one of the more versatile civil construction materials, it can adapt to various forms when in its fresh state while resisting to high compression loads in its rigid state. However, in some cases like densely reinforced concrete structures and complex geometry concrete structures can present issues to the casting, and failure in proper form filling can occur. These failures can reduce the resistance and the lifetime of the structure. This work used a simulator based on the lagrangean particle method called Moving Particle Semi-Implicit (MPS) to study the concrete behavior in its distinctive characteristics. Also, this work proposed modules to simulate the non-Newtonian viscosity, the thermal process of concrete cure and to model the turbulent flow. To model the variation of viscosity of a non-Newtonian fluid, the Herschel-Bulkley model, which relates the shear stress with the strain rate, was applied. The Herschel-Bulkley model has a singularity at low strain rates, which results in infinite viscosities. To overcome this issue, Papanastasiou (1987) proposed a modification in the model in order to eliminate the singularity. For the thermal modeling, two models for thermal dissipation were compared, the original method and other calculated from the divergence of gradient using the differential operators for the particle model and that could present improved results for the thermal dissipation calculation. Also, the thermal convection was modeled using the Boussinesq hypothesis. The hydration reaction of the concrete was modeled using a Hill type equation to reproduce the temperature elevation. In addition, a Detached Eddy Simulation (DES) based turbulence model with a simple wall model in the interaction of wall and fluid was applied in the simulations. To improve the turbulence model, an algorithm to calculate the distance between fluid and the nearest wall particle was proposed. The algorithm uses already calculated information from particles so that the execution requires much less resources than a binary search. Although the work has been focused on to the modeling of fresh concrete simulation, which is a particularly complex material, other flows found in civil engineering can also be benefited by the method, such as studies of drainage in water and sewage systems, drainage and prevention of erosion into rivers and streams, prevention and damage mitigation of landslides, reactors for sewage treatment among many others.

A presente tese trata da análise experimental do efeito da geometria da seção transversal do canal e do desempenho de refrigerantes de reduzido GWP (Global Warming Potential) durante a ebulição convectiva em canais de reduzidas dimensões. A tese inclui ainda um estudo extenso e crítico da literatura sobre métodos de previsão da perda de pressão e do coeficiente de transferência de calor, e sobre estudos experimentais em canais não-circulares e de refrigerantes com reduzido GWP na ebulição convectiva em canais de dimensões reduzidas. Resultados para o coeficiente de transferência de calor e perda de pressão durante a ebulição convectiva foram obtidos para canais com geometrias de seção circular, quadrada e triangular para o refrigerante R134a. Nos testes utilizou-se canais com perímetros internos similares obtidos a partir da conformação de um tubo com diâmetro interno igual a 1,1 mm. No caso do canal circular, dados foram também levantados para os HFOs R1234ze(E) e R1234yf e o hidrocarboneto R600a, fluidos com reduzido GWP. Ensaios foram executados para amplas faixas de fluxos de calor e velocidades mássicas, temperaturas de saturação de 31 e 41°C e títulos de vapor entre 0 e 0,95. Aspectos relacionados aos efeitos da geometria e do fluido refrigerante foram minuciosamente investigados através da análise paramétrica dos resultados. Com base na comparação do banco de dados coletado com os métodos de previsão disponíveis na literatura, constatou-se que estes proporcionam previsões satisfatórias apenas para condições experimentais especificas. Portanto, novos métodos de previsão da perda de pressão e do coeficiente de transferência de calor foram desenvolvidos com base nos dados levantados no presente estudo. Os métodos propostos preveem satisfatoriamente o banco de dados do presente estudo e resultados independentes disponíveis na literatura. Adicionalmente, com base nos resultados levantados, verificou-se que dissipadores de calor baseados em multi-microcanais com canais de seção triangular apresentam desempenho superior comparados a dissipadores com canais quadrados e circulares. / The present thesis concerns an experimental study on the effects of cross-sectional geometry and low GWP refrigerants on the thermal-hydraulic performance for convective boiling inside micro-scale channels. Experimental results for heat transfer coefficient and pressure drop gradient during convective boiling were obtained for circular, square and triangular channels for the fluid R134a. The evaluated channels present the same external perimeter and equivalent diameters of 1.1, 0.977 and 0.835 mm, respectively. In the case of the circular geometry, experimental results were also acquired for the HFOs R1234ze(E) and R1234yf and the hydrocarbon R600a (isobutane), which are fluids with low GWP and null ODP. Experiments were performed for a wide range of heat fluxes and mass velocities, saturation temperatures of 31 and 41°C and vapor qualities up to 0.95. The experimental data were carefully analyzed and discussed based on a parametrical analysis focusing on the effect of the cross-sectional geometry and the working fluid. Subsequently, the experimental data were compared against the most quoted predictive methods from literature. In general, it was verified that none of the predictive methods were able to accurately capture the experimental trends of the overall database. So, new predictive methods for the pressure drop and heat transfer coefficient were developed based on the broad database obtained in the present study. The proposed methods provided satisfactory results not only for the experimental database used for its development, but also for independent databases collected in the literature. Additionally, based on the data obtained in the present study and a performance analysis taken into account pressure drop, heat transfer coefficient and the channel packing factor, triangular cross sectional geometry is recommended for heat sinks.

Submitted by Johnny Rodrigues (johnnyrodrigues@ufcg.edu.br) on 2018-07-04T14:51:20Z No. of bitstreams: 1 NÉLIO JOSÉ LIRA PEREIRA - DISSERTAÇÃO PPGEA 2011..pdf: 15133495 bytes, checksum: d48bf3754fca4f45eda0a85dbd10f392 (MD5) / Made available in DSpace on 2018-07-04T14:51:20Z (GMT). No. of bitstreams: 1 NÉLIO JOSÉ LIRA PEREIRA - DISSERTAÇÃO PPGEA 2011..pdf: 15133495 bytes, checksum: d48bf3754fca4f45eda0a85dbd10f392 (MD5) Previous issue date: 2011-09 / Foram determinadas as propriedades reológicas e termofísicas de polpas de goiaba cv. Paluma, nas concentrações de 8, 10 e 12 °Brix. As amostras foram previamente caracterizadas quanto às características químicas, físicas e físico-químicas. As medidas viscométricas foram realizadas fazendo-se uso de um viscosímetro Brookfíeld, obtendo-se, desta forma dados de viscosidade aparente, torque, tensão de cisalhamento e taxa de deformação, nas temperaturas de 10, 20, 30, 40 e 50 °C. As propriedades termofísicas massa específica, foi determinada pelo método picnométrico; o calor específico pelo método das misturas e a difusividade térmica pelo método do cilindro infinito. A condutividade térmica foi obtida pela relação com a massa específica, calor específico e difusividade térmica, obtidos experimentalmente. Os modelos da Lei da Potência, Casson, Herschel-Bulkley e Mizrahi & Berk foram ajustados aos dados de tensão de cisalhamento em função da taxa de deformação. As propriedades termofísicas também foram estimadas por meio de equações propostas em literatura cujos resultados foram comparados com os resultados experimentais. As polpas de goiaba se caracterizaram como ácidas, apresentaram comportamento pseudoplástico e as viscosidades aparentes foram influenciadas pela temperatura e seguiram uma relação do tipo Arrhenius. O calor específico foi reduzido com o aumento da concentração; a massa específica aumentou, seguindo relação direta com a concentração e inversa, com a temperatura, e foi bem estimada por meio de equações. A condutividade térmica calculada por equações propostas em literatura não resultou em estimativas aceitáveis. / The rheological and thermophysical properties of guava cv. Paluna pulp were determined at concentrations of 8, 10 and 12 °Brix, and temperatures between 10 and 50 °C. The samples, previously characterized with respect to chemical and physical properties, were submitted to viscometric measures, using a Brookfield RVT viscometer to obtain data on viscosity, torque, shear stress and shear rate. Shear stress and shear rate data were fit to the Power Law, Casson, Herschel-Bulkley and Mizrahi & Berk models. Thermophysical properties and the respective determination methodologies were density using the pycnometer method, specifíc heat by the method of mixtures, and thermal diffusiveness using the infinite cylinder method. Thermal conductivity was obtained from the relation with density, specifíc heat and thermal diffusiveness, obtained experimentally. Thermophysical properties were also estimated by equations proposed in the literature and results were compared with experimental results. Guava pulps are characterized as acids and concentration did not affect color, luminosity, redness or yellowness. Samples exhibited pseudoplastic behavior, while viscosity was iníluenced by temperature and followed an Arrhenius-type relation. Specifíc heat was reduced with an increase in concentration and density, following a direct relation with concentration and inverse with temperature. Density was estimated by equations at ali concentrations. Thermal conductivity, calculated by equations proposed in the literature, did not result in acceptable estimates.

Petrobras / No processo de perfuração de poços de petróleo e gás o controle da pressão é uma importante tarefa. Se a pressão no interior do poço estiver abaixo da pressão de poros, um influxo da formação pode ocorrer, fenômeno denominado de kick. Se o influxo não for controlado pode ocorrer um fluxo descontrolado da formação para a superfície (blowout). As pressões de fechamento são utilizadas para calcular a pressão do reservatório, entretanto, o caráter viscoplástico do fluido de perfuração atenua as pressões medidas. No presente trabalho, é apresentada uma modelagem matemática e numérica para prever a propagação de pressão ao longo do poço durante um influxo de gás (kick). Considerou-se a compressibilidade e o comportamento viscoplástico do fluido de perfuração, modelando-o como Plástico de Bingham. O escoamento é considerado unidimensional, laminar, transiente, isotérmico e homogêneo. A solubilidade do gás no fluido de perfuração é desconsiderada e o gás é modelado como um gás real. O fluxo da formação para o interior do poço é tratado através da lei de Darcy. Os balanços de massa e de quantidade de movimento para a mistura formam um sistema de equações diferenciais parciais hiperbólicas, tendo como incógnitas a pressão e a velocidade, as quais são resolvidas pelo método das características. Consideram-se dois casos de estudo: o kick dinâmico e o kick estático. A obtenção dos campos de velocidade e pressão ao longo do poço é realizada através de um programa computacional desenvolvido em linguagem FORTRAN. Posteriormente, aferiu-se os resultados obtidos pelo modelo desenvolvido comparando-os com solução analítica e dados experimentais. Os resultados apresentam que quanto maior a tensão limite de escoamento, menor é a transmissão de pressão e que se a pressão do reservatório não for grande o suficiente, não há ganho de volume de fluido nos tanques de lama. / In the process of drilling oil and gas wells the pressure control is an important task. If the pressure inside the well is smaller than the formation pressure, an influx of the formation may occur, phenomenon denominated as kick. If the influx is not controlled there may be an uncontrolled flux from the formation to the surface (blowout). The closing pressures are used to calculate the formation pressure, however, the viscoplastic character of the drilling fluid reduces the pressures measured on the surface. This current work presents a mathematical and numerical modeling to predict the pressure propagation along the well during a gas influx (kick). The compressibility and the viscoplastic behavior of the drilling fluid were considered, modeling it as a Bingham Plastic. The flow is considered as one-dimensional, laminar, transient, isothermal and homogeneous. The solubility of the gas in the drilling fluid is disregarded and the gas is modeled as a real gas. The flux of the formation into the wellbore is treated by Darcy’s law. The balance equations of mass and momentum for the mixture result in a system of hyperbolic partial differential equations, having as unknowns the pressure and the velocity, which are solved by the method of characteristics. Two study cases were considered: the dynamic kick and the static kick. The pressure and velocity fields along the well are obtained by a computer program developed in FORTRAN language. Afterwards, the results obtained from the model developed were assessed by comparing to analytic solution and experimental data. The results present that the bigger the yield stress is, the smaller is the pressure transmission and that if the formation pressure is not high enough, there is no pit-gain in the mud tanks.