Mes travaux de recherche

Thèse de doctorat en mathématiques appliquées

Construction de conditions transparentes pour les guides d’ondes électromagnétiques, analyse et applications

Comment simuler sur un ordinateur un guide qui, dans le modèle physique, se prolonge jusqu’à l’infini ? Mon travail consiste à refermer artificiellement le domaine de calcul sans que cette frontière crée de faux échos.

ENSTA Paris · Équipe POEMS Soutenue le 17 décembre 2025 Analyse · Modélisation · Simulation

Cette thèse porte sur la propagation et la diffraction des ondes électromagnétiques dans des structures formées de guides reliés à une cavité ou comportant une perturbation : coude, obstacle, raccord ou inclusion de matériau.

Le problème en quelques mots

Un guide d’ondes canalise une onde dans une direction privilégiée. On en rencontre notamment dans les systèmes micro-ondes, les radars, la caractérisation de matériaux et certains composants photoniques. Dans le modèle mathématique, ces guides sont souvent semi-infinis. Un calcul numérique ne peut pourtant porter que sur une région bornée.

Il faut donc couper le guide. Une condition de bord ordinaire se comporterait comme un mur et renverrait artificiellement l’onde vers l’intérieur. Une condition transparente remplace la partie supprimée du guide : elle laisse sortir les ondes comme si le domaine se poursuivait réellement à l’infini.

1. 1Le guide réelUn domaine non borné, impossible à mailler entièrement.
2. 2La coupure numériqueOn conserve la zone où se trouvent les obstacles et les variations de matériaux.
3. 3La condition transparenteElle restitue exactement ou approximativement le comportement de la partie supprimée.

Le modèle mathématique : les équations de Maxwell

En régime harmonique, les champs électrique \(\mathbf E\) et magnétique \(\mathbf H\) satisfont les équations de Maxwell :

\[\begin{split}\left\{ \begin{aligned} \operatorname{curl}\mathbf E-i\omega\mu\mathbf H &= 0,\\ \operatorname{curl}\mathbf H+i\omega\varepsilon\mathbf E &= \mathbf J,\\ \operatorname{div}(\varepsilon\mathbf E) &= 0,\\ \operatorname{div}(\mu\mathbf H) &= 0. \end{aligned} \right.\end{split}\]

Ici, \(\varepsilon\) et \(\mu\) décrivent les propriétés du matériau, \(\omega\) est la fréquence et \(\mathbf J\) une source. Sur les parois métalliques, on impose une condition de conducteur parfait. La difficulté est de compléter ce modèle sur les sections artificielles où le guide a été coupé.

Deux domaines de calcul tronqués et maillés : la condition transparente est placée aux extrémités artificielles.

Pourquoi les guides hétérogènes sont plus difficiles

Dans un guide homogène, les propriétés du matériau sont constantes dans la section. Les champs peuvent alors se décomposer en modes classiques de type TE, TM ou TEM, qui jouent le rôle de briques élémentaires de la solution.

Dans un guide hétérogène, la permittivité ou la perméabilité varie dans la section. Le problème modal n’est plus auto-adjoint : les modes ne sont plus nécessairement orthogonaux et leur complétude n’est pas disponible en général. La décomposition des champs et la sélection des ondes sortantes deviennent donc nettement plus délicates.

Un phénomène particulièrement important est celui des modes inverses : leur vitesse de phase et leur vitesse de groupe ont des signes opposés. Autrement dit, les oscillations semblent avancer dans un sens alors que l’énergie se transporte dans l’autre. Cette situation peut mettre en défaut certaines couches absorbantes classiques.

Comment lire l’animation ? La courbe bleue représente l’onde et les pointillés verts son enveloppe. Le point orange suit une crête qui se déplace vers la gauche : il matérialise la vitesse de phase. Le point vert suit l’enveloppe vers la droite : il indique la vitesse de groupe, donc le sens de transport de l’énergie. Les deux vitesses ont ici des signes opposés : c’est précisément un mode inverse.

Les principales contributions

Comprendre les modes

Établissement des propriétés modales utiles dans les guides hétérogènes : caractère discret, relations de biorthogonalité, lien entre flux d’énergie et vitesse de groupe, et analyse de l’apparition des modes inverses.

Construire des frontières transparentes

Formulation rigoureuse de conditions qui raccordent la solution calculée par éléments finis à une représentation modale dans le guide. Un recouvrement entre ces deux descriptions contrôle l’erreur lorsque la série de modes est tronquée.

Valider numériquement

Calcul des modes et implémentation des opérateurs avec XLiFE++, à l’aide d’éléments finis de Nédélec. Les simulations couvrent des guides homogènes, hétérogènes, multiplement connexes et perturbés.

Deux conditions transparentes complémentaires

Les deux méthodes utilisent les modes du guide et une zone de recouvrement : pendant une courte distance, la description par éléments finis et la description modale coexistent. Ce recouvrement rend la condition compacte et amortit naturellement l’effet des modes évanescents qui ne sont pas conservés dans le calcul.

EtM

Electric-to-Magnetic

L’opérateur relie la trace tangentielle du champ électrique sur une section intérieure à la trace correspondante du champ magnétique sur la frontière artificielle.

• proche des opérateurs Dirichlet-to-Neumann classiques ;
• compact grâce au recouvrement ;
• soumis à certaines fréquences parasites, dites modes de boîte ;
• moins adapté aux matériaux anisotropes généraux.

CtM

Currents-to-Magnetic

Ce nouvel opérateur combine les traces tangentielles électrique et magnétique sur la section intérieure pour reconstruire le champ magnétique sortant.

• équivalent au problème initial pour toute fréquence ;
• valable quelle que soit l’épaisseur du recouvrement ;
• applicable à des guides anisotropes quelconques ;
• prometteur pour des solveurs itératifs efficaces.

Ce que montrent les simulations

Toutes les simulations ont été réalisées avec la bibliothèque XLiFE++ . Elles valident la reconstruction des champs dans des guides non perturbés, puis dans des géométries contenant un obstacle ou une variation de matériau. Dans le cas homogène, l’erreur liée à la condition transparente décroît exponentiellement avec l’épaisseur du recouvrement, jusqu’à atteindre l’erreur propre au maillage.

Les visuels ci-dessous correspondent à deux expériences numériques distinctes. Le schéma et les courbes montrent d’abord comment un mode inverse est identifié dans un guide hétérogène. L’animation tridimensionnelle présente ensuite le calcul d’un mode TEM dans une géométrie multiplement connexe. Ces calculs de modes constituent les données d’entrée des conditions transparentes.

1. La section du guide. Il s’agit d’une coupe perpendiculaire à l’axe de propagation. Les couronnes concentriques alternent deux matériaux de permittivités ε = 10 et ε = 0,5, tandis que μ = 1. Ce fort contraste en « peau d’oignon » favorise l’interaction entre les modes. 2. Les courbes de dispersion. Chaque branche relie la fréquence ω à la constante de propagation β d’un mode. La branche orange est décroissante pour β > 0 : sa vitesse de groupe dω/dβ est négative alors que sa vitesse de phase ω/β est positive. Elle décrit donc un mode inverse.

→

Un résultat marquant

Pour une configuration hétérogène possédant un mode inverse, la condition CtM reconstruit correctement la solution sortante, alors qu’une PML classique sélectionne le mauvais comportement et fournit une solution erronée.

Un autre calcul modal : le guide à deux trous

Dans ce second exemple, le matériau est homogène (\(\varepsilon = \mu = 1\)) mais la section circulaire est percée de deux trous. Elle n’est donc pas simplement connexe : sa topologie autorise deux modes TEM indépendants, présents même à basse fréquence.

L’animation montre l’un de ces modes reconstruit le long du guide. La paroi extérieure est rendue transparente pour laisser voir le champ autour des deux trous. L’échelle colorée code la valeur instantanée, signée, de la composante du champ représentée : bleu pour les valeurs négatives, vert au voisinage de zéro et jaune-orangé pour les valeurs positives. Au cours du temps, le profil transverse se propage suivant l’axe du guide. Ce test vérifie le calcul numérique des modes dans une géométrie complexe ; il est distinct de la configuration à mode inverse présentée juste au-dessus.

Propagation d’un des deux modes TEM du guide homogène à deux trous. Les couleurs suivent l’oscillation du champ tandis que sa structure transverse reste liée à la géométrie de la section.

Perspectives

Ces travaux ouvrent plusieurs prolongements :

• obtenir un critère théorique prédisant l’existence des modes inverses dans le problème complet ;

• achever l’analyse de convergence de la discrétisation par éléments finis en présence des conditions transparentes ;

• exploiter la structure du CtM dans des solveurs itératifs pour traiter des configurations plus réalistes ;

• concevoir des guides présentant des propriétés de diffraction particulières : réflexion nulle, conversion modale ou formes d’invisibilité.

Pour aller plus loin