Estudo do nitrogênio na combustão do gás natural

Study of nitrogen in the combustion of natural gas

natural gas; combustion; nitrogen

ARTIGO

Estudo do nitrogênio na combustão do gás natural

Study of nitrogen in the combustion of natural gas

Leonardo Paes RangelI, ^*; Pedro Maurício Büchler^II

^IEscola Politécnica, Universidade de Pernambuco, Rua Benfica, 455, 50750-410 Recife - PE

^IIDepartamento de Engenharia Química, Escola Politécnica, Universidade de São Paulo, Av. Prof. Luciano Gualberto, 380, Trav. 3, 05508-900 São Paulo - SP

ABSTRACT

Nitrogen content in natural gas was studied in experimental and computational investigations to identify its influence on the emission level of exhaust gases from combustion facilities. Changes in natural gas composition with different N₂ concentrations may result from introducing a new source gas into the system. An industrial burner fired at 75 kW, housed in a laboratory-scale furnace, was employed for runs where the natural gas/N₂ proportion was varied. The exhaust and in-furnace measurements of temperature and gas concentrations were performed for different combustion scenarios, varying N₂ content from 1-10 %v. Results have shown that the contamination of natural gas with nitrogen reduced the peak flame temperature, the concentration of unstable species, the NO_X emission level and the heat transfer rate to the furnace walls, resulting from the recombination reactions.

Keywords: natural gas; combustion; nitrogen.

INTRODUÇÃO

Até 1998 o gás natural teve uma participação coajuvante na matriz energética brasileira de apenas 2,7%. O aumento do seu consumo, através da importação de países vizinhos (Argentina e Bolívia), está causando uma diversificação na sua matriz, pois não existem, até o presente momento, grandes reservas de gás no país. De acordo com a projeção de dados de consumo, a participação do gás natural atual (2004) é de 9% e será cerca de 12% até 2010¹.

Após a inclusão das novas diretrizes energéticas no Brasil, o gás natural será utilizado também para geração de energia elétrica com 49 novas termelétricas sendo construídas em todo o país. Quase todas deverão ser operadas à gás natural. Com a finalização da construção do Gasbol (gasoduto entre a Bolívia e o Brasil, com comprimento de 3150 km) em 1999, o Brasil interligou os estados do Centro-Oeste, Sudeste e Sul, implementando um passo decisivo no suprimento de energia para as regiões mais industrializadas do país.

Ciente do consumo crescente do gás natural no mundo, o Brasil busca se beneficiar de uma alternativa mais eficiente para o processo de combustão industrial. Assim, grandes volumes de gás natural serão empregados em diferentes segmentos da indústria, como metalúrgica, cerâmica e têxtil, onde restrições ambientais têm aumentado particularmente nos últimos anos.

Tanto do ponto de vista energético como ambiental, as tecnologias relacionadas ao gás natural apresentam maior eficiência, principalmente quando comparadas com outros combustívies fósseis. Assim, é relevante a atenção sobre as técnicas de modificação no sistema de combustão, como o uso de queimadores de baixa emissão de NO_X.

O monóxido de nitrogênio (NO) e o dióxido de nitrogênio (NO₂) são os constituintes do NO_X. Este faz parte das reações de combustão de combustíveis fósseis e contribui para a poluição atmosférica. O NO_X foi identificado como o principal agente da disposição pluviométrica ácida (chuva ácida), que está relacionada com vários danos ambientais e à saúde humana, incluindo acidificação de lagos e correntes, impactos em solos de florestas e a corrosão acelerada de edifícios^2,3.

O gás natural é uma mistura de vários gases combustíveis, mas algumas substâncias inertes, como CO₂ e N₂, estão sempre presentes. De forma geral, o gás natural apresenta uma maior proporção de CH₄ e uma menor proporção composta por outros hidrocarbonetos, como alcanos e alcenos. Assim, o gás natural não possui uma composição química fixa e a proporção de hidrocarbonetos e gases não-hidrocarbonetos pode variar amplamente, dependendo de sua origem. O gás natural boliviano, por ex., contém cerca de 1%v de N₂, que ainda pode depender da contribuição de cada campo de extração na composição do gás vendido. Outras fontes apresentam quantidades elevadas de N₂, como o campo de Groningen (Holanda) com 14%v⁴.

Este trabalho apresenta o efeito da variação de concentração de N₂ na composição do gás natural, investigado em relação aos níveis de emissão de NO_X. Na prática, a concentração de N₂ no gás natural pode afetar a combustão através dos seus vários mecanismos de reação. Os resultados desse trabalho permitirão delinear uma nova metodologia na utilização de combustíveis gasosos com diferentes concentrações de N₂.

O trabalho foi dividido em duas partes: investigações teóricas e experimentais. As investigações experimentais foram conduzidas em uma fornalha de bancada. Um queimador comercial para gás pré-aquecido com capacidade térmica de 75 kW foi utilizado, permitindo a formação de chamas com temperaturas de até 1800 K. A fundamentação das investigações teóricas foi estudada através de um modelo computacional do queimador em um código comercial de CFD ("Computational Fluid Dynamics") com pré-processador para cálculo da emissão de NO_X⁵.

PARTE EXPERIMENTAL

Os equipamentos utilizados na etapa experimental foram uma fornalha de bancada, queimador à gás natural com ar pré-aquecido opcional, sistema de amostragem para gases de exaustão intra-fornalha, conjunto de analisadores de gases (NO_X, CO₂, CO e O₂) e unidade computadorizada para aquisição de dados.

O objeto do experimento foi um queimador comercial de gás natural com ar estagiado e pré-aquecido, com capacidade térmica máxima de 100 kW, projetado e fabricado por Stordy Combustion Engineering Ltd. (Inglaterra). A estratégia de queima foi configurada com pré-mistura ar/combustível no tubo externo e ar estagiado em alta velocidade através de um tubo de carbeto de silício central (Figura 1).

A chama foi considerada estável sob condições de operação com deficiência de ar. O ar estagiado foi injetado na direção do eixo central, próximo a saída do queimador. Assim, a combustão foi operada na zona rica de combustível e concluída por um jato central de ar, a 10 cm da câmara de combustão. Neste trabalho foi utilizado um pré-aquecedor elétrico de ar, pois a fornalha não contava com um sistema de recirculação de gás. As condições de operação típicas do queimador são mostradas na Tabela 1.

Uma fornalha experimental (3,00 m de comprimento, 0,60 m de largura e 0,90 m de altura) foi revestida internamente com uma cerâmica refratária de aquecimento rápido e com uma camada de 0,025 m de fibras cerâmicas do tipo Triton. A fornalha foi cimentada externamente com uma camada de tijolos tipo MPK 110. Sua camada externa acomoda ainda 0,03 m de espessura de fibras minerais e uma lâmina de aço.

Os orifícios de amostragem foram dispostos a diferentes distâncias (0,15; 0,45; 0,75; 1,05; 1,50; 2,10 e 2,70 m) da parede interna da fornalha, considerando o lado da entrada do queimador. Sete furos passantes (0,05 m de diâmetro) foram alocados na linha de simetria horizontal de uma das paredes laterais da fornalha, para coleta de amostras através de uma sonda isocinética de 1,50 m de comprimento. Do outro lado da fornalha, um furo de 0,19 m de diâmetro interno foi projetado para visualização da chama.

A zona de resfriamento inferior à água era de dimensões 0,60 x 0,18 x 1,80 m. A pressão interna dos gases da fornalha foi monitorada durante a operação do queimador através de um manômetro acoplado a um sinal digital. As temperaturas da fornalha foram monitoradas através de dois termopares tipo K, inseridos na parede interna da fornalha a uma distância de 0,90 e 2,10 m ao lado do queimador.

COMBUSTÃO TURBULENTA EM COMPUTAÇÃO FLUIDODINÂMICA

O modelo em computação fluidodinâmica (CFD) gerencia a entrada de dados do problema e adapta as sub-rotinas do módulo solver. As atividades do usuário no estágio de pré-processamento envolvem definição da geometria da região de interesse (domínio computacional mostrado na Figura 2), definição das propriedades do fluido, especificação das condições de contorno apropriadas e geração da malha para o sistema queimador-fornalha (Figura 3).

O código comercial de computação fluidodinâmica FLUENT™ resolve equações diferenciais parciais para conservação de massa, momento, energia e espécies químicas. De forma geral, essas equações podem ser mostradas por uma variável de conservação F, como mostrado na Equação 1

onde r é densidade [kg m^-3], u_i - velocidade da espécie química i, t - tempo [s], x - localização da espécie i no fluxo [m], GF - gradiente de F no termo difusivo e S_F - termo de origem.

A fase de solução envolve a aproximação inicial das variáveis desconhecidas, ou discretização, onde a substituição de aproximações nas equações algébricas por meio da técnica computacional de interações gerou resultados com precisões cada vez maiores⁶.

O código comercial FLUENT™ permite o cálculo das espécies químicas em um módulo principal e um pós-processador trata os resultados gerados, de modo a calcular as concentrações dos poluentes sem interferir nos parâmetros iniciais do modelo. Este é o caso da formação de NO_X, que é acionada após o modelo realizar todos os cálculos para as variáveis intrínsicas (velocidade, temperatura, espécies químicas).

Como o fluxo do sistema de combustão é turbulento, ocorrem flutuações temporais na temperatura e nas concentrações das espécies químicas, que possuem forte influência nas propriedades gerais da chama⁷.

A principal função do modelo de turbulência é fornecer uma expressão matemática para cálculo das propriedades do fluido em termos de "quantidade de fluxo". Para que os modelos de turbulência sejam úteis na abordagem de um código CFD, devem ter larga aplicabilidade, além de serem precisos e econômicos em termos computacionais. Para o presente trabalho foi usado o modelo de turbulência k-e, devido a suas características de precisão e baixo tempo computacional⁸. O modelo k-e representa-se através de duas equações diferenciais, uma para a energia cinética (k) [m² s^-2] e outra para a taxa de dissipação da energia (e)[m² s^-3]. Os vetores de Reynolds stress são assumidos como proporcionais ao gradiente de velocidade principal, com uma constante de proporcionalidade definida como viscosidade turbulenta (µ_t).

A intensidade turbulenta (I) e o comprimento característico da chama () são usados nas condições iniciais do modelo desenvolvido no código FLUENT™ para cálculo de e pela Equação 2. A taxa de dissipação de energia e é baseada na escala do comprimento turbulento relacionando-a com o tamanho dos grandes turbilhões que contém energia, presentes no fluxo turbulento. O valor inicial de k é calculado pela Equação 3.

onde é a magnitude da velocidade principal inicial [m s^-1] e a C_µ é a constante de turbulência que vale 0,092.

MODELAGEM DAS REAÇÕES QUÍMICAS

O código FLUENT™ permite a escolha entre três abordagens de mecanismos de reação química: formulação geral de taxa finita, frente pré-misturada e fração de mistura utilizando uma função probabilística de densidade (PDF - "probabilistic density function").

O conceito de fração de mistura foi utilizado para uma função PDF onde as equações de transporte das espécies químicas não são resolvidas individualmente, mas através de uma única equação, com variável escalar conservativa denominada fração de mistura (f). A função PDF considera os efeitos da turbulência através do modelo k-e e as espécies químicas em equilíbrio⁵. Através da fração de mistura no modelo, é possível o cálculo de espécies intermediárias, efeitos de dissociação e uma relação entre efeitos turbulentos e reações químicas. A abordagem é computacionalmente eficiente, pois não requer uma solução para um grande número de equações de transporte (uma para cada espécie).

Os mecanismos de reações químicas para combustíveis hidrocarbonetos são complexos e a maioria das espécies quimicamente reativas não possuem suas taxas de reação consideradas na análise de sistemas de combustão. O mecanismo resumido (com quatro etapas), aplicado no modelo deste trabalho permitiu uma razoável análise do comportamento da chama através das reações (R1) a (R4).

MODELO QUÍMICO PARA FORMAÇÃO DE NO_X

O NO pode ser formado através de três mecanismos químicos distintos: térmico, ativo e N₂O-intermediário, como mostrado na Tabela 2. A formação de NO pelo mecanismo N₂O-intermediário é importante apenas para condições de combustão pobre em combustível (f < 0,8) e sob temperaturas menores que 1200 K⁹. No presente trabalho esse mecanismo não apresentou importância, pois as simulações foram realizadas com altas temperaturas de chama e única condição de excesso de ar (f = 0,91). A taxa líquida da formação de NO térmico, principal contribuinte quando na combustão de gás natural, foi calculada pela Equação 4, onde as concentrações das espécies químicas estão em [mol m^-3].

Verificou-se que a taxa de formação do NO térmico é significante apenas para temperaturas maiores que 1000 K, pois a reação com o N₂ do ar requer a quebra da ligação tripla a uma energia de dissociação de 941 kJ mol^-110.

A Equação 5 indica que a formação de NO varia com a concentração de O₂ e, principalmente, com a temperatura. Bartok e Sarofim¹¹ verificaram que a taxa de produção de NO aumenta em duas vezes a cada aumento de temperatura de 90 K após 2200 K.

Para a resolução da Equação 5)foi necessário se conhecer as concentrações das espécies químicas O₂, N₂, O e OH. As concentrações de O₂ e N₂ foram medidas no sistema de combustão, porém as espécies instáveis O e OH foram apenas estimadas quando se assumiu um equilíbrio parcial, através das reações (R5) a (R8)⁷.

O NO ativo é formado pelo mecanismo de radicais derivados do combustível que reagem com o N₂ do ar. De uma forma geral, a contribuição do NO ativo para o NO total é pequena para queimadores estacionários, contudo a sua concentração é significante nos turbilhões das chamas ricas, turbulentas e não pré-misturadas¹².

METODOLOGIA PARA PESQUISA COMPUTACIONAL

Investigou-se o efeito do ar pré-aquecido de combustão na temperatura da chama e, conseqüentemente, na emissão de NO_X para o queimador Stordy. O efeito da variação da concentração de N₂ na composição do gás natural também foi investigado.

Para comparar efetivamente os resultados dos diferentes casos foram mantidos constantes o fluxo de ar e o fluxo de combustível, pois na prática, usuários de gás natural não podem variar as condições de operação com base na composição do gás por falta desta informação.

A temperatura do ar de combustão foi variada entre 300 e 1000 K considerando no modelo a densidade do ar. Na prática, os equipamentos são ajustados para temperaturas variáveis, pois utilizam o calor remanescente dos gases de combustão para aquecer o ar de entrada. Duas composições de gás natural foram utilizadas para investigar o efeito da concentração de N₂ na emissão de NO_X. As composições dos casos investigados são mostradas na Tabela 3.

Os casos 1 e 4 foram realizados com a temperatura do ar a 300 K. Nos casos 2 e 5 a temperatura de pré-aquecimento foi fixada a 600 K e nos casos 3 e 6, a 1000 K. Os fluxos de gás natural e do ar variaram de acordo com a composição apresentada na Tabela 3, onde o ar foi mantido fixo com 10%v de excesso e sofreu apenas variações de temperatura. A velocidade do gás natural variou conforme sua composição, contudo essa variação foi pequena quando comparada com a velocidade máxima da mistura.

RESULTADOS E DISCUSSÃO

A principal diferença entre os valores medidos para a temperatura e os valores calculados pela simulação em CFD para o sistema fornalha-queimador foi a exclusão da base resfriada à água no modelo computacional. Para se reduzir o número total de células na malha computacional, desprezou-se a base da fornalha e utilizou-se simetricamente apenas um quarto do domínio.

O perfil de temperatura foi calculado ao longo da fornalha, sob condições de operação apresentadas na Tabela 3 (caso 1) e mostrado na Figura 4. A figura mostra uma temperatura máxima calculada de 1921 K para um ar de combustão a 300 K. Para este ar pré-aquecido a 600 K obteve-se 2169 K e a 1000 K, 2293 K (limite máximo da temperatura de ar pré-aquecido para esse tipo de queimador).

Resultados similares foram obtidos quando o queimador foi operado com uma composição do gás natural com 10%v N₂. A temperatura máxima calculada foi de 1904, 2092 e 2187 K para as temperaturas do ar de combustão a 300, 600 e 1000 K, respectivamente.

O NO ativo e térmico foram calculados e mostrados para o caso 1 (Figuras 5 e 6). O NO combustível não foi considerado no modelo pelo fato do gás natural não conter "nitrogênio orgânico", ou seja, átomos de nitrogênio em sua estrutura molecular. Assim, o NO calculado através dos mecanismos ativo e térmico representa o NO_X, pois o modelo empregado é limitado por não considerar o mecanismo de formação de NO₂, formado principalmente em zonas de baixas temperaturas.

As concentrações de NO ativo foram calculadas quando o queimador foi operado sob as condições apresentadas na Tabela 4 para o caso 1. A concentração de NO ativo calculada foi sempre inferior à concentração de NO térmico. O NO ativo foi formado na zona rica em combustível, embora altas concentrações foram também encontradas na região da chama.

Para os casos listados na Tabela 4, o valor máximo para a concentração de NO ativo (27 ppm para o caso 4) foi encontrado próximo ao queimador, na zona de reação. Isso foi devido à existência de uma alta concentração de radicais derivados do combustível, como CH e CH₂, presentes no mecanismo de reação do NO ativo. O NO térmico apresentou baixa variação para os casos com 10%v de N₂ no gás natural. Devido à perda térmica com o aumento de N₂, os valores de emissão de NO para os casos 4, 5 e 6 foram menores que para os casos com 1%v de N₂ (1, 2 e 3).

Os resultados das simulações para os casos da Tabela 4 apresentaram níveis de diluição diferentes nos gases de exaustão. As concentrações de NO foram corrigidas para 0%v de O₂ em base seca, pois houve necessidade de se considerar a remoção do efeito causado pela diluição dos gases, permitindo uma comparação precisa entre os casos.

A Figura 7 confirma que o aumento do NO térmico é dependente da temperatura de pré-aquecimento do ar de combustão, independentemente do percentual de N₂ no gás natural.

O nível de emissão de NO para esse tipo de queimador produziu 99 ppm quando o queimador foi operado com gás natural a 1%v N₂ e à temperatura do ar na entrada de 300 K. Contudo, o nível de emissão de NO térmico aumentou para 287 ppm quando o ar foi pré-aquecido a 1000 K, superando os níveis de emissão máximos estabelecidos para o gás natural pela legislação ambiental. Resultados similares foram obtidos quando o queimador foi operado com gás natural com 10%v N₂. O nível de emissão do NO foi de 82 ppm para uma condição com o ar de combustão a 300 K, aumentando para 233 ppm com ar pré-aquecido a 1000 K. O efeito do aumento da concentração de N₂ no gás natural foi verificado comparando-se os níveis de emissão do NO térmico, NO ativo e NO total dos casos 1, 2 e 3 com os casos 4, 5 e 6.

O NO total emitido para este tipo de queimador foi reduzido em cerca de 21% quando a concentração de N₂ no gás natural aumentou de 1 para 10%v quando utilizado ar de combustão à temperatura de 300 K, e em cerca de 23% para uma temperatura do ar de combustão a 1000 K. Contudo, a contribuição do NO ativo aumentou de 5 para 26% quando a concentração de N₂ no gás natural aumentou de 1 para 10%v para o ar de combustão à temperatura a 300 K e de 2 para 12% para ar pré-aquecido à temperatura de 1000 K. Com o aumento da quantidade de N₂ no gás natural, a temperatura do gás de exaustão foi reduzida além da zona de reação. Assim, uma chama mais curta foi obtida para o caso com maior concentração de N₂ no gás natural e, conseqüentemente, maior destaque ao mecanismo ativo com espécies intermediárias.

A Figura 8 mostra os efeitos do ar pré-aquecido de combustão e da concentração de N₂ no gás natural em relação ao NO térmico para o queimador utilizado. O gráfico mostra que, a temperaturas de pré-aquecimento menores que 1000 K, a temperatura do ar pré-aquecido apresentou um efeito mínimo na emissão do NO térmico quando comparado com o seu efeito a temperaturas de pré-aquecimento mais elevadas que 1000 K. Assim, uma atenção especial deve ser dada à emissão de NO quando o ar de combustão for pré-aquecido a uma temperatura próxima de 1000 K ou eventualmente maior.

A formação do NO térmico tem grande dependência com a temperatura e isso foi verificado através do cálculo da taxa de formação de NO térmico pela Equação 4. A taxa de formação do NO térmico para o caso 3 foi a mais alta calculada, 1,33 10^-4 kg m^-3 s^-1, com base na temperatura máxima dentro da fornalha.

O pré-aquecimento do ar de combustão é um dos modos mais tradicionais para se produzir chamas a altas temperaturas, mas os níveis de emissão de NO devem ser monitorados na aplicação desta técnica.

Para os experimentos realizados com o queimador Stordy, o ar de combustão foi introduzido na zona de combustão em dois estágios (externo principal e central secundário) para reduzir as emissões de NO_X.

Os resultados obtidos pelas simulações computacionais para os cálculos do NO térmico e NO ativo mostraram que a contribuição do NO ativo foi menor que a do NO térmico quando o ar de combustão foi pré-aquecido. Contudo, a contribuição do NO ativo aumentou quando o queimador foi operado com maior percentual de N₂ do gás natural.

Os resultados de CFD confirmaram que a oxidação do nitrogênio atmosférico (N₂) pelo mecanismo térmico foi a principal rota para emissão de NO_X para este tipo de queimador.

A Figura 9 mostra que o modelo computacional subestimou a quantidade de NO_X produzido dentro da fornalha, que foi medido através de uma sonda isocinética em diferentes distâncias em relação à posição do queimador. Contudo, os perfis de concentração dos modelos computacional e experimental mostraram boa correlação sobre o efeito da variação da concentração de N₂ no gás natural para temperatura e NO_X apresentadas nas Figuras 7 e 8.

Neste trabalho foram encontradas pequenas discrepâncias entre os resultados experimentais e computacionais, permitindo certificar a validação do modelo CFD empregado. Os resultados obtidos do modelo CFD mostraram que o NO ativo teve uma contribuição significante no NO total emitido para queimadores de ar estagiado. A importância do NO ativo diminuiu com o aumento da temperatura do ar de combustão.

Os resultados de CFD também mostraram que, ao se utilizar temperaturas de pré-aquecimento do ar de combustão próximas a 1000 K, deve-se estimar altos níveis de emissões de NO_X para o processo.

CONCLUSÕES

Investigações teóricas usando CFD com uma integração de dados experimentais mostraram que o aumento da concentração de N₂ no gás natural provocou diversos efeitos em seu processo de combustão. Entre eles, redução da temperatura máxima da chama, menor emissão total de NO_X, redução na concentração máxima de O atômico e OH na zona de combustão, redução da transferência de calor da chama para a parede da fornalha causada pela perda de energia, absorvida pelas reações de recombinação das espécies instáveis O, H, OH e HO₂. Assim, do ponto de vista ambiental, a alta concentração de N₂ no gás natural é positiva, pois a emissão de NO_X será sempre menor. A menor capacidade calorífica do combustível com maior teor de nitrogênio diminui a temperatura de pico da chama, reduz as concentrações das espécies O atômico e OH e, conseqüentemente, a formação do NO térmico, maior responsável pelo NO total.

Os resultados experimentais mostraram que a emissão de NO_X foi reduzida em 21% quando a concentração de N₂ no gás natural aumentou de 1 para 10%v para um ar de combustão a temperatura de 300 K. O aumento da concentração de N₂ de 1 para 10 %v no ar de combustão pré-aquecido à temperaturas de 600 e 1000 K (investigadas computacionalmente) produziu uma redução de 22 e 23% na emissão de NO_X. A desvantagem associada ao aumento da concentração de N₂ no gás natural foi a redução da transferência de calor da chama para a parede da fornalha. Essa conclusão respalda o aspecto econômico, pois diminui a eficiência do combustível. Na prática, o consumidor de gás precisará queimar maior quantidade de combustível para obter a mesma eficiência térmica quando o gás natural apresentar uma concentração de nitrogênio mais elevada.

Investigações experimentais e modelagem com CFD confirmaram que a formação do NO térmico é fortemente dependente da temperatura de chama. Com a modelagem CFD foi possível relacionar o ar de combustão pré-aquecido, a temperatura de chama e o NO térmico formado. Os resultados mostraram que o efeito do aumento do ar pré-aquecido nas emissões de NO_X foi mais significativo quando o ar foi pré-aquecido a temperaturas próximas a 1000 K. Assim, para sistemas de combustão que demandam alta eficiência térmica com temperaturas de chama elevadas no interior da fornalha maiores que 2000 K, técnicas adicionais de redução de NO_X devem ser utilizadas, como por ex., a recirculação de gases de exaustão (FGR).

A modelagem em CFD indicou o mecanismo térmico como a principal rota para formação de NO sob temperaturas de chama maiores que 2000 K. A modelagem em CFD mostrou que a contribuição máxima do NO ativo para a formação total de NO_X, calculada para a exaustão, foi de 26% com o ar de combustão à temperatura de 300 K e concentração de N₂ no gás natural de 10%v (caso 4). Embora a concentração absoluta de NO ativo tenha sido aumentada proporcionalmente com pré-aquecimento do ar de combustão, a sua contribuição para a emissão total de NO_X no caso do sistema queimador-fornalha investigado foi reduzida com o aumento do ar pré-aquecido de combustão. Assim, constatou-se que o NO ativo não pode ser negligenciado nos sistemas de combustão, pois sua contribuição é significativa nas emissões totais de NO_X para queimadores de ar estagiado. Contudo, a importância do NO ativo diminui quando a temperatura do ar de combustão aumenta. Isto se deve à maior dependência do NO_X total com o mecanismo térmico, onde a temperatura assume o papel de agente principal.

Uma comparação entre os resultados experimentais e o modelo em CFD mostrou boa concordância entre as duas fontes de dados. Em geral, modelos em CFD para engenharia de combustão devem ser validados com resultados experimentais, pois a modelagem é uma ferramenta útil para a pesquisa aplicada. Seu uso é recomendado para pesquisas com sistemas de difícil controle, com substancial redução no custo e no tempo de execução.

AGRADECIMENTOS

À CAPES pela bolsa de doutoramento de L. P. Rangel, ao Departamento de Combustíveis e Energia da Universidade de Leeds (Reino Unido) pelo suporte laboratorial e computacional na realização deste trabalho.

Recebido em 7/6/04; aceito em 28/3/05; publicado na web em 8/9/05

• 1. Ministério de Minas e Energia do Brasil (MME); Balanço Energético Doméstico, 2004.
• 2. Agênica de Proteção Ambiental dos Estados Unidos (EPA-US); Nitrogen Oxide Control of Stationary Combustion Sources, 2000.
• 3. Missaghi, M.; Pourkashanian, M.; Williams, A.; Yap, L.; The Prediction of NO Emission From an Industrial Burner, American Flame Day: San Francisco, CA, 1991.
• 4. Geerssen, T. M.; Physical Properties of Natural Gases, N. V. Nederlese Gasunie, 1988.
• 5. Fluent User's Guide; Fluent Inc.; Lebanon - NH, 1998.
• 6. Versteeg, H. K.; Malalasekera, W.; An Introduction to Computational Fluid Dynamics, Longman Group, 1995, p. 3.
• 7. Yap, L.; Pourkashanian, M.; Williams, A.; Modelling of Natural Gas Combustion with Oxygen e Preheated Air, BFRC Flame Days, 1994.
• 8. Launder, B. E.; Spalding, D. B.; Lectures in Mathematical Models of Turbulence, Academic Press; London, 1972.
• 9. Turns, S. R.; An Introduction to Combustion: Concepts and Applications, Ed. McGraw-Hill, 1996, p. 199.
• 10. Tomeczek, J.; Goral, J.; Gradon, B.; 25^th Symposium International on Combustion, The Combustion Institute, 1994.
• 11. Bartok, W.; Sarofim, A. F.; Fossil Fuel Combustion, John Wiley & Sons Inc., 1989.
• 12. Miller, J. A.; Bowman, C. T.; Progress on Engineering Combustion Science, 1989, vol. 15, p. 287.

Datas de Publicação

• Publicação nesta coleção
  12 Dez 2005
• Data do Fascículo
  Dez 2005

Histórico

• Aceito
  28 Mar 2005
• Recebido
  07 Jun 2004