Large eddy simulation of turbulent swirl-stabilized flames using the front propagation formulation: impact of the resolved flame thickness

Este trabalho estende o modelo de combustão por formulação de propagação de frente para simulações de grandes vórtices (LES) de chamas pré-misturadas estabilizadas por turbilhão, demonstrando que a modelagem adequada da espessura da chama resolvida é crucial para capturar corretamente a dinâmica de flutuações de temperatura secundárias e a formação de bolsas de chama presas no TECFLAM.

O essencial

Imagine que você está tentando prever o comportamento de uma fogueira gigante dentro de um motor de avião. Não é uma fogueira comum; ela é alimentada por um jato de ar que gira como um redemoinho (um "vórtice"), o que ajuda a manter a chama estável e eficiente.

Este artigo de pesquisa é como um manual de instruções avançado para computadores que simulam essas chamas. Os autores, Ruochen Guo, Yunde Su e Yuewen Jiang, estão tentando resolver um problema muito específico: como fazer o computador "ver" a chama com a nitidez certa, sem cometer erros de cálculo.

Aqui está a explicação do que eles fizeram, usando analogias do dia a dia:

1. O Problema: A "Fotografia" Borrada

Quando os cientistas usam supercomputadores para simular chamas (chamado de Simulação de Grandes Vórtices ou LES), eles dividem o espaço em pequenos cubos (como pixels em uma foto). O problema é que, para economizar tempo de processamento, esses "pixels" são grandes.

• A Analogia: Imagine tentar desenhar uma linha muito fina e nítida em um quadro branco usando apenas um pincel gigante. O resultado será uma linha grossa e borrada.
• Na simulação: Quando o computador usa "pixels" grandes, a chama parece mais grossa do que realmente é. Isso causa um erro: o computador acha que a chama está se movendo sozinha (como um fantasma), quando na verdade ela deveria estar seguindo o fluxo de ar.

2. A Solução: O "Modelo de Propagação de Frente" (FPF)

Os autores usaram e melhoraram um método chamado FPF. Pense nele como um "filtro mágico" ou um "lente de aumento" para o computador.

• O que ele faz: Em vez de deixar a chama ficar borrada, o modelo força a chama a manter uma estrutura interna mais realista, mesmo que os "pixels" do computador sejam grandes. Ele garante que a química da queima (a "física" da chama) seja respeitada, mesmo na simulação.
• A Melhoria: Eles adaptaram esse modelo para situações reais, onde há perda de calor (a chama esfria nas paredes) e mistura de combustíveis. É como ajustar a receita de um bolo para que ele fique bom mesmo se a temperatura do forno variar.

3. A Descoberta Chave: A Espessura da Chama

A parte mais interessante do estudo é sobre a espessura da chama resolvida.

• A Metáfora: Imagine que a chama é uma faixa de tecido.
  • Se o computador modela essa faixa como muito grossa (porque não consegue ver os detalhes finos), ele perde a capacidade de ver o que está acontecendo nas bordas.
  • Se ele modela a faixa com a espessura correta (afinando-a quimicamente), ele consegue ver os detalhes.

Os autores descobriram que, se a "faixa" da chama ficar muito grossa na simulação, o computador erra feio em um ponto crucial: a formação de "bolsas de fogo".

4. O Fenômeno dos "Bolsões de Fogo" (Flame Pockets)

No experimento real, os cientistas viram algo curioso: além da chama principal, havia pequenos "bolsões" de fogo presos nas bordas do redemoinho de ar, criando picos de temperatura secundários.

• O que acontece na vida real: O ar girando (o vórtice) puxa e estica a ponta da chama, como se fosse um elástico. Às vezes, esse elástico se rompe e forma uma bolinha de fogo isolada que fica presa no redemoinho.
• O erro do computador (quando a chama está "gorda"): Se a simulação faz a chama parecer muito grossa, ela se afasta das bordas do redemoinho. É como se a chama fosse um casaco muito largo; o vento não consegue pegar na ponta dele para esticá-lo. Resultado: o computador não vê as bolinhas de fogo e não prevê os picos de temperatura.
• O acerto do computador (quando a chama está "fina"): Com o modelo correto (o FPF ajustado), a chama é fina o suficiente para ser "agarrada" pelo redemoinho. O computador então vê as bolinhas se formando e prevê corretamente onde o calor vai subir.

5. Conclusão: Por que isso importa?

O estudo mostra que, para projetar motores de avião mais limpos e eficientes, não basta apenas ter um computador rápido. É preciso ter o modelo matemático certo para definir a "espessura" da chama.

Se você modelar a chama como muito grossa, você perde a física real dos pequenos redemoinhos que criam bolsões de calor. Isso pode levar a motores que não funcionam como esperado ou que poluem mais.

Resumo em uma frase:
Os autores criaram um "óculos de realidade aumentada" para computadores, permitindo que eles vejam os detalhes finos de como o fogo e o ar giratório interagem, garantindo que a simulação preveja corretamente onde o calor vai aparecer, evitando erros que poderiam custar caro na engenharia de motores reais.

Resumo técnico

Título: Simulação de Grandes Vórtices (LES) de Chamas Turbulentas Estabilizadas por Vórtice usando a Formulação de Propagação de Frente: Impacto da Espessura da Chama Resolvida

Autores: Ruochen Guo, Yunde Su e Yuewen Jiang.

---

1. Problema e Contexto

A combustão pré-misturada estabilizada por vórtice é fundamental em turbinas a gás e motores aeroespaciais modernos devido à sua alta eficiência e estabilidade. No entanto, a interação complexa, não linear e multiescala entre o escoamento turbulento e a chama apresenta desafios significativos para simulações numéricas de alta fidelidade.

A Simulação de Grandes Vórtices (LES) é uma abordagem promissora, mas enfrenta o "problema de resolução da chama". Em aplicações práticas, o tamanho da malha computacional é frequentemente comparável ou maior que a espessura da chama laminar. Quando a malha é usada como o tamanho do filtro (LES implícito), a zona de reação filtrada fica sub-resolvida. Isso leva a:

• Previsões imprecisas da taxa total de reação.
• Propagação espúria (fictícia) da frente de chama filtrada.
• Dificuldade em capturar a interação correta entre a chama e estruturas vorticais, especialmente em combustores complexos com estratificação de combustível e perdas de calor.

Embora métodos como a Front Propagation Formulation (FPF) tenham sido desenvolvidos para mitigar a propagação espúria em chamas pré-misturas canônicas, sua aplicação a combustores com vórtices complexos, estratificação e efeitos não adiabáticos ainda não havia sido explorada. Além disso, o impacto da espessura da chama resolvida (determinada pelo balanço entre difusão turbulenta e afinamento químico) sobre a dinâmica da chama em LES permanece pouco estudado.

2. Metodologia

Extensão do Modelo FPF:
Os autores estenderam a formulação de propagação de frente (FPF) para LES de chamas turbulentas estratificadas e não adiabáticas. O modelo baseia-se na formulação da taxa de reação filtrada como o produto da densidade não queimada, a velocidade de consumo da chama sub-filtro e uma função delta regularizada.

• Função de Estrutura da Frente ($\psi$): Introduzida para eliminar erros de propagação em chamas sub-resolvidas e reproduzir os efeitos de "afinamento químico" (chemical steepening). Os parâmetros do modelo ($\alpha$ e $\gamma$) controlam a espessura da chama resolvida.
• Modelagem Não Adiabática: Foi incorporado um banco de dados (library) para velocidade de chama laminar ($S_L$) e espessura ($l_F$) que considera perdas de calor, baseado em simulações 1D de chamas pré-misturadas com diferentes fatores de perda de calor.
• Velocidade de Chama Sub-filtro ($S_{c,t}$): Utilizou-se o modelo de Peters-Pitsch para o fator de enrugamento sub-filtro ($\Xi_\Delta$). Para resolver inconsistências causadas pela liberação de calor na turbulência local, os autores propuseram uma melhoria: estimar o fator de enrugamento dentro da escova de chama filtrada estendendo os valores da região não queimada (onde a liberação de calor é desprezível) ao longo da normal da superfície, impondo uma condição de gradiente nulo.

Configuração Numérica e Experimental:

• Objeto de Estudo: Queimador swirl-stabilized TECFLAM (sem confinamento), com mistura metano-ar pré-misturada, número de Reynolds de 10.000 e número de swirl geométrico de 0,75.
• Condições de Contorno: Chama estratificada, perdas de calor no corpo cego (bluff body) resfriado a 353 K.
• Esquema Numérico: Diferenças finitas em malha não uniforme cilíndrica. Esquema WENO de 5ª ordem para transporte de escalares e esquema central conservador de energia de 2ª ordem para conservação de massa e momento.
• Simulações Comparativas: Foram realizadas duas simulações principais para investigar o impacto da espessura da chama:
  1. Caso Base ($\gamma = 2,5$): Parâmetros determinados fisicamente para capturar o afinamento químico.
  2. Caso de Controle ($\gamma = 1,0$): Efeitos de afinamento químico desativados (sem steepening), resultando em uma espessura de chama resolvida artificialmente maior.

3. Resultados Principais

Validação A Posteriori:
As simulações com o modelo FPF estendido ($\gamma = 2,5$) apresentaram excelente concordância com dados experimentais para:

• Perfis de velocidade média e flutuante (axial, radial e azimutal).
• Fração de mistura e temperatura.
• Captura bem-sucedida da zona de recirculação central e das camadas de cisalhamento interna e externa.

Mecanismo dos Picos Secundários de Temperatura:
A simulação revelou que estruturas vorticais de grande escala na camada de cisalhamento externa esticam a frente de chama, causando o enrolamento (roll-up) e a formação de bolsas de chama isoladas (flame pockets) ou penínsulas que se destacam da superfície principal.

• Essas bolsas ficam presas na região de baixa velocidade da camada de cisalhamento externa.
• A presença dessas bolsas é a causa física dos picos secundários de temperatura (médios e flutuantes) observados experimentalmente na camada de cisalhamento externa.

Impacto Crítico da Espessura da Chama Resolvida:
A comparação entre os casos $\gamma = 2,5$ e $\gamma = 1,0$ demonstrou que a espessura da chama resolvida é um fator determinante:

• Caso $\gamma = 1,0$ (Sem afinamento químico): A espessura da chama resolvida foi superestimada (cerca de 3 vezes maior que o caso base na região próxima à saída). Isso levou a uma subestimação da área superficial da chama e da taxa de consumo total.
• Consequência: A escova de chama "desacoplou" (detachou-se) da camada de cisalhamento externa. Sem o contato correto entre a frente de chama e os vórtices externos, a formação de bolsas de chama não ocorreu.
• Falha na Predição: O caso com espessura superestimada falhou em capturar os picos secundários de temperatura, pois o mecanismo de formação das bolsas de chama foi suprimido.

4. Contribuições Chave

1. Extensão do Modelo FPF: Adaptação bem-sucedida do método FPF para combustão pré-misturada turbulenta complexa, incluindo estratificação de combustível e efeitos não adiabáticos (perdas de calor).
2. Correção de Inconsistência no Modelo de Velocidade: Desenvolvimento de uma abordagem para estimar o fator de enrugamento sub-filtro que resolve a inconsistência entre os parâmetros turbulentos definidos no modelo e os parâmetros locais afetados pela liberação de calor.
3. Descoberta do Mecanismo de Formação de Bolsas: Identificação de que a interação entre vórtices de grande escala na camada de cisalhamento externa e a frente de chama é o mecanismo responsável pelos picos secundários de temperatura no queimador TECFLAM.
4. Importância do Afinamento Químico: Demonstração de que, em LES de chamas com vórtices, modelar corretamente os efeitos de afinamento químico (para manter a espessura da chama resolvida próxima à física real) é essencial. Ignorar esse efeito leva a erros graves na previsão da dinâmica da chama e na interação chama-vórtice.

5. Significado e Conclusão

Este trabalho destaca que, em simulações LES de combustores com vórtices, a espessura da chama resolvida não é apenas um detalhe numérico, mas uma variável física crítica. Se a espessura da chama for superestimada (devido à falta de modelagem do afinamento químico), a taxa de consumo de combustível é subestimada, alterando a posição relativa da frente de chama em relação às estruturas turbulentas.

Isso impede a captura de fenômenos físicos importantes, como a formação e aprisionamento de bolsas de chama, resultando em previsões incorretas de perfis de temperatura e estabilidade da chama. O estudo valida que a formulação FPF, quando devidamente estendida e calibrada, é uma ferramenta robusta para simular combustão turbulenta complexa com alta fidelidade, desde que os efeitos de afinamento químico sejam adequadamente representados.