Investigação numérica e experimental do escoamento em válvulas de compressores herméticos

Pós-graduação em Engenharia Mecânica - FEIS

Análise numérica de condensadores do tipo arame-sobre-tubo usados em refrigeradores domésticos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise numérica de evaporadores do tipo Roll-Bond usados em refrigeradores domésticos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise numérica e experimental na determinação da potência térmica dissipada em componentes eletrônicos

Pós-graduação em Engenharia Mecânica - FEIS

Avaliação do impacto do reservatório de Três Irmãos sobre a superfície potenciométrica do aquífero livre na cidade de Pereira Barreto (SP): uma abordagem numérica e geoestatística

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Simulação de perfis nucleares de poço em formações complexas

A identificação e descrição dos caracteres litológicos de uma formação são indispensáveis à avaliação de formações complexas. Com este objetivo, tem sido sistematicamente usada a combinação de ferramentas nucleares em poços não-revestidos. Os perfis resultantes podem ser considerados como a interação entre duas fases distintas: • Fase de transporte da radiação desde a fonte até um ou mais detectores, através da formação. • Fase de detecção, que consiste na coleção da radiação, sua transformação em pulsos de corrente e, finalmente, na distribuição espectral destes pulsos. Visto que a presença do detector não afeta fortemente o resultado do transporte da radiação, cada fase pode ser simulada independentemente uma da outra, o que permite introduzir um novo tipo de modelamento que desacopla as duas fases. Neste trabalho, a resposta final é simulada combinando soluções numéricas do transporte com uma biblioteca de funções resposta do detector, para diferentes energias incidentes e para cada arranjo específico de fontes e detectores. O transporte da radiação é calculado através do algoritmo de elementos finitos (FEM), na forma de fluxo escalar 2½-D, proveniente da solução numérica da aproximação de difusão para multigrupos da equação de transporte de Boltzmann, no espaço de fase, dita aproximação P₁, onde a variável direção é expandida em termos dos polinômios ortogonais de Legendre. Isto determina a redução da dimensionalidade do problema, tornando-o mais compatível com o algoritmo FEM, onde o fluxo dependa exclusivamente da variável espacial e das propriedades físicas da formação. A função resposta do detector NaI(Tl) é obtida independentemente pelo método Monte Carlo (MC) em que a reconstrução da vida de uma partícula dentro do cristal cintilador é feita simulando, interação por interação, a posição, direção e energia das diferentes partículas, com a ajuda de números aleatórios aos quais estão associados leis de probabilidades adequadas. Os possíveis tipos de interação (Rayleigh, Efeito fotoelétrico, Compton e Produção de pares) são determinados similarmente. Completa-se a simulação quando as funções resposta do detector são convolvidas com o fluxo escalar, produzindo como resposta final, o espectro de altura de pulso do sistema modelado. Neste espectro serão selecionados conjuntos de canais denominados janelas de detecção. As taxas de contagens em cada janela apresentam dependências diferenciadas sobre a densidade eletrônica e a fitologia. Isto permite utilizar a combinação dessas janelas na determinação da densidade e do fator de absorção fotoelétrico das formações. De acordo com a metodologia desenvolvida, os perfis, tanto em modelos de camadas espessas quanto finas, puderam ser simulados. O desempenho do método foi testado em formações complexas, principalmente naquelas em que a presença de minerais de argila, feldspato e mica, produziram efeitos consideráveis capazes de perturbar a resposta final das ferramentas. Os resultados mostraram que as formações com densidade entre 1.8 e 4.0 g/cm³ e fatores de absorção fotoelétrico no intervalo de 1.5 a 5 barns/e^-, tiveram seus caracteres físicos e litológicos perfeitamente identificados. As concentrações de Potássio, Urânio e Tório, puderam ser obtidas com a introdução de um novo sistema de calibração, capaz de corrigir os efeitos devidos à influência de altas variâncias e de correlações negativas, observadas principalmente no cálculo das concentrações em massa de Urânio e Potássio. Na simulação da resposta da sonda CNL, utilizando o algoritmo de regressão polinomial de Tittle, foi verificado que, devido à resolução vertical limitada por ela apresentada, as camadas com espessuras inferiores ao espaçamento fonte - detector mais distante tiveram os valores de porosidade aparente medidos erroneamente. Isto deve-se ao fato do algoritmo de Tittle aplicar-se exclusivamente a camadas espessas. Em virtude desse erro, foi desenvolvido um método que leva em conta um fator de contribuição determinado pela área relativa de cada camada dentro da zona de máxima informação. Assim, a porosidade de cada ponto em subsuperfície pôde ser determinada convolvendo estes fatores com os índices de porosidade locais, porém supondo cada camada suficientemente espessa a fim de adequar-se ao algoritmo de Tittle. Por fim, as limitações adicionais impostas pela presença de minerais perturbadores, foram resolvidas supondo a formação como que composta por um mineral base totalmente saturada com água, sendo os componentes restantes considerados perturbações sobre este caso base. Estes resultados permitem calcular perfis sintéticos de poço, que poderão ser utilizados em esquemas de inversão com o objetivo de obter uma avaliação quantitativa mais detalhada de formações complexas.

Estabilidade de métodos numéricos LCR para a simulação de escoamento

Pós-graduação em Matematica Aplicada e Computacional - FCT

Desenvolvimento de uma metodologia numérica para escoamentos viscoelásticos não-isotérmicos

Pós-graduação em Matematica Aplicada e Computacional - FCT

Simulação tridimensional adaptativa da separação das fases de uma mistura bifásica usando a equação de Cahn-Hilliard

Simulamos a separação dos componentes de uma mistura bifásica com a equação de Cahn-Hilliard. Esta equação contém intrincados termos não lineares e derivadas de alta ordem. Além disso, a delgada região de transição entre os componentes da mistura requer muita resolução. Assim, determinar a solução numérica da equação de Cahn-Hilliard não é uma tarefa fácil, principalmente em três dimensões. Conseguimos a resolução exigida no tempo usando uma discretização semi-implícita de segunda ordem. No espaço, obtemos a precisão requerida utilizando malhas refinadas localmente com a estratégia AMR. Essas malhas se adaptam dinamicamente para recobrir a região de transição. O sistema linear proveniente da discretização é solucionado por intermédio de técnicas multinível-multigrid.

Impacto da temperatura da superfície do mar na simulação da Zona de Convergência do Atlântico Sul

Durante o verão austral, diversas regiões do Brasil são afetadas por precipitação intensa, geralmente associada à Zona de Convergência do Atlântico Sul (ZCAS). O objetivo geral deste trabalho foi investigar a influência da resolução espacial e temporal dos dados de temperatura da superfície do mar (TSM) na simulação da precipitação associada à ZCAS. Foram realizadas simulações com o modelo BRAMS (Brazilian developments on the Regional Atmospheric Modeling System) para dois eventos de ZCAS, ocorridos em 1998 (El Niño) e 1999 (La Niña). A quantidade de precipitação acumulada na parte oceânica da ZCAS foi maior nos experimentos com TSM mais quente. Índices estatísticos foram utilizados para verificação do desempenho do modelo na simulação de precipitação nas regiões que compõem a ZCAS (oceânica, costeira e amazônica), com diferentes dados de TSM. A resolução espaço-temporal dos dados de TSM influencia de forma pouco significativa na representação da ZCAS pelo modelo BRAMS. O modelo é mais eficiente em identificar a ocorrência/não ocorrência de chuva do que em localizar núcleos mais intensos e seu desempenho foi superior (inferior) na região amazônica (oceânica) da ZCAS.

Análise experimental e numérica de quedas de pressão em canais estenosados

O interesse nos estudos relacionados com o escoamento de fluidos em microcanais tem assumido um lugar de destaque na área de biomédica, uma vez que que estes estudos permitem analisar escoamentos sanguíneos evitando questões de ética que se colocam quando se trabalha fisicamente com o sangue. Numa primeira fase deste trabalho foram fabricados modelos anatómicos tridimensionais, mais concretamente canais circulares, um sem patologia e dois com ateromas de alturas diferentes, onde foram efetuadas medições de quedas de pressão. Com os resultados obtidos foi possível concluir que a presença de ateromas influencia o escoamento sanguíneo, verificando-se que as quedas de pressão aumentam com o aumento do grau da estenose. Numa segunda fase os escoamentos estudados experimentalmente foram alvo de um estudo numérico com o objetivo de comparar os resultados da simulação com os obtidos experimentalmente. As simulações numéricas foram efetuadas recorrendo a um software comercial que implementa o método dos volumes finitos e os resultados obtidos revelaram-se concordantes com os obtidos experimentalmente, o que permitiu concluir que a dinâmica de fluidos computacional é de facto uma boa estratégia a ser utilizada. Na fase final do trabalho, procedeu-se à identificação da melhor modelação matemática para caracterizar a diferença de valores obtidos nos estudos experimental e numérico. Após diversas estratégias conclui-se que a diferença de valores pode ser caracterizada por uma função polinomial.

Otimização de vazão de poços injetores em projeto de injeção de água utilizando simulação por linhas de fluxo

Until the early 90s, the simulation of fluid flow in oil reservoir basically used the numerical technique of finite differences. Since then, there was a big development in simulation technology based on streamlines, so that nowadays it is being used in several cases and it can represent the physical mechanisms that influence the fluid flow, such as compressibility, capillarity and gravitational segregation. Streamline-based flow simulation is a tool that can help enough in waterflood project management, because it provides important information not available through traditional simulation of finite differences and shows, in a direct way, the influence between injector well and producer well. This work presents the application of a methodology published in literature for optimizing water injection projects in modeling of a Brazilian Potiguar Basin reservoir that has a large number of wells. This methodology considers changes of injection well rates over time, based on information available through streamline simulation. This methodology reduces injection rates in wells of lower efficiency and increases injection rates in more efficient wells. In the proposed model, the methodology was effective. The optimized alternatives presented higher oil recovery associated with a lower water injection volume. This shows better efficiency and, consequently, reduction in costs. Considering the wide use of the water injection in oil fields, the positive outcome of the modeling is important, because it shows a case study of increasing of oil recovery achieved simply through better distribution of water injection rates

Simulação do escoamento monofásico em um estágio de uma bomba centrífuga utilizando técnicas de fluidodinâmica computacional

Oil production and exploration techniques have evolved in the last decades in order to increase fluid flows and optimize how the required equipment are used. The base functioning of Electric Submersible Pumping (ESP) lift method is the use of an electric downhole motor to move a centrifugal pump and transport the fluids to the surface. The Electric Submersible Pumping is an option that has been gaining ground among the methods of Artificial Lift due to the ability to handle a large flow of liquid in onshore and offshore environments. The performance of a well equipped with ESP systems is intrinsically related to the centrifugal pump operation. It is the pump that has the function to turn the motor power into Head. In this present work, a computer model to analyze the three-dimensional flow in a centrifugal pump used in Electric Submersible Pumping has been developed. Through the commercial program, ANSYS® CFX®, initially using water as fluid flow, the geometry and simulation parameters have been defined in order to obtain an approximation of what occurs inside the channels of the impeller and diffuser pump in terms of flow. Three different geometry conditions were initially tested to determine which is most suitable to solving the problem. After choosing the most appropriate geometry, three mesh conditions were analyzed and the obtained values were compared to the experimental characteristic curve of Head provided by the manufacturer. The results have approached the experimental curve, the simulation time and the model convergence were satisfactory if it is considered that the studied problem involves numerical analysis. After the tests with water, oil was used in the simulations. The results were compared to a methodology used in the petroleum industry to correct viscosity. In general, for models with water and oil, the results with single-phase fluids were coherent with the experimental curves and, through three-dimensional computer models, they are a preliminary evaluation for the analysis of the two-phase flow inside the channels of centrifugal pump used in ESP systems

Modelo dinâmico de simulação e otimização da força normal de reação do solo para um mecanismo de corte basal

The physical model was based on the method of Newton-Euler. The model was developed by using the scientific computer program Mathematica®. Several simulations where tried varying the progress speeds (0.69; 1.12; 1.48; 1.82 and 2.12 m s-1); soil profiles (sinoidal, ascending and descending ramp) and height of the profile (0.025 and 0.05 m) to obtain the normal force of soil reaction. After the initial simulations, the mechanism was optimized using the scientific computer program Matlab® having as criterion (function-objective) the minimization of the normal force of reaction of the profile (FN). The project variables were the lengths of the bars (L1y, L2, l3 and L4), height of the operation (L7), the initial length of the spring (Lmo) and the elastic constant of the spring (k t). The lack of robustness of the mechanism in relation to the variable height of the operation was outlined by using a spring with low rigidity and large length. The results demonstrated that the mechanism optimized showed better flotation performance in relation to the initial mechanism.