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

Resumen Técnico: Tomografía de Quimioluminiscencia de Ultra Alta Velocidad de Ondas de Detonación en Modo de Rotación

Planteamiento del Problema
Las detonaciones gaseosas son inherentemente multidimensionales e inestables, y suelen exhibir estructuras complejas como tallos de Mach, choques incidentes y ondas transversales que convergen en puntos triples. Aunque las teorías clásicas (Chapman–Jouguet y Zel'dovich–von Neumann–Döring) proporcionan una base para predecir las velocidades medias de propagación, no pueden capturar completamente la disposición tridimensional (3D) instantánea de las ondas de choque y las estructuras reactivas observadas en geometrías confinadas del mundo real. Los diagnósticos convencionales —como las láminas de hollín, los sensores de presión y las imágenes planas— están limitados dimensionalmente, proporcionando únicamente trazas de pared, historias puntuales o proyecciones integradas de línea de visión. Estas limitaciones dificultan la resolución de la dinámica volumétrica, particularmente para modos de propagación complejos como las detonaciones rotatorias en tubos redondos, donde los efectos de la capa límite en canales estrechos pueden alterar la física. Además, las aplicaciones tomográficas existentes a las detonaciones se han visto limitadas por ángulos de visión dispersos o por la falta de corrección por refracción a través de elementos ópticos curvos, lo que distorsiona los campos reconstruidos.

Metodología
Los autores presentan una campaña de tomografía de quimioluminiscencia con resolución temporal diseñada para reconstruir estructuras reactivas 3D en mezclas basadas en etileno a 1 atm. El montaje experimental utiliza un tubo redondo con acceso óptico con una sección de prueba de zafiro cilíndrica (DI 66.70 mm, DE 88.90 mm).

• Sistema de Imagen: Cinco cámaras de alta velocidad (frecuencias de cuadro de 2–5 MHz) están posicionadas alrededor de la sección de prueba en ángulos específicos (0°, 40°, 90°, 130°, 200°) para capturar la quimioluminiscencia de banda ancha (dominada por la emisión de CH*).
• Calibración: Se desarrolló un procedimiento de calibración personalizado para abordar la refracción causada por la ventana de zafiro cilíndrica. A diferencia de los métodos de objetivo plano estándar, este enfoque impone la ley de Snell en las interfaces del cilindro. Se tomó la imagen de un objetivo de rejilla de puntos de 9×9 a través del cilindro, y un algoritmo de Levenberg–Marquardt optimizó las posiciones de las cámaras, los parámetros del cilindro (radios, índice de refracción) y las posiciones de los puntos para minimizar los errores de reproyección (logrando un error RMS de ~0.15–0.3 px).
• Algoritmo de Reconstrucción: El campo de emisión 3D se reconstruye utilizando una rejilla cilíndrica ajustada a la pared con una base trilineal para resolver con precisión las estructuras adheridas a la pared del tubo. El problema inverso se formula como una optimización de mínimos cuadrados no negativos de gran escala. Minimiza un término de pérdida de datos (diferencia entre las proyecciones medidas y las predicciones del modelo) regularizado por una penalización de Tikhonov de segundo orden para imponer suavidad y suprimir artefactos. La optimización se resuelve utilizando un algoritmo de Broyden–Fletcher–Goldfarb–Shanno de memoria limitada (L-BFGS).

Contribuciones Clave

1. Tomografía Corregida por Refracción: El desarrollo de un marco de calibración personalizado que tiene en cuenta la refracción a través de una sección de prueba óptica cilíndrica, permitiendo un trazado de rayos preciso para la reconstrucción tomográfica en tubos redondos.
2. Imagen Volumétrica de Alta Velocidad: La aplicación de la tomografía de quimioluminiscencia de múltiples vistas a frecuencias de megahercios para resolver la evolución macroscópica de las ondas de detonación, yendo más allá de las limitaciones de las proyecciones 2D.
3. Extracción Cinemática Cuantitativa: Una metodología para extraer parámetros geométricos y cinemáticos clave (velocidades axiales/azimutales, ángulos de trayectoria, ángulos del tallo de Mach) directamente de los campos reconstruidos 3D.

Resultados
El método se aplicó a tres casos representativos:

• Caso A (Detonación Rotatoria Persistente): En una mezcla de C₂H₄–aire–Ar casi estequiométrica, la reconstrucción visualizó claramente una onda de rotación de un solo frente. El campo 3D reveló un fuerte frente adherido a la pared (tallo de Mach) y una estructura reactiva oblicua (onda transversal) que se encuentra en un punto triple. El choque incidente no fue visible, lo cual es consistente con el retraso en la liberación de calor detrás de él. El seguimiento de la trayectoria del punto triple arrojó una velocidad axial de 1388 m/s (coincidiendo con los datos de la sonda de iones) y una velocidad azimutal de 1468 m/s. El ángulo de trayectoria derivado ($\beta_1$) fue de 46.6°, y el paso fue de 2.97 diámetros de tubo, consistente con la literatura previa. El ángulo del tallo de Mach local ($\beta_2$) se estimó entre 80° y 90°, indicando un tallo fuerte y casi normal.
• Caso B (Detonación Rotatoria Fallida): En una mezcla rica en combustible, la onda se fragmentó. La tomografía aisló con éxito características en forma de anillo que emanaban de la pared y que estaban superpuestas sobre la emisión distribuida en las proyecciones 2D. La reconstrucción permitió el seguimiento de estos anillos, revelando una tasa de expansión azimutal de ~1720 m/s.
• Caso C (Frentes Transversales de Contrapropagación): En una mezcla pobre, dos frentes transversales colisionaron. La reconstrucción tomográfica separó las estructuras post-colisión, que aparecían de forma de convolución en las proyecciones 2D, revelando un bloque axial de movimiento rápido y uno de movimiento lento.

Significancia y Reclamaciones
El artículo sostiene que la tomografía de quimioluminiscencia con resolución temporal sirve como una herramienta no intrusiva capaz de resolver la dinámica volumétrica de la detonación que es difícil de inferir mediante diagnósticos puntuales, planos o de línea de visión. Al reconstruir el campo del término de fuente de emisión efectiva, los autores demuestran la capacidad de visualizar las morfologías de las ondas y determinar con precisión las velocidades de las ondas y los parámetros geométricos. Los autores señalan que, si bien la técnica proporciona un marcador para las estructuras reactivas, no mide directamente las variables de estado (densidad, temperatura) ni la estructura del choque en sí; por consiguiente, las regiones de emisión débil (como la estela del choque incidente) pueden permanecer indetectables. Este trabajo establece una base para comparaciones más estrictas para simulaciones y descripciones de orden reducido de la estabilidad y los límites de propagación de la detonación.