In this Note, we introduce the Fourier Singular Complement Method, for solving Maxwell equations in a 3D prismatic domain. The numerical implementation of this method provides a continuous approximation of the electromagnetic field. It can be applied to the computation of propagating and evanescent modes in prismatic stub filters, thus generalizing 2D approaches.
RésuméRésolution des équations de Maxwell dans des domaines prismatiques tridimensionnels. Dans cette Note, nous introduisons la Méthode du Complément Singulier avec Fourier, pour résoudre les équations de Maxwell dans des domaines prismatiques tridimensionnels. La mise en oeuvre numérique de cette méthode permet de calculer une approximation continue du champ électromagnétique. Elle peut être appliquée à la détermination des modes propagatifs ou bloquants dans un filtre à stubs prismatique, ce qui constitue une généralisation des méthodes applicables en domaine bidimensionnel.
Version française abrégéeLes nombres entre parenthèses renvoient à la version anglaise. L'opérateur rotationnel est noté curl. Nous nous intéressons à la résolution des équations de Maxwell dans un domaine Ω prismatique, i.e. Ω = ω × ]0, L[ (voir la Fig. 1), à l'aide d'une généralisation de la méthode du Complément Singulier, bien connue en domaine 2D [2,1,10,8,11]. Elle fournit une discrétisation continue du champ électromagnétique, avec prise en compte explicite des conditions aux limites. On souhaite simuler la propagation d'ondes électromagnétiques autour E-mail addresses: Patrick.Ciarlet@ensta.fr (P. Ciarlet, Jr.), garciaemm@yahoo.fr (E. Garcia), zou@math.cuhk.edu.hk (J. Zou). Ciarlet, Jr., et al. / C. R. Acad. Sci. Paris, Ser. I 339 (2004) [721][722][723][724][725][726] d'un conducteur parfait. On se ramène à un domaine borné, en introduisant une frontière artificielle, notée Γ A , sur laquelle on impose une condition aux limites de type Silver-Müller. Il est important de noter que le champ est régulier au voisinage de cette frontière artificielle [6]. La frontière du domaine de calcul Ω est donc décomposée en ∂Ω = Γ C ∪ Γ A ; on appelle n la normale unitaire extérieure à ∂Ω. Soit (E, H) le champ électromagnétique, fonction de x et t. Les équations de Maxwell dans le vide sont relatées en (1) et (2). Les densités de charge ρ et de courant J vérifient l'équation de conservation de la charge ∂ t ρ + div J = 0, et la connaissance de e permet de modéliser les ondes planes incidentes (e = 0), ou la frontière absorbante (e = 0). A la suite de [3,4], on peut reformuler de façon équivalente (1), (2), en remplaçant les deux premières équations de (1) par (3), (4), ce qui permet d'obtenir deux systèmes « découplés » du second ordre en temps, en E et en H. Dans la suite, on considère uniquement le système en E, sachant que celui en H peut être résolu similairement. On peut prouver (cf.[4]) que, sous des hypothèses de régularités des données ad hoc, le champ E est l'unique solution de la Formulation Variationnelle (6), de régularité en temps (7). Si on suppose que la densi...