Ultra-high-speed chemiluminescence tomography of spinning-mode detonation waves

Resumo Técnico: Tomografia de Quimioluminescência de Ultra-Alta Velocidade de Ondas de Detonação em Modo de Rotação

Definição do Problema
As detonações gasosas são inerentemente multidimensionais e instáveis, frequentemente exibindo estruturas complexas, como hastes de Mach, choques incidentes e ondas transversais que se encontram em pontos triplos. Embora as teorias clássicas (Chapman–Jouguet e Zel'dovich–von Neumann–Döring) forneçam uma base para prever velocidades médias de propagação, elas não conseguem capturar totalmente o arranjo tridimensional (3D) instantâneo das ondas de choque e das estruturas reativas observadas em geometrias confinadas do mundo real. Os diagnósticos convencionais — como folhas de fuligem, sensores de pressão e imagens planares — são dimensionalmente limitados, fornecendo apenas traços de parede, históricos pontuais ou projeções integradas de linha de visão. Essas limitações dificultam a resolução da dinâmica volumétrica, particularmente para modos de propagação complexos, como detonações rotativas em tubos redondos, onde efeitos de camada limite em canais estreitos podem alterar a física. Além disso, as aplicações tomográficas existentes às detonações foram limitadas por ângulos de visualização esparsos ou pela falta de correção para refração através de elementos ópticos curvos, o que distorce os campos reconstruídos.

Metodologia
Os autores apresentam uma campanha de tomografia de quimioluminescência resolvida no tempo, projetada para reconstruir estruturas reativas 3D em misturas de etileno a 1 atm. A configuração experimental utiliza um tubo redondo acessível opticamente com uma seção de teste de safira cilíndrica (DI 66,70 mm, DE 88,90 mm).

• Sistema de Imagem: Cinco câmeras de alta velocidade (taxas de quadros de 2–5 MHz) estão posicionadas ao redor da seção de teste em ângulos específicos (0°, 40°, 90°, 130°, 200°) para capturar a quimioluminescência de banda larga (dominada pela emissão de CH*).
• Calibração: Um procedimento de calibração personalizado foi desenvolvido para abordar a refração causada pela janela de safira cilíndrica. Diferente dos métodos padrão de alvo planar, esta abordagem impõe a lei de Snell nas interfaces do cilindro. Um alvo de grade de pontos 9×9 foi imageado através do cilindro, e um algoritmo de Levenberg–Marquardt otimizou as poses das câmeras, os parâmetros do cilindro (raios, índice de refração) e as posições dos pontos para minimizar os erros de reprojeção (alcançando um erro RMS de ~0,15–0,3 px).
• Algoritmo de Reconstrução: O campo de emissão 3D é reconstruído usando uma grade cilíndrica ajustada à parede com uma base trilinear para resolver com precisão as estruturas aderidas à parede do tubo. O problema inverso é formulado como uma otimização de mínimos quadrados não negativos de grande escala. Ele minimiza um termo de perda de dados (diferença entre as projeções medidas e as previsões do modelo) regularizado por uma penalidade de Tikhonov de segunda ordem para impor suavidade e suprimir artefatos. A otimização é resolvida utilizando um algoritmo de Broyden–Fletcher–Goldfarb–Shanno de memória limitada (L-BFGS).

Principais Contribuições

1. Tomografia Corrigida por Refração: O desenvolvimento de um framework de calibração personalizado que considera a refração através de uma seção de teste óptica cilíndrica, permitos o rastreamento de raios preciso para a reconstrução tomográfica em tubos redondos.
2. Imagem Volumétrica de Alta Velocidade: A aplicação de tomografia de quimioluminescência de múltiplas visões em taxas de megahertz para resolver a evolução macroscópica de ondas de detonação, indo além das limitações de projeções 2D.
3. Extração Cinemática Quantitativa: Uma metodologia para extrair parâmetros geométricos e cinemáticos fundamentais (velocidades axiais/azimutais, ângulos de trajetória, ângulos de haste de Mach) diretamente dos campos 3D reconstruídos.

Resultados
O método foi aplicado a três casos representativos:

• Caso A (Detonação Rotativa Persistente): Em uma mistura de C₂H₄–ar–Ar próxima à estequiométrica, a reconstrução visualizou claramente uma onda de rotação de cabeça única. O campo 3D revelou uma forte frente aderida à parede (haste de Mach) e uma estrutura reativa oblíqua (onda transversal) encontrando-se em um ponto triplo. O choque incidente não foi visível, consistente com o atraso na liberação de calor atrás dele. O rastreamento da trajetória do ponto triplo resultou em uma velocidade axial de 1388 m/s (correspondendo aos dados de sonda de íons) e uma velocidade azimutal de 1468 m/s. O ângulo de trajetória derivado ($\beta_1$) foi de 46,6°, e o passo foi de 2,97 diâmetros de tubo, consistente com a literatura prévia. O ângulo da haste de Mach local ($\beta_2$) foi estimado entre 80° e 90°, indicando uma haste forte e quase normal.
• Caso B (Falha de Detonação Rotativa): Em uma mistura rica em combustível, a onda fragmentou-se. A tomografia isolou com sucesso características em forma de anel emanando da parede que estavam sobrepostas à emissão distribuída em projeções 2D. A reconstrução permitiu o rastreamento desses anéis, revelando uma taxa de expansão azimutal de ~1720 m/s.
• Caso C (Frentes Transversais Contra-Propagantes): Em uma mistura pobre, duas frentes transversais colidiram. A reconstrução tomográfica separou as estruturas pós-colisão, que pareciam convolutas em projeções 2D, revelando um "blob" axial de movimento rápido e outro de movimento lento.

Significância e Alegações
O artigo afirma que a tomografia de quimioluminescência resolvida no tempo serve como uma ferramenta não intrusiva capaz de resolver a dinâmica volumétrica de detonação que é difícil de inferir por meio de diagnósticos pontuais, planares ou de linha de visão. Ao reconstruir o campo do termo de fonte de emissão efetivo, os autores demonstram a capacidade de visualizar morfologias de ondas e determinar velocidades de onda e parâmetros geométricos precisos. Os autores observam que, embora a técnica forneça um marcador para estruturas reativas, ela não mede diretamente variáveis de estado (densidade, temperatura) ou a própria estrutura do choque; consequentemente, regiões de emissão fraca (como o rastro do choque incidente) podem permanecer indetectáveis. O trabalho estabelece uma base para benchmarks mais rigorosos para simulações e descrições de ordem reduzida de estabilidade e limites de propagação de detonação.