Anisotropic Cylindrical Waves in a Square Lattice of Acoustic Waveguides

Resumo Técnico: Ondas Cilíndricas Anisotrópicas em uma Rede Quadrada de Guias de Onda Acústicas

Definição do Problema
Embora a dinâmica de ondas não lineares em metamateriais acústicos aéreos unidimensionais (1D) — especificamente guias de onda carregados com ressonadores de Helmholtz — tenha sido extensivamente estudada, a exploração de tais fenômenos em configurações bidimensionais (2D) permanece limitada. Estudos 2D anteriores frequentemente utilizaram elementos ressonantes para suprimir a anisotropia inerente às redes quadradas, impondo assim uma simetria radial. Este trabalho investiga o cenário oposto: uma rede quadrada de guias de onda acústicas sem ressonadores de Helmholtz. O problema central é caracterizar como a anisotropia intrínseca da geometria da rede quadrada influencia a propagação de ondas cilíndricas, particularmente na interação entre a dispersão dependente da direção e a não linearidade.

Metodologia
O estudo emprega uma abordagem multifacetada combinando derivação teórica, validação experimental e simulação numérica:

1. Estrutura Teórica:

   • Análise de Dispersão: Os autores derivam a relação de dispersão para a rede quadrada utilizando o Método da Matriz de Transferência (TMM) sob uma aproximação monomodal. Isso estabelece a estrutura de bandas e identifica as direções de alta simetria ($\Gamma X$, $\Gamma M$ e $XM$).
   • Modelagem Efetiva: Para modelar a dinâmica de ondas de longo comprimento de onda e fraqueza não linear, os autores desenvolvem uma analogia eletroacústica (EA) aprimorada. Eles introduzem um esquema de discretização de "supercélula" refinado para capturar a dispersão anisotrópica com precisão. Isso leva a uma equação diferencial-diferença discreta (DDE) para o campo de pressão.
   • Limite Contínuo: Ao aplicar uma aproximação contínua à DDE, os autores derivam uma equação de Boussinesq 2D efetiva. Crucialmente, este modelo inclui derivadas espaciais de quarta ordem para dar conta da simetria quadrada, resultando em coeficientes dispersivos explicitamente dependentes do ângulo.
   • Redução para cKdV: No limite de grandes raios e pequenas amplitudes, a equação de Boussinesq 2D é reduzida a uma equação de Korteweg–de Vries cilíndrica (cKdV). Mostra-se que o coeficiente de dispersão nesta equação é uma função do ângulo de propagação $\theta$.

2. Validação Experimental:

   • Uma rede acústica quadrada física 2D ($15 \times 15$ células unitárias) foi construída utilizando obstáculos quadrados rígidos dispostos em uma rede com espaçamento $d = 12,5$ cm.
   • Ondas acústicas foram geradas por uma fonte de câmara de compressão central.
   • Campos de pressão foram medidos em vários nós utilizando um microfone portátil para capturar sinais temporais.
   • Relações de dispersão foram extraídas via Transformada Rápida de Fourier (FFT) 2D dos dados espaço-temporais.

3. Simulação Numérica:

   • Simulações numéricas diretas da dinâmica total da rede foram realizadas.
   • Adicionalmente, os autores utilizaram a equação de Westervelt 2D para simulações não lineares para corroborar as previsões teóricas.

Principais Resultados

• Dispersão Anisotrópica: O estudo confirma que a rede quadrada exibe uma forte dispersão anisotrópica. A propagação ao longo da direção diagonal ($\theta = \pi/4$, ou $\Gamma M$) é indispersiva, enquanto a propagação ao longo dos eixos principais ($\theta = 0$, ou $\Gamma X$) é altamente dispersiva. Este comportamento é atribuído ao espalhamento múltiplo em fase ao longo das diagonais devido à simetria quadrada.
• Precisão do Modelo: A equação de Boussinesq 2D aprimorada derivada, e subsequentemente a equação cKdV dependente do ângulo, apresentam excelente concordância tanto com as previsões analíticas do TMM quanto com as medições experimentais através da primeira banda de propagação. O modelo captura com sucesso a transição do comportamento dispersivo para o indispersivo conforme o ângulo de propagação muda.
• Evolução da Forma de Onda:
  • Regime Linear: A rede suporta ondas cilíndricas com formas de onda que variam significativamente com a direção. Ao longo de direções dispersivas ($\theta \neq \pi/4$), as ondas evoluem para pacotes do tipo Airy. Ao longo da diagonal indispersiva ($\theta = \pi/4$), as ondas mantêm uma forma de pulso quase indispersiva.
  • Regime Não Linear: A interação entre anisotropia e não linearidade produz regimes de propagação distintos. Para $\theta \neq \pi/4$, o sistema suporta ondas solitárias suaves (solitons cilíndricos). No entanto, conforme $\theta \to 0$ (aproximando-se do limite indispersivo ao longo da diagonal não é o caso aqui; em vez disso, a redução para a equação de Burgers radial não viscosa ocorre quando a dispersão desaparece em limites específicos, embora o texto note que para $\theta \to 0$ no contexto da cKdV, a equação se reduz para suportar estruturas do tipo choque), as formas de onda transitam para frentes agudas, do tipo choque. Inversamente, ao longo da diagonal indispersiva, a dinâmica é governada pela equação de Burgers generalizada, suportando ondas de choque cilíndricas.
• Corroboração Numérica: Simulações da rede completa e da equação de Westervelt 2D confirmam a coexistência desses distintos regimes não lineares dependendo da direção de propagação.

Significância e Alegações
O artigo alega estabelecer um arcabouço geral para prever fenômenos de ondas dependentes da direção em metamateriais acústicos 2D. Sua principal significância reside em demonstrar que a anisotropia intrínseca de uma simples rede quadrada, mesmo sem elementos ressonantes, influencia profundamente a evolução de ondas radialmente simétricas. O trabalho revela uma rica família de ondas cilíndricas anisotrópicas onde características como a forma de onda (de pacotes tipo Airy a estruturas tipo choque) e propriedades de dispersão variam continuamente com o ângulo de propagação. Ao validar essas descobertas experimental e teoricamente, os autores fornecem um modelo robusto para compreender e manipular a propagação de ondas em redes acústicas estruturadas 2D, estendendo a compreensão da dinâmica de ondas não lineares além dos ambientes 1D.