Nouvelles formulations intégrales pour les problèmes de diffraction d'ondes

We present an integral equation method for solving boundary value problems of the Helmholtz equation in unbounded domains. The method relies on the factorisation of one of the Calderón projectors by an operator approximating the exterior admittance (Dirichlet to Neumann) operator of the scattering o...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:ESAIM. Mathematical modelling and numerical analysis 2004-01, Vol.38 (1), p.157-175
Hauptverfasser: Levadoux, David P., Michielsen, Bastiaan L.
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:We present an integral equation method for solving boundary value problems of the Helmholtz equation in unbounded domains. The method relies on the factorisation of one of the Calderón projectors by an operator approximating the exterior admittance (Dirichlet to Neumann) operator of the scattering obstacle. We show how the pseudo-differential calculus allows us to construct such approximations and that this yields integral equations without internal resonances and being well-conditioned at all frequencies. An implementation technique is elaborated, where again reasonings from pseudo-differential calculus play an important rôle. Some numerical examples are presented which appear to confirm that the new integral equation leads to linear systems which are much better conditioned than the classical ("direct") integral equations and hence have much better behaviour when solved with iterative techniques and matrix sparsification. On présente un formalisme intégral destiné à la résolution de l'équation de Helmholtz en domaine non borné. Ce formalisme repose sur la factorisation de l'un des projecteurs de Calderón par un opérateur dont la vocation est d'approcher au mieux l'admittance externe de l'objet diffractant. On montre alors comment le calcul pseudo-différentiel permet d'envisager la construction d'approximations conduisant à des équations intégrales sans résonance, bien posées à toutes les fréquences. Une technique de mise en œuvre est ensuite exposée, où de nombreux arguments provenant du calcul pseudo-différentiel interviennent encore. Enfin, on présente quelques résultats numériques venant conforter la démarche. Principalement, on constate que les systèmes issus du nouveau formalisme sont très bien conditionnés comparés à ceux provenant d'équations plus classiques.
ISSN:0764-583X
1290-3841
DOI:10.1051/m2an:2004008