Etude de la diffraction E.M par la méthode itérative basée sur le ...

SETIT2005
• a= 0.1?
Fig.8: convergences du module de courant pour les
deux paramètres ( z0 et z1), pour a= 0.1*?,pixel
sur le domaine du métal
Fig. 9: convergences du module de courant pour
a= 0.1?
D´après les deux figures précédentes, on constate :
• sur le domaine de la source la différence entre les
vitesses de convergence pour les paramètres z1 et
z 0 n´est pas assez important (V (z1)=
3
1 V (z0))
• sur le domaine du métal la figure 8 explique bien
le choix de z1 par rapport à z0, en effet pour
avoir la convergence, il faut 250 itérations pour
z0, mais pour z1 30 itérations sont suffisantes.
3-2-2 Variation de la densité de courant :
On discrétise le contour circulaire au nombre M=
64 pixels, on localise le domaine de la source sur P
pixels et le domaine du métal sur (M- P) pixels.
Ainsi, on a bien définit les fonctions indicatrices
H M et HS.
Pour calculer la densité totale du courant. Nous
allons effectuer les premières itérations a fin de
dégager la forme générale du courant à l’itération
N.
La densité total de courant est donnée par :
Fig. 10: variation du module de courant pour a=
0.1?
Pour a= 0.1?, on remarque que le courant est
maximal sur le domaine de la source, il démunit
sur le métal jusqu’il s’annule sur le« dos » du
cylindre l et on voie qu’il y a une cassure sur le
pic (domaine de la source), celle ci est interprétée
par l’apparition de l’effet de bords sur les limites
des deux domaines (source –métal).
• a= ?
( N)
B = HS
( N)
J
E
2
0
- M
A
( N)
= Γˆ FFT ( B )
z
( N ) ( N )
( )
A − B
N
= [10]
z
−1
( N −1 )
H FF T ( A )
[9]
Les courbes suivantes sont étudiées pour une source
localisée sur 3 pixels de la structure de diffraction.
Les courbes de convergence en fonction de nombre
d’itérations sont calculées au pixel numéro 10 (pris
sur le domaine du métal).
Fig. 11 : convergences du module de courant pour
a= ?
4