Approximations unidirectionnelles de la propagation acoustique en guides d’ondes irréguliers

Dans l’environnement urbain, l’essentiel des nuisances sonores est produit par les moyens de transport. La forte concentration des axes de communication à l’intérieur des villes a pour effet d’exposer en permanence les citadins à de forts niveaux sonores. Perçu comme une nuisance dans le confort de vie, le bruit des transports présente également un risque sanitaire. Un récent rapport de l’Organisation Mondiale de la Santé [1] établit un lien direct entre l’exposition continue aux bruits environnementaux et le risque de maladies cardio-vasculaires, les troubles du sommeil et les troubles cognitifs chez les enfants. Depuis plusieurs années, la législation  impose une réduction des nuisances sonores en milieu urbain. À cet effet, une directive européenne  exige une évaluation de l’exposition aux bruits des transports pour toutes les villes de plus de 25000 habitants. Cette évaluation a pour support des cartographies de bruit, qui sont utilisées pour mettre en œuvre les plans d’actions dans le traitement de ces nuisances. La réalisation de ces cartes de bruit est donc un enjeu majeur de la lutte contre les nuisances sonores.

La taille étendue des axes de communication rendant irréalisables des campagnes de mesures à l’échelle d’une ville, les cartographies de bruit sont générées par des logiciels de prédictions. Les renseignements importés dans ces logiciels sont les données géométriques et topographiques de la ville, ainsi que la description du trafic routier (nombre et type de véhicules) en fonction du type de voiries (voiries locales, artères principales, réseaux secondaire et autoroutier). Ces données sont traitées par un modèle de bruit routier avant d’être intégrées dans un logiciel de prédictions de la propagation acoustique. Un recalage de ce modèle avec des mesures ponctuelles in situ permet de générer une cartographie des nuisances sonores. Les logiciels de prédictions communément utilisés sont basés sur une approche énergétique de la propagation des ondes acoustiques. Cette approche ne prend pas en compte la phase des ondes acoustiques, ce qui limite sa validité pour les basses fréquences.

Généralités sur les équations paraboliques

Très souvent associée à l’approximation unidirectionnelle, l’équation parabolique est une forme d’équation d’onde unidirectionnelle couramment utilisée dans de nombreux domaines de l’acoustique, tels que l’acoustique sous-marine [22,44,91], la propagation en milieu extérieur [14, 43] ou encore la propagation d’ondes sismiques [94]. Ce type d’équation d’onde est largement utilisé pour sa facilité de résolution et son aptitude à prendre en compte les paramètres du milieu de propagation (hétérogénéité du milieu de propagation [9], milieu en mouvement [17,72], géométries des frontières [27,59]). Aussi connu sous le nom de Beam Propagation Method dans le domaine de l’optique [34,48,96], ce type d’équation est apparu en premier lieu dans le domaine de l’électromagnétisme pour résoudre des problèmes de propagation d’ondes le long de la surface de la Terre [57].

Formulation multimodale de la propagation acoustique unidirectionnelle 

La méthode multimodale permet de modéliser la propagation acoustique dans des guides d’ondes de géométrie complexe [15, 40, 41, 74, 75]. Basée sur le développement du champ de pression sur les modes propres d’un guide d’ondes, cette méthode permet de prendre en compte tous les modes, propagatifs et évanescents, qui contribuent de façon significative au champ acoustique afin de caractériser avec précision la propagation. Nous proposons d’utiliser cette méthode dans le cadre de l’approximation unidirectionnelle afin d’étudier les couplages entre les modes induits par des variations de section transverse du guide d’ondes. Cela permettra par la suite de mieux appréhender la limite de validité de cette approche en fonction de la fréquence, de la géométrie du guide d’ondes ou encore de la nature de la source. Après avoir présenté la méthode multimodale pour la propagation acoustique en guide de section variable, une décomposition du champ de pression en fonction du sens de propagation est réalisée. Cela permettra de formuler une approximation unidirectionnelle de cette méthode.

La méthode multimodale ainsi que son approximation unidirectionnelle ont été présentées. Issues de la décomposition du champ de pression selon son sens de propagation, deux formes d’équations vectorielles unidirectionnelles ont été proposées. Par rapport à la méthode multimodale sous sa forme traditionnelle, il n’est pas nécessaire de calculer de matrice d’admittance. Contrairement aux méthodes de modes couplés utilisées en acoustique sous marine (par exemple : [3,32,66]), qui sont basées sur la projection de l’équation parabolique, les équations (1.36) et (1.37) sont issues de la projection de l’équation d’onde sans approximation sur les modes transverses. Ces équations d’ondes unidirectionnelles ont donc l’avantage de ne pas comporter de limitation angulaire dans leurs formulations. A notre connaisance, cette forme d’approximation unidirectionnelle du champ de pression vectoriel est originale.

Propagation unidirectionnelle dans une rue droite

Le fait de réaliser une approximation unidirectionnelle de la propagation acoustique en milieu urbain est motivé par la limitation du coût numérique de la résolution de l’équation de Helmholtz à basses fréquences (pour des longueurs d’ondes de même ordre de grandeur que les dimensions de la rue). C’est pourquoi la méthode de résolution de l’équation parabolique doit être construite de manière à être la moins coûteuse possible en terme de temps de calcul. Cette première partie présente, pour le cas particulier des rues droites, la méthode de résolution de l’équation parabolique tridimensionnelle en guide d’ondes ouvert.

Validation de la méthode par comparaison avec des résultats expérimentaux

Présentation du dispositif expérimental 

La méthode de résolution de l’équation parabolique tridimensionnelle, précédemment présentée, est validée par comparaison avec des résultats expérimentaux. Ces résultats proviennent de mesures réalisées sur une maquette de rue à l’échelle 1/35ème, placée dans une salle semianéchoïque, caractérisée par un sol lisse ainsi que des murs et un plafond traités avec de la mousse de mélamine (efficace jusqu’à 1 kHz à pleine échelle). Cette maquette est composée par deux rangées de cubes en plexiglas afin de reproduire les bâtiments d’une rue. Les deux rangées de cubes mesurent 30 cm de haut et sont espacées de 20 cm, ce qui correspond à pleine échelle à une rue de 10.5 m de hauteur et de 7 m de largeur. Une grande plaque de plexiglas sert de support à cette maquette afin de représenter le sol  . La matière constituant les éléments de la maquette permet d’approcher des conditions de paroi rigide au sol et sur les façades des bâtiments. Enfin, une terminaison anéchoïque est placée en sortie de la maquette afin d’éviter les réflexions en sortie de la rue .

La mesure du champ de pression à l’intérieur de la rue est réalisée par une antenne de huit microphones quart de pouce B&K 4190 . Ces huit microphones sont connectés à des pré-amplificateurs (B&K 2669) avant d’être reliés à un amplificateur de conditionnement (B&K Nexus 2693). Cette antenne microphonique est fixée sur le bras d’un robot comportant trois axes de déplacement, ce qui permet de cartographier le champ de pression dans plusieurs plans de la maquette. Le positionnement des cubes est renseigné dans un logiciel pilotant le robot et gérant l’acquisition des données. Les signaux acoustiques captés par les microphones sont transmis à une carte d’acquisition multi-voies.

La source de la maquette est constituée d’un haut-parleur raccordé à l’entrée de la rue par l’intermédiaire d’un tube carré (de 5 cm de côté). Un microphone de référence est placé à l’entrée du tube afin de capter l’amplitude du signal source. Un bafflage est réalisé à l’extrémité du tube par une plaque en bois , ce qui permet d’orienter le rayonnement acoustique de la source dans la direction de la rue. Le haut-parleur est relié à un amplificateur de puissance alimenté par un générateur basses fréquences produisant un signal sinusoïdal à une fréquence donnée.

Le signal source étant monochromatique, l’interprétation des signaux mesurés est assurée par une technique de détection synchrone réalisée à partir du microphone de référence. Cette technique permet de détecter avec précision un signal de faible amplitude en milieu bruité. A chaque point de mesure, le champ de pression est alors caractérisé par une amplitude et une phase. Avant chaque cycle de mesure, les huit microphones sont calibrés par rapport au microphone de référence .

Table des matières

Introduction
1 Approximation unidirectionnelle de la propagation acoustique
1.1 Généralités sur les équations paraboliques
1.1.1 Formulation de l’équation parabolique petit angle
1.1.2 Décomposition de l’équation de Helmholtz
1.1.3 Splitting de l’équation d’onde conduisant aux équations paraboliques petit angle
1.1.4 Equation parabolique grand angle
1.1.5 Conclusion
1.2 Formulation multimodale de la propagation acoustique unidirectionnelle
1.2.1 Généralités sur la méthode multimodale en guide de section variable
1.2.2 Décomposition du champ de pression
1.2.3 Approximation unidirectionnelle
1.2.4 Conclusion
Annexe 1.A Projection de l’équation d’onde sur les modes transverses
I Utilisation de l’équation parabolique pour une application à l’acoustique urbaine
2 Higher order Padé approximation for the parabolic equation in a varying cross-section waveguide
2.1 Introduction
2.2 Modified parabolic equation method
2.2.1 Formulation
2.2.2 Boundary condition
2.2.3 Numerical scheme
2.3 Validation of the modified parabolic equation
2.4 Range of validity
2.5 Conclusion
Appendix 2.A Finite Difference implementation of the boundary condition
Appendix 2.B Rotation angles in the complex plane for complex Padé coefficients
3 Résolution de l’équation parabolique tridimensionnelle dans des guides d’ondes ouverts
3.1 Propagation unidirectionnelle dans une rue droite
3.1.1 Formulation d’une équation parabolique 3D en guide d’ondes ouvert
3.1.2 Résolution numérique de l’équation parabolique 3D
3.1.3 Validation de la méthode par comparaison avec des résultats expérimentaux
3.1.4 Conclusion
3.2 Variation brusque de la section du guide d’ondes
3.2.1 Approximation de Kirchhoff
3.2.2 Validation par comparaison avec des résultats expérimentaux
3.2.3 Conclusion
3.3 Variation continue de la section du guide d’ondes
3.3.1 Formulation de l’équation parabolique 3D modifiée
3.3.2 Modification de la méthode de résolution de l’équation parabolique 3D
3.3.3 Application de la méthode
3.3.4 Conclusion
Annexe 3.A Intégration des PML dans la discrétisation par différences finies
II Approche multimodale de l’approximation unidirectionnelle de la propagation acoustique
4 Coarse-grid computation of the one-way propagation of coupled modes
4.1 Introduction
4.2 Multimodal admittance method
4.2.1 Coupled mode equations
4.2.2 Numerical integration
4.3 One-way approximation
4.4 Accuracy of the one-way method for coarse discretization
4.5 One-dimensional analytical solutions
4.6 Conclusion
5 Approche multimodale de la décomposition du champ de pression en série de
Bremmer
5.1 Multidiffusion dans un guide d’ondes discontinu
5.1.1 Méthode multimodale pour les guides avec discontinuités de section
5.1.2 Utilisation des matrices de diffusion pour la décomposition du champ de pression
5.1.3 Multidiffusion dans un créneau
5.1.4 Convergence de la série
5.1.5 Sens de propagation des modes couplés
5.1.6 Conclusion
5.2 Série de Bremmer dans des guides d’ondes bidimensionnels à variation de section continue
5.2.1 Décomposition en ordres : introduction
5.2.2 Décomposition en ordres : résolution itérative
5.2.3 Méthode de résolution numérique
5.2.4 Approximation unidirectionnelle
5.2.5 Reconstruction du champ de pression jusqu’à un ordre donné
5.2.6 Modélisation du champ de pression transmis à travers un guide d’ondes
irrégulier
5.2.7 Analogie avec l’approximation de Born appliquée à une équation d’onde
unidimensionnelle non-homogène
5.2.8 Conclusion
Conclusions et perspectives
Bibliographie

Télécharger le rapport complet