Détecteurs électromagnétiques basés sur la méthode FDTD
1. Introduction
Dans une simulation FDTD, les détecteurs agissent comme des sondes de mesure virtuelles placées à l’intérieur du domaine de calcul. Ils sont utilisés pour les simulations optiques et radiofréquences (RF) : les mesures sous-jacentes des champs électriques et magnétiques sont les mêmes, même si les longueurs d’onde et les applications peuvent être très différentes. Une analogie utile est celle du voltmètre dans ??. Alors qu’un voltmètre mesure le potentiel électrique, un détecteur FDTD échantillonne les champs électromagnétiques lorsque des ondes traversent une région sélectionnée de la simulation.
Les détecteurs échantillonnent les valeurs du champ en des points déterminés par le maillage FDTD. Un détecteur peut être utilisé pour suivre l’évolution d’un signal dans le temps ou, dans une simulation tridimensionnelle, pour mesurer la distribution de l’énergie électromagnétique sur une surface. Dans OghmaNano, les surfaces des détecteurs apparaissent sous la forme de plans en forme de grille dans la vue 3D, comme illustré dans ??. En fonction de l’analyse, le champ mesuré peut être utilisé pour calculer la transmission, la réflexion ou la puissance diffusée par un objet.
Les champs variant dans le temps sont accessibles au détecteur pendant l’exécution de la simulation FDTD. OghmaNano peut également accumuler directement des champs complexes dans le domaine fréquentiel, ce qui permet d’effectuer des mesures spectrales aux fréquences choisies sans avoir à enregistrer l’intégralité de l’historique temporel et à appliquer une transformée de Fourier par la suite. Cela rend le même système de détection adapté aux composants optiques, aux antennes et aux problèmes de diffusion radiofréquence.
Pour obtenir des informations générales sur le positionnement et la configuration des détecteurs, veuillez consulter le document suivant. détecteurs optiques Section. Cette page est consacrée aux mesures FDTD et à la manière dont ces mesures se traduisent en résultats spectraux et de puissance.
2. Configuration du détecteur
Ouvrez l’éditeur de détecteurs à partir du ruban « Optique », comme indiqué dans ??, ou faites un clic droit sur un plan de détecteur dans la vue 3D et choisissez « Modifier l’objet », comme indiqué dans ??. Une simulation peut utiliser plusieurs détecteurs : par exemple, un détecteur placé avant un dispositif et un autre après, ou un détecteur orienté vers un objet qui diffuse les ondes radio.
L’éditeur de détecteur est présenté dans ??. Son onglet « Configuration » permet de contrôler la géométrie, la position et les options du détecteur, telles que les spectres normalisés. La normalisation d’une mesure de sortie par rapport à une mesure de référence ou d’entrée est particulièrement utile pour les études de transmission, car le spectre de la source peut ne pas être uniforme. Pour les mesures de la puissance traversant une surface de détecteur, l’orientation de cette surface est également importante : le signe du résultat est défini par rapport à sa direction normale.
L’onglet « Maillage de fréquences » dans ?? indique les fréquences auxquelles les informations complexes du champ sont accumulées. Dans les applications optiques, ces points sont souvent exprimés en longueurs d’onde ; pour les simulations RF, ils sont généralement spécifiés en MHz ou en GHz. Cet échantillonnage de fréquences est un paramètre du détecteur et ne doit pas être confondu avec le maillage spatial FDTD utilisé pour résoudre la géométrie et les champs.
3. Accumulateurs dans le domaine fréquentiel de la méthode FDTD.
La méthode FDTD fait évoluer les champs électromagnétiques au fil de différentes étapes temporelles, mais de nombreuses caractéristiques utiles des dispositifs sont mieux décrites en fonction de la fréquence. Au lieu de stocker chaque étape temporelle et de la traiter ultérieurement, un détecteur peut construire la réponse aux fréquences sélectionnées au fur et à mesure que la simulation se déroule. Cela est réalisé grâce à des accumulateurs de Fourier complexes.
Pour un composant de champ échantillonné aux instants \(t_n=n\Delta t\), la quantité accumulée à la fréquence angulaire \(\omega=2\pi f\) est :
\[\widetilde{E}(\omega)=\sum_{n=0}^{N-1} E(t_n)\,e^{-i\omega t_n}\,\Delta t\]
À chaque étape temporelle, une contribution complexe est ajoutée à la fréquence sélectionnée. Une composante oscillant à cette fréquence s’accumule de manière cohérente, tandis que les contributions aux autres fréquences tendent à s’annuler sur une période d’observation suffisamment longue. Le symbole « tilde » indique une intégrale de Fourier, et non encore l’amplitude normalisée d’un champ sinusoïdal. Le même calcul peut être effectué pour les trois composantes du champ électrique \((E_x,E_y,E_z)\) et du champ magnétique \((H_x,H_y,H_z)\).
L’accumulateur conserve à la fois l’amplitude et la phase. Cette distinction est importante : les spectres de transmission peuvent être basés sur les intensités de champ, tandis que le calcul de la puissance électromagnétique directionnelle nécessite la relation de phase entre les champs électrique et magnétique. Pour une réponse continue en onde stationnaire accumulée sur une durée \(T\), l’amplitude complexe du champ maximal est approximativement \(2\widetilde{E}/T\) ; l’expression analogue s’applique à \(\widetilde{H}\). Les temps d’observation finis et les phénomènes transitoires peuvent affecter cette approximation, il est donc important de choisir l’intervalle d’accumulation de manière appropriée en fonction de l’excitation.
Un accumulateur distinct est maintenu pour chaque fréquence demandée, ce qui permet d’obtenir des informations spectrales sans avoir à stocker un grand ensemble de données dans le domaine temporel. Les points de fréquence sont configurés dans Grille de fréquences tab (??). Pour en savoir plus sur la méthode numérique elle-même, consultez Théorie FDTD.
4. Signaux de sortie du détecteur
Une fois la simulation terminée, les résultats du détecteur peuvent être consultés directement à partir de l’onglet « Résultats » d’OghmaNano. Les fichiers disponibles dépendent des mesures du détecteur qui ont été activées et comprennent les mesures dans le domaine temporel, les spectres dans le domaine fréquentiel et, pour les simulations d’éparpillement calibrées, les données sur la puissance diffusée.
Une fois la simulation terminée, les données du détecteur s’affichent sous la forme d’une petite icône de type CCD dans l’onglet des résultats (??). En double-cliquant sur cette icône, les fichiers de résultats présentés dans ?? sont affichés.
Les fichiers lam_E.csv et lam_E_norm.csv sont générés à partir des accumulateurs dans le domaine fréquentiel. Le premier présente un spectre brut lié au champ électrique, tandis que le second présente un spectre normalisé, tel qu’une mesure de sortie référencée à une mesure d’entrée. Ces fichiers sont utiles pour étudier les dispositifs de transmission ou les dispositifs à sélection de longueur d’onde. Les exemples correspondants sont présentés dans ?? et ??.
Le fichier power.csv contient un signal de détecteur dans le domaine temporel, illustré dans ??. Bien que les résultats obtenus dans les domaines temporel et fréquentiel décrivent la même simulation électromagnétique, les spectres du domaine fréquentiel ne sont pas obtenus en appliquant une FFT à la courbe de puissance représentée. Au lieu de cela, ils sont accumulés directement à partir des champs électriques et magnétiques pendant la simulation, comme expliqué dans la section 3.
Dans une simulation de diffusion des ondes radio, il peut également y avoir une scattered_power.csv en sortie. Il ne s’agit pas simplement de la différence entre les puissances totales mesurées avec et sans un objet. Elle est calculée à l’aide des champs électriques et magnétiques diffusés complexes, obtenus en comparant l’exécution principale à une exécution de référence servant à l’étalonnage. La méthode de mesure est expliquée ci-dessous.
5. Mesure de la puissance électromagnétique diffusée.
Lorsqu’une onde électromagnétique rencontre une balle de golf, un objet métallique ou une autre structure tridimensionnelle, le champ résultant est la combinaison de l’onde incidente d’origine et de l’onde diffusée par l’objet. Pour mesurer la diffusion, OghmaNano utilise une simulation d’étalonnage dans laquelle l’objet est absent, tandis que la source, le domaine de calcul, les positions du détecteur, l’échantillonnage de fréquence et le temps d’accumulation restent les mêmes. Une deuxième simulation est ensuite effectuée avec l’objet présent. Cette comparaison permet de séparer le champ électromagnétique diffusé du champ incident enregistré par le détecteur.
À chaque position et fréquence du détecteur, les composantes du champ calibrées sont soustraites de celles enregistrées lorsque l’objet est présent :
\[\widetilde{\mathbf E}_{\mathrm{sc}}=\widetilde{\mathbf E}_{\mathrm{object}}-\widetilde{\mathbf E}_{\mathrm{cal}}\qquad \widetilde{\mathbf H}_{\mathrm{sc}}=\widetilde{\mathbf H}_{\mathrm{object}}-\widetilde{\mathbf H}_{\mathrm{cal}}\]
Ici, \(\widetilde{\mathbf E}\) et \(\widetilde{\mathbf H}\) sont les accumulateurs de Fourier complexes décrits dans la section 3. En soustrayant les champs complexes plutôt que leurs amplitudes, on préserve leurs phases, ce qui est nécessaire pour isoler l’onde diffusée. Les six composantes du champ cartésien sont toutes prises en compte dans le calcul.
Le flux d’énergie électromagnétique est décrit par le vecteur de Poynting. Pour les amplitudes de champ sinusoïdal complexes, sa valeur moyenne dans le temps est :
\[\langle\mathbf S_{\mathrm{sc}}\rangle=\frac{1}{2}\operatorname{Re}\!\left(\mathbf E_{\mathrm{sc}}\times\mathbf H_{\mathrm{sc}}^{*}\right)\]
Le produit vectoriel prend en compte la direction du flux d’énergie, tandis que le conjugué complexe \((*)\) prend en compte la phase relative des champs électrique et magnétique. Le vecteur de Poynting local est mesuré en W/m². Pour les champs sinusoïdaux stables, l’utilisation directe des accumulateurs de Fourier et l’intégration sur la surface du détecteur donnent
\[\boxed{P_{\mathrm{sc}}(f)\simeq\frac{2}{T^2}\sum_{k} \operatorname{Re}\!\left(\widetilde{\mathbf E}_{\mathrm{sc},k}\times \widetilde{\mathbf H}_{\mathrm{sc},k}^{*}\right) \cdot\hat{\mathbf n}\,\Delta A_k}\]
Dans cette expression, \(T\) représente la durée d’accumulation du champ, \(k\) désigne un point d’échantillonnage du détecteur, \(\Delta A_k\) est la zone représentée par ce point, et \(\hat{\mathbf n}\) est la normale à la surface du détecteur. Le facteur \(2/T^2\) est obtenu en convertissant les deux champs accumulés en amplitudes sinusoïdales, puis en calculant leur produit moyen dans le temps. La somme de la composante normale du vecteur de Poynting sur la surface du détecteur donne le spectre de la puissance diffusée, en watts, qui est enregistré sous la forme de scattered_power.csv.
💡 Interprétation du signal : la puissance diffusée est un flux signé mesuré par le détecteur, et non une mesure non négative de l’intensité du champ. Une valeur négative indique que la puissance diffusée nette circule dans la direction opposée à la normale du détecteur choisi ; pour un détecteur orienté vers un objet, cela peut indiquer une diffusion inverse vers la source. La valeur indiquée n’est pas non plus la même que celle obtenue en soustrayant la puissance mesurée dans l’air seul de la puissance totale : ces deux puissances incluent des termes d’interférence entre les ondes incidentes et diffusées, qui sont absents lorsque les champs sont d’abord soustraits.
Cette méthode d’étalonnage est particulièrement utile pour comparer la dépendance en fréquence de la diffusion de différentes formes, tout en maintenant les paramètres de la simulation environnante inchangés. Pour des exemples concrets, voir : diffusion des ondes radio par une balle de golf à une fréquence de 10 GHz et diffusion des ondes radio par le site de StonehengePour obtenir des résultats quantitatifs fiables, utilisez une résolution spatiale et temporelle suffisante, des paramètres d’étalonnage cohérents et placez le détecteur de manière appropriée.