Anisotropic Cylindrical Waves in a Square Lattice of Acoustic Waveguides

Résumé technique : Ondes cylindriques anisotropes dans un réseau carré de guides d'ondes acoustiques

Énoncé du problème
Bien que la dynamique des ondes non linéaires dans les métamatériaux acoustiques aériens unidimensionnels (1D) (spécifiquement les guides d'ondes chargés de résonateurs de Helmholtz) ait été largement étudiée, l'exploration de tels phénomènes dans des contextes bidimensionnels (2D) reste limitée. Les études antérieures en 2D ont souvent utilisé des éléments résonants pour supprimer l'anisotropie inhérente aux réseaux carrés, afin d'imposer une symétrie radiale. Ce travail étudie le scénario opposé : un réseau de guides d'ondes acoustiques carrés sans résonateurs de Helmholtz. Le problème central est de caractériser comment l'anisotropie intrinsèque de la géométrie du réseau carré influence la propagation des ondes cylindriques, en particulier dans l'interaction entre la dispersion dépendante de la direction et la non-linéarité.

Méthodologie
L'étude emploie une approche multidimensionnelle combinant dérivation théorique, validation expérimentale et simulation numérique :

1. Cadre théorique :

   • Analyse de la dispersion : Les auteurs dérivent la relation de dispersion pour le réseau carré en utilisant la méthode de la matrice de transfert (TMM) sous une approximation monomodale. Cela établit la structure de bandes et identifie les directions de haute symétrie ($\Gamma X$, $\Gamma M$, et $XM$).
   • Modélisation effective : Pour modéliser la dynamique à grande longueur d'onde et faible non-linéarité, les auteurs développent une analogie électroacoustique (EA) améliorée. Ils introduisent un schéma de discrétisation par « supercellule » raffiné pour capturer précisément la dispersion anisotrope. Cela conduit à une équation de différence-différentielle (DDE) discrète pour le champ de pression.
   • Limite continue : En appliquant une approximation continue à la DDE, les auteurs dérivent une équation de Boussinesq 2D effective. Crucialement, ce modèle inclut des dérivées spatiales d'ordre quatre pour rendre compte de la symétrie carrée, ce qui se traduit par des coefficients de dispersion explicites dépendant de l'angle.
   • Réduction vers la cKdV : Dans la limite de grands rayons et de faibles amplitudes, l'équation de Boussinesq 2D est réduite à une équation de Korteweg–de Vries cylindrique (cKdV). Il est démontré que le coefficient de dispersion dans cette équation est une fonction de l'angle de propagation $\theta$.

2. Validation expérimentale :

   • Un réseau acoustique carré physique 2D ($15 \times 15$ cellules unitaires) a été construit en utilisant des obstacles carrés rigides disposés sur un réseau avec un espacement $d = 12,5$ cm.
   • Les ondes acoustiques ont été générées à partir d'une source de chambre de compression centrale.
   • Les champs de pression ont été mesurés en divers nœuds à l'aide d'un microphone portable pour capturer les signaux temporels.
   • Les relations de dispersion ont été extraites via une transformée de Fourier rapide (FFT) 2D des données spatio-temporelles.

3. Simulation numérique :

   • Des simulations numériques directes de la dynamique complète du réseau ont été effectuées.
   • De plus, les auteurs ont utilisé l'équation de Westervelt 2D pour les simulations non linéaires afin de corroborer les prédictions théoriques.

Résultats clés

• Dispersion anisotrope : L'étude confirme que le réseau carré présente une forte dispersion anisotrope. La propagation selon la direction diagonale ($\theta = \pi/4$, ou $\Gamma M$) est sans dispersion, tandis que la propagation selon les axes principaux ($\theta = 0$, ou $\Gamma X$) est hautement dispersive. Ce comportement est attribué à la diffusion multiple en phase le long des diagonales due à la symétrie carrée.
• Précision du modèle : L'équation de Boussinesque 2D améliorée, et par la suite l'équation cKdV dépendante de l'angle, montrent un excellent accord avec les prédictions analytiques de la TMM ainsi qu'avec les mesures expérimentales à travers la première bande de propagation. Le modèle capture avec succès la transition du comportement dispersif vers le comportement sans dispersion lorsque l'angle de propagation change.
• Évolution de la forme d'onde :
  • Régime linéaire : Le réseau supporte des ondes cylindriques dont les formes d'onde varient considérablement selon la direction. Le long des directions dispersives ($\theta \neq \pi/4$), les ondes évoluent en paquets de type Airy. Le long de la diagonale sans dispersion ($\theta = \pi/4$), les on waves maintiennent une forme d'impulsion quasi sans dispersion.
  • Régime non linéaire : L'interaction entre l'anisotropie et la non-linéarité produit des régimes de propagation distincts. Pour $\theta \neq \pi/4$, le système supporte des ondes solitaires lisses (solitons cylindriques). Cependant, à mesure que $\theta \to 0$ (l'approche de la limite sans dispersion le long de la diagonale n'est pas le cas ici ; plutôt, la réduction à l'équation de Burgers radiale non visqueuse se produit lorsque la dispersion disparaît dans des limites spécifiques, bien que le texte note que pour $\theta \to 0$ dans le contexte de la cKdV, l'équation se réduit pour supporter des structures de type choc), les formes d'onde transitent vers des fronts de choc abrupts. Inversement, le long de la diagonale sans dispersion, la dynamique est gouvernée par l'équation de Burgers généralisée, supportant des ondes de choc cylindriques.
• Corrélation numérique : Les simulations du réseau complet et de l'équation de Westervelt 2D confirment la coexistence de ces différents régimes non linéaires selon la direction de propagation.

Signification et affirmations
L'article affirme établir un cadre général pour prédire les phénomènes d'ondes dépendant de la direction dans les métamatériaux acoustiques 2D. Sa principale importance réside dans la démonstration que l'anisotropie intrinsèque d'un simple réseau carré, même sans éléments résonants, influence profondément l'évolution des ondes radialement symétriques. Ce travail révèle une riche famille d'ondes cylindriques anisotropes où les caractéristiques telles que la forme d'onde (de type Airy à type choc) et les propriétés de dispersion varient continuellement avec l'angle de propagation. En validant ces découvertes de manière expérimentale et théorique, les auteurs fournissent un modèle robuste pour comprendre et manipuler la propagation des ondes dans les réseaux acoustiques structurés en 2D, étendant ainsi la compréhension de la dynamique des ondes non linéaires au-delà des contextes 1D.