Anisotropic Cylindrical Waves in a Square Lattice of Acoustic Waveguides

Resumen Técnico: Ondas Cilíndricas Anisotrópicas en una Red Cuadrada de Guías de Onda Acústicas

Planteamiento del Problema
Si bien la dinámica de ondas no lineales en metamateriales acústicos aéreos unidimensionales (1D) (específicamente guías de onda cargadas con resonadores de Helmholtz) ha sido ampliamente estudiada, la exploración de tales fenómenos en entornos bidimensionales (2D) sigue siendo limitada. Estudios previos en 2D utilizaron frecuentemente elementos resonantes para suprimir la anisotropía inherente de las redes cuadradas, imponiendo así una simetría radial. Este trabajo investiga el escenario opuesto: una red cuadrada de guías de onda acústicas sin resonadores de Helmholtz. El problema central es caracterizar cómo la anisotropía intrínseca de la geometría de la red cuadrada influye en la propagación de ondas cilíndricas, particularmente en la interacción entre la dispersión dependiente de la dirección y la no linealidad.

Metodología
El estudio emplea un enfoque multifacético que combina derivación teórica, validación experimental y simulación numérica:

1. Marco Teórico:

   • Análisis de Dispersión: Los autores derivan la relación de dispersión para la red cuadrada utilizando el Método de la Matriz de Transferencia (TMM) bajo una aproximación monomodal. Esto establece la estructura de bandas e identifica las direcciones de alta simetría ($\Gamma X$, $\Gamma M$, y $XM$).
   • Modelado Efectivo: Para modelar la dinámica de longitudes de onda largas y débilmente no lineales, los autores desarrollan una analogía electroacústica (EA) mejorada. Introducen un esquema de discretización de "supercelda" refinado para capturar con precisión la dispersión anisotrópica. Esto conduce a una ecuación diferencial-de diferencia discreta (DDE) para el campo de presión.
   • Límite Continuo: Al aplicar una aproximación de continuo a la DDE, los autores derivan una ecuación de Boussinesq 2D efectiva. Crucialmente, este modelo incluye derivadas espaciales de cuarto orden para dar cuenta de la simetría cuadrada, lo que resulta en coeficientes dispersivos explícitos que dependen del ángulo.
   • Reducción a cKdV: En el límite de radios grandes y amplitudes pequeñas, la ecuación de Boussinesq 2D se reduce a una ecuación de Korteweg–de Vries cilíndrica (cKdV). Se demuestra que el coeficiente de dispersión en esta ecuación es una función del ángulo de propagación $\theta$.

2. Validación Experimental:

   • Se construyó una red acústica cuadrada física de 2D ($15 \times 15$ celdas unitarias) utilizando obstáculos cuadrados rígidos dispuestos en una red con un espaciamiento de $d = 12.5$ cm.
   • Las ondas acústicas fueron generadas mediante una fuente de cámara de compresión central.
   • Los campos de presión se midieron en varios nodos utilizando un micrófono portátil para capturar las señales temporales.
   • Las relaciones de dispersión se extrajeron mediante la Transformada Rápida de Fourier (FFT) 2D de los datos espaciotemporales.

3. Simulación Numérica:

   • Se realizaron simulaciones numéricas directas de la dinámica completa de la red.
   • Adicionalmente, los autores utilizaron la ecuación de Westervelt 2D para simulaciones no lineales para corroborar las predicciones teóricas.

Resultados Clave

• Dispersión Anisotrópica: El estudio confirma que la red cuadrada exhibe una fuerte dispersión anisotrópica. La propagación a lo largo de la dirección diagonal ($\theta = \pi/4$, o $\Gamma M$) es libre de dispersión, mientras que la propagación a lo largo de los ejes principales ($\theta = 0$, o $\Gamma X$) es altamente dispersiva. Este comportamiento se atribuye a la dispersión en fase de la propagación múltiple a lo largo de las diagonales debido a la simetría cuadrada.
• Precisión del Modelo: La ecuación de Boussinesq 2D mejorada derivada, y posteriormente la ecuación cKdV dependiente del ángulo, muestran una excelente concordancia tanto con las predicciones analíticas de TMM como con las mediciones experimentales a través de la primera banda de propagación. El modelo captura con éxito la transición del comportamiento dispersivo al comportamiento libre de dispersión a medida que cambia el ángulo de propagación.
• Evolución de la Forma de Onda:
  • Régimen Lineal: La red soporta ondas cilíndricas con formas de onda que varían significativamente según la dirección. A lo largo de las direcciones dispersivas ($\theta \neq \pi/4$), las ondas evolucionan hacia paquetes tipo Airy. A lo largo de la diagonal libre de dispersión ($\theta = \pi/4$), las ondas mantienen una forma de pulso casi sin dispersión.
  • Régimen No Lineal: La interacción entre la anisotropía y la no linealidad produce distintos regímenes de propagación. Para $\theta \neq \pi/4$, el sistema soporta ondas solitarias suaves (solitones cilíndricos). Sin embargo, a medida que $\theta \to 0$ (el texto señala que para $\theta \to 0$ en el contexto de cKdV, la ecuación se reduce para soportar estructuras tipo choque; es decir, cuando la dispersión desaparece en límites específicos, ocurre la reducción a la ecuación de Burgers radial no viscosa), las formas de onda transicionan hacia frentes agudos tipo choque. Por el contrario, a lo largo de la diagonal libre de dispersión, la dinámica está gobernada por la ecuación de Burgers generalizada, soportando ondas de choque cilíndricas.
• Corroboración Numérica: Las simulaciones de la red completa y de la ecuación de Westervelt 2D confirman la coexistencia de estos distintos regímenes no lineales dependiendo de la dirección de propagación.

Significancia y Reivindicaciones
El artículo afirma establecer un marco general para predecir fenómenos de ondas dependientes de la dirección en metamateriales acústicos 2D. Su principal significancia radica en demostrar que la anisotropía intrínseca de una simple red cuadrada, incluso sin elementos resonantes, influye profundamente en la evolución de las ondas radialmente simétricas. El trabajo revela una rica familia de ondas cilíndricas anisotrópicas donde las características, tales como la forma de la onda (desde tipos Airy hasta tipos choque) y las propiedades de dispersión, varían continuamente con el ángulo de propagación. Al validar estos hallazgos experimental y teóricamente, los autores proporcionan un modelo robusto para comprender y manipular la propagación de ondas en redes acústicas estructuradas en 2D, extendiendo la comprensión de la dinámica de ondas no lineales más allá de los entornos 1D.