Simulação numérica da combustão de metano em câmara cilíndrica utilizando o CFX 10.0

As ferramentas de simulação numérica têm evoluído consideravelmente nos últimos vinte anos e aproximam-se cada vez mais do dia a dia das soluções de Engenharia, fornecendo diversas funcionalidades que o Engenheiro pode implementar para a previsão de fenômenos, antes exclusivamente compreendidos experimentalmente. Mesmo sendo muito úteis e razoavelmente compreensíveis para os usuários, estas ferramentas não compreendem toda a gama de adaptações necessárias para resolver casos específicos, e podem trazer prejuízos à compreensão do problema e até uma solução errônea do caso em questão. O objetivo do presente trabalho é utilizar uma ferramenta comercial de simulação numérica para a solução de dois casos de escoamentos reativos. Visa também discutir criticamente a solução obtida e apresentar uma maneira de aplicar os conhecimentos adquiridos na simulação de uma câmara de combustão real. O primeiro estudo de caso é uma câmara cilíndrica experimentada e apresentada em 1994 por Garreton e Simonin. Trata-se de um caso semelhante a fornalhas de aquecimento d’água e em processos com caldeiras fumotubulares. O queimador é de jatos concêntricos utilizando gás natural e ar atmosférico como reagentes. Uma análise das soluções obtidas com diferentes modelos de turbulência é realizada e discutida para uma segunda abordagem, onde a influência das condições de contorno do domínio e dos modelos de radiação são abordadas. O segundo estudo de caso é a câmara de combustão experimentada por Meier, 2000. A geometria e condições do caso assemelham-se a de combustores de turbinas a gás, incineradores e fornos industriais. O caso apresenta um escoamento muito turbulento e instável, de gás natural e ar atmosférico. Uma discussão sobre a metodologia de solução deste problema é feita para explorar os pontos comuns a soluções de combustores industriais, como por exemplo, o de um forno secador de minério rotativo. Os diferentes resultados numéricos obtidos nos dois casos mostram que todos os modelos geram resultados pelo menos qualitativamente coerentes, mas apresentam erros numéricos razoáveis. A partir da discussão dos resultados apresenta-se uma proposta de modelagem para um caso prático envolvendo uma geometria mais complexa, e apontam-se os pontos de risco para análise dos resultados.Numerical simulation tools have reasonably evolved in the last twenty years and are getting increasingly closer to routine Engineering solutions, providing features that the Engineer can use to predict phenomena that were previously only experimentally understood. Despite being useful and easy-to-understand, these tools do not include all adaptations required to solve specific cases, leading to misunderstanding of the problem and even to wrong solutions. The present study aimed at using a commercial CFD code for the solution of two reacting flow cases, critically discussing the solutions and presenting an application of acquired the knowledge in simulations of real combustion chambers. The first case study was a cylindrical chamber, tested and presented in 1994 by Garreton and Simonin, which is similar to water-heating boilers and parallel tube furnaces. The burner uses concentric jets, with natural gas and atmospheric air as reactants. An analysis of the obtained solutions applying different turbulence models was carried out, and the solutions were discussed in a second approach on boundary conditions and radiation transfer models. The second case study was the combustion chamber tested by Meier, 2000. Its geometry and boundary conditions are similar to those of gas turbine combustors, incinerators, and industrial furnaces. The chamber presents a strongly turbulent and unstable flow of natural gas and atmospheric air. The methodology to solve this problem was discussed, considering issues also found in solutions for industrial burners, such as those used in mineral rotatory dryers. The different numerical results obtained in both cases show that all models provide at least qualitatively coherent results, but present numerical errors. Based on the discussion of the results, a proposal for modeling a practical case involving a much more complex geometry is presented, and risk issues for further result analysis are pointed out