Conditions aux limites pour la méthode FDTD

1. Introduction

Dans les simulations par la méthode des différences finies dans le domaine temporel (FDTD), les conditions aux limites définissent la manière dont les champs électromagnétiques se comportent aux limites extérieures du domaine de calcul. En termes physiques, elles déterminent ce qui se passe lorsqu’une onde atteint le bord de la cellule de simulation. Il s’agit d’un des aspects les plus importants d’un modèle FDTD, car le domaine numérique est toujours fini, alors que le problème physique représente souvent une propagation dans l’espace ouvert, des structures répétitives ou des cavités confinées. Le rôle de la condition aux limites est donc d’imposer la continuation mathématique correcte des champs au-delà du maillage tronqué.

Si une condition aux limites est mal choisie, les ondes se réfléchissent et reviennent dans la zone de calcul, ce qui perturbe la solution physique et peut entraîner la formation d’artefacts importants liés aux ondes stationnaires. Dans le pire des cas, l’énergie est piégée à l’intérieur du volume de simulation et la distribution du champ calculée ne représente plus le dispositif prévu. C’est particulièrement important dans les structures photoniques ouvertes, les problèmes de diffusion, les systèmes rayonnants et les calculs de réponse impulsionnelle, où les ondes sortantes doivent quitter le domaine de manière nette. À l’inverse, dans les structures périodiques, le comportement correct n’est pas l’absorption, mais la répétition, et dans certains problèmes de test, il est souhaitable de fixer explicitement les champs à la limite. Le choix de la condition aux limites doit donc correspondre à la physique du modèle.

D’un point de vue mathématique, les équations de Maxwell, exprimées sous forme de rotationnel, sont résolues sur une grille de Yee finie. Pour une région isotrope sans source, ces équations peuvent être formulées comme suit :

\[ \mu \frac{\partial \mathbf{H}}{\partial t} = - \nabla \times \mathbf{E}, \qquad \varepsilon \frac{\partial \mathbf{E}}{\partial t} + \sigma \mathbf{E} = \nabla \times \mathbf{H}. \]

Les équations de mise à jour des valeurs intérieures découlent directement de ces relations, mais aux cellules les plus externes, le modèle dépasse la zone de maillage définie. La condition aux limites fournit les informations manquantes nécessaires pour compléter la mise à jour. En pratique, cela signifie que chaque face externe de la boîte de simulation doit être associée à une règle numérique qui détermine si les champs tangents sont fixés, absorbés, répétés périodiquement ou atténués par le biais d’une couche d’absorption adaptée.

Dans OghmaNano, les conditions aux limites sont définies à partir de Conditions aux limites appuyer sur le bouton Optique ruban (voir) Figure ??). L’éditeur correspondant est indiqué ci-dessous. Figure ??) permet d’attribuer une condition différente à chacune des faces externes de la région FDTD. Les types de conditions aux limites disponibles sont les suivants : Dirichlet, Mur ABC, Périodique.et PML.

Les six limites sélectionnables correspondent aux six faces de la boîte de simulation cartésienne. Dans la configuration d’affichage par défaut utilisée par OghmaNano, \(y_0\) représente la face supérieure de la zone de simulation. \(y_1\) c’est la face inférieure. \(x_0\) c’est la face gauche, \(x_1\) est-ce le bon visage ? \(z_0\) est la face la plus proche de l’utilisateur, et \(z_1\) est la face la plus éloignée de l’utilisateur. La représentation schématique de l’orientation est illustrée dans Figure ??.

L’éditeur des conditions aux limites attribue une règle à chacune de ces six faces, de manière indépendante. Cela permet de créer des domaines mixtes, par exemple avec des limites périodiques dans les directions latérales et des limites absorbantes dans la direction de propagation, ou avec une couche PML sur toutes les faces pour la propagation dans l’espace libre. Lorsqu’une couche PML est sélectionnée, un paramètre de profondeur supplémentaire est disponible pour cette face. Cette profondeur est exprimée en cellules du maillage et détermine le nombre de couches FDTD externes réservées au milieu absorbant.

2. Éditeur de conditions aux limites

L’éditeur des conditions aux limites propose une fonction de contrôle distincte pour les faces. \(y_{\min}\), \(y_{\max}\), \(x_{\min}\), \(x_{\max}\), \(z_{\min}\)et \(z_{\max}\), correspondant respectivement à \(y_0\), \(y_1\), \(x_0\), \(x_1\), \(z_0\)et \(z_1\) dans le sens indiqué par la convention de représentation. Figure ??Pour chaque face, l’option sélectionnée détermine la manière dont les champs électriques et magnétiques tangents sont mis à jour au niveau du bord extérieur de la grille.

3. Condition aux limites de Dirichlet.

Une condition aux limites de Dirichlet fixe la valeur du champ sur la frontière à une valeur prédéfinie. Dans l’utilisation la plus courante de la méthode FDTD, cette valeur est égale à zéro, ce qui signifie que la frontière impose...

\[ \mathbf{E}_{\mathrm{tan}} = 0 \qquad \text{on the selected face,} \]

Par conséquent, la limite se comporte comme une paroi parfaitement réfléchissante pour le champ électrique, et toute onde sortante est renvoyée vers la zone de simulation.

3. Condition aux limites de la paroi ABC.

Les conditions aux limites absorbantes sont conçues pour réduire la réflexion en simulant le comportement d’une onde se propageant vers l’extérieur au niveau de la limite du domaine de calcul. L’idée est de satisfaire une équation d’onde unidirectionnelle normale à la limite, de sorte que l’énergie quittant la grille continue de se propager vers l’extérieur plutôt que de se réfléchir.

Pour une onde unidimensionnelle se propageant dans la direction \(+x\), la relation exacte décrivant l’onde sortante est la suivante :

\[ \frac{\partial u}{\partial x} + \frac{1}{c}\frac{\partial u}{\partial t}=0. \]

La condition limite absorbante du premier ordre de Mur est obtenue en discrétisant cette équation sur la grille FDTD. Sur la limite gauche, une forme standard est utilisée.

\[ u_0^{n+1} = u_1^{n} + \frac{c\Delta t-\Delta x}{c\Delta t+\Delta x} \left( u_1^{n+1}-u_0^{n} \right), \]

une expression analogue se retrouve sur le côté opposé. Ici, \(u_i^n\) désigne le champ à l’indice spatial \(i\) et à l’indice temporel \(n\), \(\Delta x\) représente l’espacement du maillage, et \(\Delta t\) est le pas de temps de la méthode FDTD.

Les conditions aux limites de type « mur absorbant » sont peu coûteuses en termes de calcul et nécessitent très peu de mémoire supplémentaire, ce qui les rend intéressantes pour les simulations légères ou les calculs exploratoires. Elles visent à approximer une onde sortante de manière à ce que l’énergie puisse quitter le domaine de simulation plutôt que de se réfléchir sur la limite extérieure.

Étant donné que la condition de Mur n’est qu’une approximation, une certaine réflexion persiste toujours. La précision diminue pour les angles d’incidence obliques, les impulsions à large bande et les champs de diffusion multidimensionnels complexes. Les conditions aux limites de Mur sont donc les plus appropriées pour les problèmes de propagation simples, où les ondes atteignent la limite approximativement perpendiculairement à la surface et où une très grande précision n’est pas requise.

Pour les simulations optiques complexes, les structures rayonnantes ou les cas où les ondes atteignent la limite sous plusieurs angles, il est préférable d’utiliser des couches parfaitement adaptées (PML). Les limites des PML introduisent une région absorbante qui élimine les ondes sortantes de manière beaucoup plus efficace et avec une réflexion nettement inférieure à celle des conditions de Mur.

4. Condition aux limites périodique

Une condition aux limites périodique permet d’identifier une face de la région de simulation avec la face opposée, de sorte que le champ qui sort d’un côté réentre par l’autre. Cette approche est appropriée lorsque la structure physique se répète indéfiniment dans l’espace et que la cellule de calcul représente une seule unité de cette géométrie répétitive.

Sous sa forme la plus simple, la périodicité impose...

\[ \mathbf{E}(x+L_x,y,z,t)=\mathbf{E}(x,y,z,t), \qquad \mathbf{H}(x+L_x,y,z,t)=\mathbf{H}(x,y,z,t), \]

De même, dans les directions \(y\) ou \(z\), lorsque ces limites sont définies comme périodiques. Sur le plan numérique, les valeurs du champ requises juste au-delà d’un côté du domaine sont prélevées du côté opposé. Par conséquent, les faces opposées doivent être traitées comme une paire appariée. Sur le plan physique, cela signifie que la solution correspond non pas à un objet isolé, mais à un réseau infini de copies de la cellule simulée.

Les conditions aux limites périodiques sont fréquemment utilisées pour les structures qui se répètent dans l’espace, comme les réseaux de diffraction, les cristaux photoniques, les cellules élémentaires de métamatériaux et les géométries périodiques de guides d’ondes. En modélisant uniquement une seule cellule élémentaire, on peut réduire considérablement les coûts de calcul tout en représentant le comportement de la structure répétée dans son ensemble, à condition que la source et la géométrie soient cohérentes avec la périodicité supposée.

Dans le cas de frontières périodiques, les champs qui sortent d’un côté du domaine de simulation réapparaissent de l’autre côté. L’énergie n’est donc pas absorbée, mais elle est recyclée dans la cellule de calcul. Par conséquent, les frontières périodiques ne doivent pas être utilisées pour représenter un espace ouvert. Si des frontières périodiques sont appliquées sur toutes les faces de la région de simulation, aucune énergie ne peut quitter le domaine, et l’énergie électromagnétique totale peut s’accumuler au fil du temps. Dans de nombreuses simulations FDTD pratiques, des frontières périodiques sont donc appliquées dans les directions transversales, tandis que des frontières absorbantes, telles que les PML, sont utilisées dans la direction de propagation.

Dans le cas le plus général, les champs périodiques peuvent comporter un déphasage entre les cellules adjacentes. Ce phénomène est décrit par la condition de Bloch.

\[ \mathbf{E}(x+L_x,y,z,t)=\mathbf{E}(x,y,z,t)e^{ik_xL_x}, \]

il existe une relation analogue pour le champ magnétique. La condition aux limites périodique standard mise en œuvre ici correspond au cas particulier \(k_x=0\), ce qui signifie que les champs se répètent de manière identique sur les faces opposées du domaine de simulation.

5. Condition aux limites de PML

La couche parfaitement adaptée (CPA) est la méthode standard utilisée pour créer une frontière absorbante performante dans les simulations FDTD modernes. Son objectif est d’éliminer les ondes électromagnétiques sortantes en minimisant la réflexion sur une large plage de fréquences et pour une grande variété d’angles. Au lieu d’imposer une formule de frontière unique au bord extérieur, la CPA introduit un milieu absorbant artificiel qui entoure la région de simulation physique. Les ondes qui pénètrent dans cette couche subissent une atténuation sans rencontrer de discontinuité d’impédance marquée à l’interface.

L’idée maîtresse est que la couche PML est adaptée au milieu intérieur au niveau du plan d’entrée, de sorte qu’une onde incidente ne rencontre pas de discontinuité qui générerait elle-même une réflexion. Une fois à l’intérieur de la couche, les champs sont atténués de manière exponentielle. Sous forme continue, cela peut être interprété par une dilatation complexe des coordonnées, par exemple.

\[ x \;\rightarrow\; \int_0^x s_x(\xi)\,d\xi, \qquad s_x(\xi)=\kappa_x(\xi)+\frac{\sigma_x(\xi)}{j\omega\varepsilon_0}, \]

avec des facteurs d’étirement correspondants dans les autres directions. Ici, \(\sigma_x\) représente un profil de conductivité artificielle et \(\kappa_x\) est un facteur d’échelle utilisé dans les formulations de PML à coordonnées étirées. L’effet est que les champs qui se propagent sont atténués lorsqu’ils traversent la couche, tout en restant approximativement sans réflexion à l’interface avec le domaine physique.

Dans OghmaNano, la sélection de la méthode PML pour une face donnée active cette couche absorbante sur ce côté de la zone de simulation. Le paramètre associé, qui définit la profondeur de la couche PML, détermine le nombre de cellules du maillage allouées à cette couche. Une couche PML plus épaisse assure généralement une atténuation plus importante, car le champ a plus d’espace pour s’affaiblir. Si la couche PML est trop fine, une partie de l’énergie peut atteindre la limite extérieure et être réfléchie. Si elle est suffisamment épaisse, le champ sortant est atténué jusqu’à une valeur négligeable avant d’atteindre le bord extérieur.

Une façon pratique de comprendre la dépendance à l’épaisseur consiste à noter qu’une composante du champ à l’intérieur d’un milieu absorbant diminue approximativement selon la loi suivante :

\[ E(d)\sim E(0)\exp(-\alpha d), \]

où \(d\) représente la distance parcourue dans la couche PML et \(\alpha\) est un coefficient d’atténuation effectif déterminé par le gradient de conductivité. Par conséquent, l’augmentation du nombre de cellules PML réduit l’amplitude résiduelle atteignant la limite extérieure. En termes pratiques, la profondeur de la couche PML doit être suffisamment importante pour que les réflexions provenant de l’arrière de la couche soient négligeables par rapport au signal physique d’intérêt.

La méthode PML est généralement la condition aux limites privilégiée pour les simulations ouvertes en optique photonique et en électromagnétisme, notamment pour la modélisation du rayonnement des guides d’ondes, de la diffusion par des objets isolés, de l’émission de type antenne et de la propagation d’impulsions transitoires. Elle est beaucoup plus robuste que les simples conditions aux limites absorbantes lorsque les champs atteignent la limite de manière oblique ou contiennent un large spectre de fréquences. Pour cette raison, la méthode PML est le choix par défaut dans de nombreux modèles FDTD où l’on s’attend à ce que les ondes quittent le domaine de calcul.

Même en présence d’une couche PML, il est essentiel de respecter les bonnes pratiques de modélisation. La couche absorbante ne doit pas être placée trop près d’une région où le champ est fortement décroissant ou très résonant, car les champs proches non propagatifs peuvent interagir avec le matériau absorbant et modifier la solution. En général, la structure physique doit être séparée de la couche PML par un espace libre suffisant pour permettre à l’onde sortante de se former de manière nette. Utilisée de cette manière, la couche PML offre l’approximation la plus précise, dans le cadre de la méthode FDTD, d’un domaine extérieur ouvert et non réfléchissant.